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Introduction 


Nous  considerons  le  probleme  de  filtrage  non  lineaire  dans  une  situation  pjurticu- 
liere,  celle  de  I’observation  faiblement  bruitee.  ,En  pratique,  plusieurs  problem, es  de 
filtrage,  en  particulier  en  genie  aerospatial,  peuvent  etre  modelises  de  fagon  conven- 
able  par  des  situations  de  ce  t3q)e. 

Ce  travail  est  consacre  a  I’etude  de  problemes  asymptotiques.  Nous  considerons 
des  filtres  sous-optimaux,  dependant  d’un  parametre  e  et  nous  etudions  leur  com- 
portement  lorsque  e  converge  vers  zero. 

II  est  bien  connu  que  Ic  probleme  de  filtrage  non  lineaire  est  un  probleme  de 
dimension  infinie,  sa  solution  etant  donnee  par  exemple  par  Pequation  de  Zakai. 
Quelques  tentatives  ont  ete  faites  pour  trouver  des  classes  de  problemes  pour  lesquels 
un  filtre  de  dimension  finie  existe  (cf.  Benfis  [1],  Levine  [19]  et  Wong  [30],  entre 
autres).  Cependant  ces  tentatives  se  sont  montr^  d'un  interev  limite,  n’ayant  pas 
abouti  a  des  classes  de  problemes  “interessants”.  Des  travaux  posterieurs  viendraient 
prouver  qu’on  ne  peut  construire  un  filtre  de  dimension  finie  que  pour  une  classe  de 
problemes  negligeable.  Les  conditions  que  problemes  doivent  verifier  sont  li^s 
a  I’existence  d’une  algebre  de  Lie  de  dimension  finie  (cf.  par  exemple  Chaleyat- 
Michel  [5]).  On  regarde  done  avec  interet  la  possibilite  d’obtenir  de  “bonnes”  ap¬ 
proximations  du  filtre  optimal  qui  puissent  etre  calculees  a  I’aide  d’un  nombre  fini 
de  statistiques  donnees  comme  des  solutions  d’un  systeme  d’equations  differentielles 
stocheistiques  conduites  par  le  processus  d’observation. 

De  nombreux  travaux  ont  ete  tffectuw  dans  ce  domaine  -  voir  par  exemple 
Fleming-Ji-Pardoux  [9],  Fleming-Ji-Salame-Zhang  [10]  et  Roubaud  [28]  pour  le 
probleme  lineaire  par  morceaux,  Fleming-Pardoux  [11]  pour  le  probleme  monotone 
par  morceaux,  Katzur-Bobrovsky-Schuss  [17],  Picard  [22]  et  [23]  et  Bensoussan 
[2]  pour  le  probleme  asymptotiquement  observe,  Yaesh-Bobrovsky-Schuss  [31]  et 
Picard  [26]  pour  le  probleme  asymptotiquement  non  observe,  entre  autres. 

Cette  th^e  peut  etre  vue  comme  composee  de  trois  chapitres  independants. 
Cependant  tous  les  trois. s’inscrivent  dans  un  memesujet:  le  comportement  asymp- 
totique  de  problemes  de  filtrage. 

Dans  un  premier  chapitre  nous  etudions  un  probleme  de  filtrage  non  lineaire 
unidimensionnel  en  temps  discret,  lorsque  le  bruit  d’observation  est  faible.  11  se 
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caracterise  par  une  fonction  d’observation  injective. 

Le  deuxieme  chapitre  est  consacre  a  un  probleme  de  iiltrage  lineaire  pcir  mor- 
ceaux  en  temps  discret.  Nous  nous  interessons  aux  aspects  de  mise  en  eeuvie  d’al- 
gorithmes  pour  resoudre  numeriquement  ce  probleme,  Ce  travail  a  ete  fait  en  col¬ 
laboration  avec  M.C.  Roubaud. 

Dans  le  dernier  chapitre  nous  etudions  un  probleme  en  temps  cpntinu,  ayant  des 
dimensions  differentes  pour  le  signal  et  I’observation.  Plus  exactement  nous  consi- 
derons  un  signal  de  dimension  2  dont  seulement  une  des  composantes  est  bruits  et 
un  processus  d’observation  de  dimension  1  observant  I’autre  composante. 


Filtres  approchds  pour  un  probleme  de  filtrage  non  lineaire  en  temps 
discret  avec  petit  bruit  d’observation 

[1]  Filtres  approches  pour  un  probleme  de  filtrage  non  lindairc  en  temps  discret 
avec  petit  bruit  d’observation,  RR  1142,  INRIA,  Dec.  1989. 

[2]  Filtres  approches  pour  un  probleme  de  filtrage  non  liniaire  en  temps  discret 
avec  petit  b^itit  d’observaiion,  Lecture  N.  in  Cont.  and  Inf.  Sci.,  Proc.  of  the 
9th  Internat.  Conf.  on  Analysis  and  Optim.  of  Systems,  Antibes,  .June  12-15, 
1990,  pp.  198-207. 


On  considere  le  systeme  differentiel  stochastique 


[  Vk 


Xk  b{Xit)  At  +  <T  y/At  Wk+i , 


oil  {lyjfc}  et  {%}  sont  des  bruits  blancs  gaussiens  independants,  e  est  un  parametre 
suppose  “petit”  et  At  depend  dc  t,  Lc  processus  {X*}  est  le  signal  a  estimer  et  {y*} 
est  le  processus  d ’observation  engendrant  la  filtration  yk- 

On  suppose  que  la  foJiction  ft  est  injective,  i.e.  le  signal  est  asymptotiquement 
observe.  On  cherche  des  approximations  Mk  de 

Xk  =  E[Xk\yk] 


qui  soient  donnees  de  fa^on  recursive  et  qui  puissent  etre  calcul^es,  a  I’instant  k,  a 
partir  de  I’observation  Pour  ces  processus  {Mk},  on  veut  pouvoir  estimer  Xk—Mk 
lorsque  e  converge  vers  0. 

Nous  commengons  par  I’etude  du  cas  lineaire  et  nous  y  trouvons  I’inspiration 
necessaire  pour  construire  les  filtres  approches  dans  le  cas  non  lindaire.  Nous  ex- 
ploitons  I’idee  de  Picard  [22]  et  [23]  pour  le  probleme  en  temps  continu  et  nous 
I’adaptons  au  probleme  en  temps  discret.  Nous  faisons  appel  a  des  outils  comme 
les  changements  de  probabilite  et  une  version  discrete  du  Theoreme  de  Girsanov, 


la  derivation  par  rapport  a  la  condition  initiale  et  une  formule  d’integration  par 
parties.  Nous  demontrons  quCj  si  At  =  e®  avec  a  >  1,  le  filtre 

Mk+i  =  Ah  +  HMk)  +  iVk+r  -  h{Mk)) , 

avec  condition  initiale  Mo  =  mo  quelconque,  approche  Xk  a  I’ordre  0(e  V 
si  a  =■  1,  le  filtre  de  dimension  2 

Afk+i  —  Aik  +  i(Mfc)  e  +  0k  ivk+i  —  h{Mk)) 

^  h'iAik)[Sk-t  +  anmk-i)] 
h^AIk)[  dk-i  +  h»{Mk-i)]  +  h'{AIk-i)  ’ 

avec  condition  initiale  AIq  —  toq,  Oq  —  6q  quelconque  telle  que  cBq  >  0 ,  approche 
le  filtre  optimal  a  I’ordre  O^e) ;  si  a  <  1,  iJ  cst  suffisant  d’inverser  la  fonction 
d’observation. 

Les  estimations  obtenues  verifient  la  propriete  d’uniformite  par  rapport  au  temps 
et  montrent  que  ces  filtres  approchcs  out  peu  de  memoire. 

Nous  etudions  ensuite  un  cas  particulier  oil  la  fonction  d’obervation  est  supposee 
non  injective.  II  s’agit  du  cas  d’un  systeme  lin^ire  par  morceaux. 

Filtrage  lineaire  par  morceaux  d’un  systhme  en  temps  discret  avec  petit 
bruit  d’observation 

Nous  considerons  le  systeme  lineaire  par  morceaux  suivant: 

f  Xk+i  =  Xk  +  eb{Xk)  +  y/e(T{Xk)uk 
\  ijk  =  h(Xk)  +  \^Vk, 

oil  {u;^}  et  {uit}  sont  des  bruits  blancs  gaussiens  independapts  et 

6(x)  =  5_xl{x<o>  d'  B+xl{x>o} 

<7(x)  <7_l{ar<0)  +  ^{x>0) 

h{x)  ~  -f  i/+Xl{*>o}. 

Si  H+H-  >  0,  i.e.  la  fonction  h  est  monotone,  le  problem-'  n’uffre  aucune  diffi- 
culte.  Nous  nouf  interessons  au  cas  H^.H-  <  0,  i.e.  h  n’est  pas  monotorie.  Dans  ce 
C2is,  bien  que  h\Xt)  jiuisse  etre  facilement  estime  il  n’eii  est  pas  de  ineme  pour  Xt. 
L’idee  de  base  est  la  suivante:  des  que  le  signal  reste  dans  une  region  ou  la  fonc¬ 
tion  d’obervation  est  injective  (voir  lin wire),  on  peut  appliquer  on  filtre  classique, 
type  Kalman-Bucy.  Le  probleme  se  reduit  done  a  la  recherche  d’une  methode  pour 
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detecter  les  changements  de  signe  du  processus  {X*},  ainsi  que  determiner  le  signe 
pris  dans  chacun  des  intervjilles  ainsi  separes. 

Cette  idee  avait  deja  ete  exploitee  dans  Fleming-Ji-Pardoux  [9],  pour  un  sys- 
teme  lineaire  par  morceaux,  sous  une  certciine  hypoth^e  de  “detectabilite^  c  t,  plus 
tard,  par  Fleming-Pardoux  [11]  pour  un  systeme  non  lineaire  par  morce^•  'X.  Le 
probleme  en  temps  discret  a  ete  aborde  dans  Fleming-Ji-Salame- Zhang  [10].  Nous 
etudions  ici  un  probleme  plus  general  que  celui  traite  pan-  ces  auteurs.  En  efFev,  nous 
considerons  une  hypothec  plus  large  que  I’hypothese  de  “detectabilite”  emise  dans 
I'etude  precedente.  Nous  supposons  que 

’  Hial  ^  Hlal  {HD\) 

i  ou 

Hlal  =  Hlal  et  B+  7^  B.  {HD2) . 

Ceci  permet  de  traiter  en  particulier  le  cas  d’une  fonction  d’observation  symetrique 
par  rapport  a  I’axe  des  ordonnees. 

Les  aspects  de  mise  en  oeuvre  d’algorithmes  pour  resoudre  ce  probleme  sont 
abordes. 

Nous  pr^entons  des  tests  permettant  de  separer  les  intervalles  ou  le  sigTial  est 
positif  ou  negatif  et  de  determiner  ce  signe.  Un  premier  type  de  tests  detecte  les 
changements  de  signe  de  {A'fc}.  Deux  tests  sont  proposes:  un  test  sur  les  observations 
et  un  test  sur  les  sorties  des  filtres  de  Kalman-Bucy.  Le  premier  de  ces  tests  a 
deja  ete  etudie  dans  [10],  sous  I’hypothese  (HDl).  Le  deuxieme  type  de  tests  vise 
la  determination  du  signe  de  {A’^}  dans  I’intervalle  de  monotonie  detecte  par  un 
test  du  premier  type.  Sous  I’hypothese  (HDl),  deux  tests  on  ete  decrits  dans  [10]: 
un  test  de  variation  quadratique  et  un  test  de  rapport  de  vraisemblance  sur  les 
accroissements  des  observations.  Sous  i’hypothese  (HD2),  cependant,  le  premier  de 
ces  tests  ne  permet  pas  de  prendre  une  decision  tet  on  se  contente  done  d’adapter 
le  deu  <ieme.  Nous  decrivons  un  nouveau  test  de  rapport  de  vraisemblance  sur  les 
sorties  des  filtres  de  Kalman-Bucy  (cf.  Roubaud  [28],  pour  le  probleme  en  temps 
continu). 

L’application  de  ces  tests  repose  sur  une  methode  efficace  de  determination  des 
bornes.  Nous  presentons  une  etude  purement  heuristique  permettant  d’obtenir  des 
formules  pour  ces  quantites,  ainsi  que  pour  les  temps  d’attente  esper&  dans  les  tpsts 
de  decision  sur  le  s;gne. 

Nous  comparons  les  differents  tests  et  etudions  I’influence  des  parametres  du 
systeme  sur  les  r&ultats  obtenus  sur  des  exemples.  En  particulier  nous  remarquons 
que  les  temps  d’attente  d’un  test  sous  I’hypothese  (HD2)  sont  beaucoup  plus  longs 
que  ceux  sous  I’hypothese  (HDl). 

Nous  etudions  ensuite  une  autre  situation  ou  la  fonction  d’observation  n’est 
pas  injective.  II  s’agit  d’un  cas  ou  le  signal  est  un  processus  bidimensionnel  et 
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I’observation  est  un  processus  unidimensionnel. 

Fiftres  approch4s  pour  un  probleme  de  filtrage  non  Irnifcaire  d’une  diffu¬ 
sion  de  dimension  2  mesur^e  par  un  processus  de  dimension  1  faiblement 
bruite 

Nous  nous  intertssons  au  probleme  de  filtrage  associe  au  systeme  difFerentiel  stochas- 
tique,  au  sens  d’lto 

-  dx2{t)  =  f2{xi{t),X2it))dt  +  a{xi{t),Xi(t))dw{t), 
dy{t)  =  h{xi{t))  dt  +  e  dibit)  ^ 

Nous  cherchons  done  a  estimer  le  signal  {Xt  =  (a;i(t),i2(0)}>  ^  valeurs  dans  IR^, 
a  I’aide  du  processus  d’observation  yt,  processus  a  valeurs  dans  IR..  Notre  but  est 
d’etudier  le  comportement  asymptotique  d’un  tel  estimateur,  lorsque  s  tend  vers  0. 
Comme  dans  Picard  [26],  nous  privilegions  des  methodes  permettant  de  determiner 
I’ordre  de  grandeur  de  I’erreur  de  filtrage,  voir  sa  vitesse  de  convergence.  Nous  y 
trouvons  I’inspiration  necessaire. 

Nous  traitons  d’abord  le  cas  ou  la  function  k  est  lineaire  et  /j  est  lineaire  par 
rapport  a  X2.  Nous  considerons  un  filtre  approcb'  {A/*}  du  type  propose  dans 
Yaesh-Brobovsky-Schuss  [31]: 

=  /j(mi(t).m2(t))  dt  +  [dy{t)  -  Hmi{t)dt] 

dm2{i)  =  f2{xi{t),X2{t))dt  +  ^[dy{t)  -  Hmi{t)dt], 

avec  condition  initiaie  Mo  et  ou  a  est  une  constante  strictement  positive.  Nous 
d.  '^ortrons  que  le  moment  d’ordre  2p  de  Xk  —  Mk  est  asymptotiquement  d’ordre 
pour  la  premiere  composante,  et  d’ordre  pour  la  .seconde.  Nous 

den  itrons  que  la  condition  initiaie  disparait  exponcntiellement  vite. 

Nous  traitons  ensuite  le  probleme  general.  Par  une  transformation  de  coordon- 
nees  nous  nous  ramenons  au  cas  etudie  precedemment.  Nous  obienons  les  memes 
resultats  pour  la  covariance  de  Xy,  —  Mk- 
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Filtres  approches  pour  un 
probleme  de  filtrage  non  lineaire 
en  temps  discret  avec  petit  bruit 
d’observation 
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FILTRES  APPROCHES  POUR  UU  PROBLEME  DE 

FILTRAGE  NON  LINBAIRE  DISCRBT  AVEC  PETIT 
BRUIT  D’OBCERVATION 

APROXIMATE  FILTERS  FOR  A  NONLINEAR  DISCRETE 
TIME  FILTERING  PROBLEM  WITH  SMALL 
OBSERVATION  NOISE 


Resume 

On  fait  une  etude  asymptotique  d’un  probleme  de  filtrage  non  lineaire  en  temps  discret 
unidimensionnel,  lorsqu’un  certain  parametre  e  tend  vers  0.  On  traite  le  cas  d’un  pro’oleme 
de  filtrage  non  lineaire  en  temps  discret  issu  d’un  probleme  en  temps  continu  avec  petit 
bruit  d’obervation.  Des  filtres  approchds  de  dimension  finie  sont  proposes  et  les  estimations 
sur  les  performances  de  ces  filtres  sont  enoncees  et  demontrecs.  On  commence  par  faire 
une  etude  dans  le  cas  lineaire,  ce  qui  vanous  inspirer  les  filtres  approchfe  h.  prendre  dans  le 
cas  general.  Pour  le  resultat  concernant  I’erreur  entre  le  filtre  approche  et  le  filtre  optimal, 
la  demonstration  presentee  utilisera  des  changements  de  probability  et  la  ddrivation  par 
rapport  a  la  condition  initiale.  Pour  terminer,  on  presente  des  r^sultats  d’application 
de  ces  filtres  a  un  exemple  et  on  constate  que  les  propositions  enoncees  sont  verifiyec 
numeriquement. 


Abstract 

We  study  the  asymptotic  behaviour  of  a  nonlinear  one-dimensional  filtering  discrete 
time  problem,  as  some  parameter  e  tends  to  0.  We  treate  the  case  of  a  nonlinear  discrete 
time  problem  coming  from  a  continuous  time  one  with  small  observation  noise.  Finite 
dimensional  aproximate  filters  are  proposed  and  results  concerning  estimations  of  their 
performance  are  stated  and  proved.  We  present,  at  first,  a  brief  study  of  the  linear  case 
and  get  the  necessary  inspiration  to  construct  the  approximate  filters  in  the  general  case. 
The  proof  of  the  result  we  give  concerning  the  error  between  the  approximate  filter  and 
the  optimal  one  will  make  use  of  probability  changes  and  differentiation  with  respect  to 
the  initial  condition.  Finally,  we  present  the  results  obtained  when  applying  those  filter.*'  to 
an  example  and  we  notice  that  the  propositions  previously  stated  are  verified  numerically. 


1  Introduction 


On  considere  le  probleme  de  filtrage  non  lineaire  en  temps  continu  suivant : 

On  a  un  signal  X,  solution  de  I’equation  differentielle  stochastique  unidi- 
mensionnelle 

dXt  =  b  (Xt)  dt  +  c  dw\  ,Xo  =  ^  (1) 

et  on  dispose  de  I’observation  unidimensionnelle  verifiant 

dYt  =  h{Xt)  dt  +  edwf  ,Yo  =  0  (2) 

oil  sont  des  processus  de  Wiener  standards  independants  et  ^  est  une 

variable  al^toire  independante  de  et  w^.  Le  parametre  e  est  suppose  petit. 

Le  probleme  de  filtrage  consiste  a  calculer,  pour  toute  fonction  mesurable 
d>,  I’esperance  conditionnelle  de  ^(Xt),  etant  donnee  I’observation  Y,  jusqu’a 
I’instant  i. 

On  aimerait  r«oudre  numeriquement  ce  probleme.  Une  possibilite  est  donn^ 
par  la  resolution  de  I’equation  de  Zakai,  une  fois  discretis^,  ce  qui  a  ete  etudie 
dans  [Le  Gland].  Cependant  un  tel  calcul  presente  une  grande  complexite,  puis- 
qu’il  s’agit,  dans  le  cas  general,  d’une  ^nation  de  dimension  iniinie.  Une  autre 
possibilite  consiste  a  utiliser  les  filtres  approcha  propose  dans  [Picard]  et  les 
discretiser  en  temps.  Ces  filtres  ctwit  de  dimension  finie,  le  probleme  devient 
plus  simple  a  r«oudre.  De  plus,  on  y  trouve  des  resultats  concemant  I’erreur  de 
filtrage  et  en  particulier  ceux  qui  correspondent  au  filtre  de  Kalman  etendu.  La 
methode  de  demonstration  utilise  fait  appel  a  des  changements  de  probabilite, 
a  la  theorie  de  retournement  du  temps  pour  les  diffusions  et  a  des  resultats 
sur  les  flots  definis  par  des  equations  differentielles  stochastiques.  Dans  [KBS] 
on  peut  trouver  les  premiers  resultats  sur  ce  genre  de  filtres,  cependant  leur 
justification  reste  purement  formelle.  Notons  ^alement  que  dans  [Bensoussan], 
I’auteur  considere  le  meme  probleme  de  filtrage  en  temps  continu  et  utilise  une 
methode  autre  que  celle  developpee  dans  [Picard]  pour  demontrer  les  r^ultats 
asymptotiques. 

En  pratique,  neanmoins,  on  ne  dispose  souvent  que  de  ^observation  a  des 
instants  discrets  et,  de  plus,  pour  mener  des  calculs  sur  ordinateur,  il  est  n&es- 
saire  de  manipuler  des  processus  discrets.  Au  lieu  de  discretiser  le  filtre  on  peut 
penser  a  discretiser  le  probleme  lui  meme.  L’avantage  de  cette  approche  reside 
dans  le  fait  que,  d^  qu’on  se  limite  a  des  processus  en  temps  discret,  on  n’a 
plus  besoin  de  la  th^rie  des  flots  ou  du  retournement  du  temps  et  meme  pas 
de  I’integrale  d’lto. 

Cette  idee  a  motive  I’etude  qui  suit,  dans  laquelle  on  recherche  des  filtres 
approches  pour  un  probleme  de  filtrage  non  lineaire  discret  avec  petit  bruit 
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d’observation. 

On  considere  la  discretisation  la  plus  simple  des  equations  (1)  et  (2): 

Xfc+i  =  Xk  +  f>(Xk)  At  +  (T  \/Ki  ,  Xq  =  (  (3) 

yk  =  h{Xk)  +  -^^^Wk  (4) 

oil  Xk  est  une  approximation  de  Xtj,  (tk  =  fc  At),  {wk)  et  {tn*}  sont  des  bruits 
blancs  gaussiens  standards  independants  et  ^  une  v.a.  independante  de  {t/?*}  et 
{iDfc}.  On  se  propose  de  resoudre,  de  fagon  approchee,  le  probleme  de  filtrage 
associe  au  systeme  discret  (3  -  4),  ou  At  est  pris  comme  etant  un  parametre 
dependant  de  e. 

Soit  {A'jt}  le  filtre  optimal  pour  le  probleme  discret,  i.e., 

=  E[Xk\Y>‘]  ,  r*  =  i  =  0, 1, . . . ,  fc). 

Puisque,  dans  le  cas  general,  la  determination  de  Xk  presente  une  grande  com- 
plexite,  notre  but  sera  de  construire  une  “bonne  approximation”  de  ce  filtre 
representant  un  compromis  entre  les  couts  en  temps  de  calcul  et  la  precision 
des  resultats.  Si  s  =  0  et  A  est  inversible,  le  signal  est  observe  exactement: 
Xk  =  h~^{yk).  Pour  At  fixe,  quand  e  — »  0, 1’observation  est  faiblement  bruitee 
et  un  “bon”  filtre  consistera  a  prendre  tout  simplement  Af*  =  h~^{yk).  Cepen- 
dant,  notre  interet  portera  sur  des  resultats  asymptotiques:  on  fera  tendre  At 
vers  0  en  meme  temps  que  £  et  on  donnera  une  approximation  {Mk}  du  filtre 
optimal  en  temps  discret.  On  prendra  At  =  e“  et  on  s’interessera  a  I’estimation 
de  jA'"*  —  Mk\  quand  e  devient  “petit”.  Quand  e  — »  0,  on  aura  a  distinguer  3  cas, 
selon  la  valeur  de  a.  Les  cas  a  >  1,  a  <  1  et  le  ceis  intermediaire  a  =  1  presentent 
un  comportement  different.  Quand  a  est  “grand”  on  doit  pouvoir  approcher  le 
probleme  de  filtrage  en  temps  continu  (1  -  2).  En  revanche,  quand  a  est  “petit”, 
le  systeme  (3  -  4)  ne  sera  plus  necesseiirement  une  “bonne  discretisation”  du 
systeme  en  temps  continu  done  le  filtre  n’aura  pas  le  meme  comportement. 

On  commencera  par  etudier  le  cas  lineaire  dans  la  section  2,  on  proposera 
des  filtres  approches  unidimensionnels  possibles  et  des  resultats  asymptotiques 
sur  leur  “qualite”,  Cette  etude  nous  suggerera  le  filtre  approche  a  considerer 
dans  le  cas  non  lineaire,  ce  qui  est  fait  dans  la  section  3,  et  la  generalisation  des 
resultats  asymptotiques  de  “qualite”  du  filtre  est  obtenue. 

On  s’interessera  en  premier  a  I’estimation  de  Xk  —  Mk,  ce  qui  nous  ramene 
simplement  a  I’etude  d’une  serie  recurrente.  Ensuite,  on  procedera  a  I’estimation 
de  Xk  —  Mk-  Dans  le  cas  lineaire,  puisqu’une  equation  pour  Xk  est  disponible  (le 
filtre  de  Kalman-Bucy),  I’etude  d’une  serie  recurrente  permet  encore  d’aboutir 
au  resultat  voulu.  Par  centre,  dans  le  cas  non  lineaire,  la  methode  sera  toute 
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autre:  on  commencera  par  introduire  des  changements  de  probabilites,  on  utilis- 
era  une  version  discrete  du  Theoreme  de  Girseinov  (voir  annexe  A)  et,  par  deriva¬ 
tion  par  rapport  a  la  condition  initiale,  on  se  ramenera  a  I’etude  d’une  expression 
asymptotique  pour  Xk  —  Mk  dont  on  prendra  I’esperance  conditionnelle. 

La  section  3.4  contient  un  bref  resume  sur  les  filtres  approches  a  considerer 
dans  le  cas  non  lineaire,  selon  la  valeur  de  a,  et  des  resultats  sur  leur  efBcacite. 
On  termine  cette  etude  par  un  essai  numerique,  presente  dans  la  section  4. 

Notation  1.1  Dans  la  suite,  on  notera  c  ou  C  n’importe  quelle  constante  pos¬ 
itive  independant e  de  At  et  de  e. 


Remerciements:  J’aimerais  remercier  E.  Pardoux,  Professeur  a  I’Universite 
de  Provence,  pour  m’avoir  suggere  ce  travail  et  J.  Picard,  chcirge  de  recherche 
a  I’INRIA,  pour  les  suggestions  et  discussions  fructueuses  qui  m’ont  permis 
d’avancer  dans  ce  sujet. 


2  Cas  lineaire 


On  considere  le  systeme  lineaire: 


oil 


f  Xfc+i  =  (1  +  bAt)Xk  +  (Ty/At  t/Jjt+i ,  Xo  —  ^ 


yk  =  hXk  + 


VKi 


Wk 


(5) 


(Hi)  {lOfc}  et  {«)*:}  sont  des  bruits  blancs  gaussiens  standards  independants, 

(H2)  (  est  une  v.a.  gaussienne  independante  de  et  {ujjt},  de  variance  po, 
oil  Po  est  une  constante  ne  dependant  pzis  de  e  ou  de  At. 


(H3)  dh  >  0. 

Remarque  2.1  L’hypothese  (H3)  peut,  en  fait,  etre  remplacee  par  ah  ^0  et, 
dans  ce  cas,  on  remplacera,  dans  les  expressions  du  gain,  ah  et  a/h  par  leurs 
valeurs  absolues.  Neanmoins  I’hypothese  (H3)  n’est  pais  plus  restrictive  que  cette 
derniere,  puisque  si  on  remplace  par  —a  dans  le  systeme  (5)  on  obtient  des 
modHes  equivalents. 


II  est  bien  connu  que,  dans  le  cas  lin^ire,  la  loi  conditionnelle  de  Xk  sachant 
est  gaussienne  et  I’estimation  optimale  est  done  donnee  par  des  equations 
en  dimension  finie,  les  equations  du  filtre  de  Kalman-Bucy: 


Pk+l\k 


(1  +  bAt)Xk-i  + 


hAipk\k-i 

e'^  +  h'^  At  pk\k-i 


{Vk  —  /*(!  +  bAt)Xk-i) 


(1  +  6At)Vpt|fc_i 

Atpk\k-i 


(6) 

(7) 


avec  les  notations: 


Pk  =  E{{Xk  -  Xk)^]  ,  Xk\k.i  =  ElXklV^-^]  et  pkik-i  =  E[{Xk  -  Xqfc-i)']  • 


En  outre, 

_ ^  Pqfc-l  /j;,\ 

*  +  h'^Atpk\k-i  ' 

Notre  objectif  est  de  construire  un  processus  {Mk]  qui  approche  {X^}  et  qui 
est  plus  simple  a  calculer  que  ce  dernier.  Dans  ce  but  on  commencera  par  etudier 
le  gain  stationnaire,  pour  passer  en  suite  a  la  recherche  d’une  approximation  de 
ce  meme  gain,  selon  les  differentes  valeurs  de  a.  Le  gain  du  filtre  ayant  un 
rapport  avec  la  variance  de  I’erreur  d’estimation  on  deduira  simtiltanement  des 
resultats  sur  cette  deuxieme  quantite  et  ses  approximations. 
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2.1  Calcul  de  quelques  expressions  utiles 

Dans  cette  section  on  calcule  certaines  expressions  qui  nous  seront  necessaires 
dans  les  sections  suivantes. 


Calcul  du  gain  stationnaire 

On  commence  par  determiner  la  covariance  de  I’erreur  de  prediction  pjt|jfc_i  dans 
une  situation  stationnaire,  ce  qui  est  equivalent  a  calculer  la  valeur  stationnaire 
de  la  covariance  de  I’erreur  d’estimation  pk  (not^  p+),  puisqu’on  passe  de  Tune 
a  I’autre  par  I’expression  (8). 

On  notera  6k  le  gain  a  I’instant  tk  : 

A  hAtpk\k-i  _  hAi 

*  h}Atpk\k-\ 


Soit  p,  la  valeur  stationnaire  de  Pfc|jk_i,  i.e,  la  solution  de 

(l+6At)Vp,  2 


Ps  = 


+  a^At . 


On  obtient 


ou 


Ps  = 


£2  +  At  p3 

-[e^  —  (1  +  bAtYe^  —  Ai^]  +  Atp{e,At) 
___  ^ 


p(e,A<)=  [[(26  +  6^A<)  £^  +  ^  . 


Ainsi 


(26  +  6^  At)  +  cr^/i^  At  +  p{e,  At) 

_ 


Ps  = 

Le  gain  stationnaire  sera  done, 

_  hp,At  _  h  At  ^ 

+  h?  Atpg 


6, 


Le  filtre  approche 

On  considere  un  schema  qui  resulte  de  (6)  en  remplagant  Ok  par  sa  valeur  station¬ 
naire  O3  ou,  d’une  fa^on  plus  generate,  par  une  approximation  de  cette  valeur. 
Dans  la  suite,  6  designera  done  soit  6,  soit  une  approximation  de  0,, 

On  considere  le  processus  {Mk}  donne  par  I’expression 

Mk+i  =  (1  -f  6  At)Mk  -f  0  {rjk+x  -  h  {I +  b  At)Mk) ,  Mo  =  mo  .  (9) 

On  obtient  les  resultats  qui  suivent. 

Le  premier  resultat  concerne  les  differences  0,  et  pk\k-\  -  Ps- 
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Lemme  2.2  Soil  rjig  =  Ok  —  9a  et  Hk  =  Pk\k-\  —  Pb-  On  a  Iss  expressions 


Vk+l 


hAte' 


et 


[e^  +  At  p,] 

C2 


Pk+l 


H  =  -) 
C2 


Cl/^0  +  ^^2 

i=0  C2 


Po 


ou 


Co  =  (l  +  6Af)^e‘‘ 

Cl  =  At  {e^  +  At  Pa) 
C2  =  (£■*  +  Ai  p,)*  . 


Preuve  De  la  clefaiition  cle  rjk  vient  que 

hAtpk+i\k  hpaAt 


Vk+\  = 


+  h'^  At  pk+i\k  e'^  +  h'^Atpa 

hAt[e^Pk+i\k  +  h?Atpk+iikPt-e'^pB-h'^  Afpfc.n|fep,] 

(£2  +  ^2  Afp*4.i|ii..][£2  +  /i2  p,] 

d’ou  I’expression  de  rjk+i- 

Quand  a  pk,  on  a  la  formule  recursive 

(1  +  b  At)^  pkik-i  _  {l  +  bAt)^c'^Pa 
e'^  +  Atpk\k-i  e'^  +  h'^Atpa 
il+bAt)^£* 


Pk+l  = 


(£2  +  /l2  At  Pk  +  h}  At  P,](£2  +  h?  At  Pa] 
Co 

- Pk  • 

C\Pk  +  Cj 


Pk 


Done 


Pk^rl  = 


Co  1 

C,  .  ^2 
/2jk  +  — 
Cl 


Suivant  un  raisonnement  par  recurrence  on  trouve  I’expression: 

C2 


Pk  = 


_ ^ _ 

Cl/20  S  C^"^"‘<<,  +  <4 
(=0 


/^o  =  (-) 
C2 


ClA^O  ]C(“)’  +  ^2 
1=0  *^2 


Po 


(10) 
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Lemme  2.3  Soil  rik  —  9k  —  0.  On  a  Vegaliti 

E[{Xk^,-Mk^rY\  =  {\-^ht^tf{\-hOfE[{Xk-Mkf] 

£2 

Wk+i  +  , 


oii 


I 


\-h6 


c-2 

•«> 

+  h^^tps 


+  k{6,-  9). 


Preuve  D’apr^  les  expressions  (6)  et  (9), 

Xk+i  —  Mk^i  =  {l  +  b  Ai)Xk  +  Ok+i{yk+i  —  h{l  +  b  At)Xk)  —  {I  +  b  At)Mk 
-9  {Vk+i  -h{li-  bAt)Mk) 

=  {I  +  bAt){Xk  —  Mk)  +  {9k+i  —  9)yk+i 
-h{l  +  bAt)i9k,t^k-dMk) 

Etant  donne  que  9k+i  =  9  +  —  5),  on  a 

=  (l  +  6A0(l-/i^)(A’*-Mfc) 

+(^i+i  —  ^)(y*+i  —  h  {I +  b  At)Xk), 

ou  j/fcj-i  ~  h(l  +  bAt)Xk  =  Vk  est  I’innovation  {uk  independant  de  K*)  et 

£2 

Evl  =  A.2pfc+i|t  +  —  , 

d’oii  le  lemme.  ■ 


Expressions  asymptotiques  de  et  de  E[{Xk  -  A/*)*]  quand 

le  pas  de  temps  est  e" 

On  prendra  dans  la  suite  At  =  ,  a  >  0.  Notre  but  est  de  faire  une  discussion 

sur  I’ordre  de  grandeur  de  I’erreur  associee  au  schema  (9),  pour  les  differentes 
valeurs  de  a. 


Remarque  2.4  On  rappelle  que  : 


oil 

Done 


(26  +  6^  At)  +  <7^  At  +  Pc(£,  At) 
__ 


Pai^)  =  [  ({26  +  6^  At)  £^  +  At]*  +  4  o*  h“  £*  ]*/* 


c{e^y  e)<Pf<C  (£®  V  £) 
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Alors,  puisque  /io  ^  0  (i-e-  la  valeur  initiale  de  la  covariance,  po>  est  une 
constante  independante  de  e),  des  calculs  precedents  il  vient  que  : 

0  <  /iit  <  (— )  Ho  , 

C2 

oil 

Co  _  (1  +  6e°')^£'* 

C2  [e^  +  £“  ps]^  ^ 

i.e.  Co  <  C2,  puisque 

Co  _  [g^  +  h'^e^’pa  -  (1  +  6e°)e^][£^  +  h}  p^  +  (1  +  &e°)g^] 

C2  {g^  +  g“  p»Y 

e°  [h'^Pa  -  6g^][£^  +  k^e^ps  +  (1  +  6g*)g*] 


g^  +  g“  paY 


pour  £  “assez  petit”. 


On  etablit  la  proposition  suivante: 


Proposition  2.5  On  a  I’estimation 


Cexp{-c-^}  +  Cg  , 


$i  a>  1 


C exp{  [4  ^^(1  —  Of)  +  2(1  -  a)]  logg  }  ,  st  a  <  1  . 


Preuve 
Par  definition, 


E[{Xk+^  -  ^  p,+, 


(Hk+i  +  Ps) 

g2  +  A2e?  {pk+l  +Ps) 


<  £  ■  £^5 

~  g2  +  ^2  gO  ^  ^2  ^  f,2  gor  p^ 


C  {colciY'^^  +  C  £  ,  si  a  >  1 

Cg^l*“®l(co/c2)*'^^  +  C g^““  ,  si  a  <  1 


d’ou  la  proposition. 


D’autre  part,  du  lemme  2.3,  on  obtient  le  lemme  suivant: 
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Lemme  2.6  On  a  Vestimation 


/VT* 

E[{Xk^,  -  <  £(1  -  , 


A  =  i-(i  +  bs“)^(l-k0y 

„  A  _2  e“  PJt+i|it 

^fc+i  -  Vk+i  ^ 

Bo  =  El(Xo-Afon 


Preuve  Le  lemme  2.3  devK  i»; 


Ei(x,+r  -  =  (1  -  i)  -  M,y] + p,+, 


2.2  Filtre  utiJisant  le  gain  stationnaire 

On  prend  comme  approximation  du  gain  0;.  le  gain  stationnaire  0,,  i.e.  0  =  0g 
(voir  le  schema  (9)).  On  deduira  par  la  suite  des  resultats  pour  d’autres  types 
d ’approximations  0  a  preciser. 

Dans  ce  cas,  on  obtieat  la  proposition  suivcinte: 

Proposition  2.7  Pour  le  schema  (9)  avec  0  =  0,  on  a  I ’estimation 


E[iX,+,  -  < 


C exp{— , 


si  a  >  1 


C  exp{4  ^^(1  —  a)  logf:}  ,  s:  a  <  1  . 


Preuve  Suivant  les  notations  du  lemme  2.6,  on  a 

1-.4  =  {l^he°f{\-h0,f 
(1  +  6£°)^£^ 

(£2  +  h?  £°  p,y 

_  Co 


Bi+i  — 


2  + 


/ie"+* 

[£24.  /i2  gor  p,.^j|.](e2  +  h}  e“ p,] 

_ _ 2 

[g2  ^  /l2e«p;^l|,][£2  T-  42e“p,]2  ^‘+1 

et,  puisque  Pi+i\i  >  p*,  on  a  la  majoratk:. 


^•■+1  -  [g2  4.  Fgap,]3^«+1  -  Bf^i+1  > 


£2  +  ^2  gQ  p.^^n 


A 

(£2  -  ^2  gO  p^j3 


= 


Ol  po  Z)jso(^/ <^2)^  +  C2 


et  i/o  est  tel  que 


D’aprfej  ce  lemme, 


E[{Xo-Mon  =  BHopl. 


E[{Xk+i  -  <  {C0/C2)  -  Mk?]  +  ^[(co/c2)^]''+^  pI 


<  E(A)/c2)''+‘-'5[(Q,/c2)*]‘i/,p2 

1=0 


On  cherche  maintenant  un  majorant  pour  la  serie 


= 


l-(Co/C2p  .  f’ 

C]  po — ; - 7 - !■  C2 

1  -  C0/C2 


^  ^  1 
avec  —  <  1 


=  (^) 
Cl  Po 


(1  +  -  Mc,Y*' 

Cl  Po 
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‘  ^  C,  /io  ^  h  [1  +  ic2  -  Co)/{c,  /Xo)  -  (Cb/C2)f 


On  commence  pax  calculer  I’ordre  de  grandeur  de 


C2-CQ 


C2-C0  _  (g^  4-  ptf  -  (1  4-  6e“)^£‘* 

Ci/xq  /Xo  (g^  + /i2£“P,) 

Soil.  Num  le  numerateur  du  membre  a  droite  dans  cette  expression: 

Nwn  = 

=  hY  e°  (p,  -  ^  g^)[g^  +  h'^  £<=’  ^  +  be'')  g^] 

Or  p,  —  —  g^  =  0(g"  V  g),  ce  qui  entraine  Num  sr  (P(£^“  V  g“+^)  et  done 


C2-C0 


C?(g“  V  g)  . 


De  plus 


C2  -  Co 


On  peut  faire  la  majoration  suivante: 


<-i 


^ _ (^0/^2)* _ 

1=1  (1  +  (C2  —  Co)/(cx  Po)  —  (co/c2)‘]^ 

r _ (coMr _ 

•A)  (1  +  (C2  -  Co)/(ci  po)  —  (Co/c2)^]^ 


log  (C0/C2)  [1  +  (C2  -  Co)/(Ci  Po)](C2  -  Cd)/(Ci  po)  ‘ 

Done 

it  4*1 

r-r-/"i  rr? — : + »•>  ■ 

i=0  Cl  Po  log  (coM)  1  +  (C2  -  Co)/{Ci  /Xo) 

De  I’expression  (11),  il  vient  alors  que 

Ei{X,„-MMy\<4B{cc/c]''*'i‘^-’^  ^ - L_ - - 

Cl  Po  log  (C0/C2)  1  +  (C2  -  Co)/{Ci  Po) 


+A. , 
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Etant  donne  que 


et 


B 


_  f  , 


si  a  >  1 
si  O'  <  1 


log  (co/cj) 
on  obtient  la  majoration: 

-1  1 


1 _ ^  1  _  r  ,  si  Q  >  1 

0/C2)  ~  1  -  (co/cj)  ~  \  0(1) ,  si  o  <  1 


B 


log  (co/ca)  1  +  (c2  -  co)/(ci  fio) 


4.  H  ^  '  sia>l 


et 


1 


BHo  =  -^E{{Xo-Moy]^ 

/^o 


Done, 


si  a  >  1,  puisque  1  —  Co/ca  =  C7(e""‘),  on  a  I’estimation 
E[{Xk+i  -  <  Cexpi-c^}  , 

si  a  <  1,  puisque  cq/ci  =  on  a  I’estimation 

E[{Xui  -  Mfc+i)')  <  C7exp{(4(^  +  1)(1  -  a)]loge} 


(12) 


Remarque  2.8  On  peut  constater  que,  quel  que  soit  a  >  0,  les  valeurs  sta- 
tionnaires  de  E[{Xk+i  -  Mk+i)^]  et  E[{Xk+i  —  sont  du  meme  ordre  de 

grandeur  done  I’utilisation  du  filtre  approche  (voir  le  seheina  (9))  est  justifiee. 


2.3  Filtres  utilisant  une  approximation  du  gain  station- 
naire 

On  eherche  maintenant  les  expressions  plus  generales  qu’on  obtient  quand  on 
utilise  une  approximation  5  du  gain  asymptotique  0,. 

Pour  la  eonstruction  de  eette  approximation  on  va  d’abord  considerer  des 
approximations  pour  la  covariance  p,.  Elies  conduisent  naturellement  a  des  ap¬ 
proximations  pour  le  gain  9^. 

On  va  estimer  I’erreur  Xk  —  Mu- 
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Proposition  2.9  Supposons  que  Pa  —  P  =  0{e^)  ,  m  >  1,  i.e.  p  est  une 
approximation  de  p,  d’ordre  e"'. 

Alors, 


Cexp{-c^}  +  Ce2’"-^  , 

£ 

Cexp{4^(l-a)loge}  +  C7e2”>+<-5“  , 


51  a  >  1 

St  O'  <  1  . 


Preuve  D’aprfe  le  lemme  2.6,  on  a  I’estimation 

k^X 

E[{X^^^  -  Mt+i)')  <  E(1  -  , 


oil 

{l  +  be°ye^ 

[£2  +  /i2£»p]2 

et 

=  rii  <  Epi  , 

avec 

[s*  +  h'^  €“  p,](e^  +  h'^  £“  p]* 

et 

Pi  =  fii  +  {Ps  -  P)  • 

Done 

fc+1 

1=0 

Ar+l  D 

i=0 

En  utilisant  les  majorations  dans  la  demonstration  de  la  proposition  2.7,  on 
obtient  la  majoration: 


k-X-X 

E[{X,+r  -  <  2BpI  (1  -  ^ 


'(co/cajT 

1 


Hi  +  Ce^”' 


B 
A  ■ 


Maintenant,  si  p,  —  p  >  0  alors  1  —  i4  >  Co/c2 ;  par  contre,  si  p,  —  p  <  0  alors 
I  ~  A  <  co/c2.  Par  consequence,  on  a 
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CQ 


+  Ce2'":| 

A 

Quelques  calculs  nous  donnent: 

Q  (C2  -  Co)/(Ci  Ho) 

1  +  (C2  -  co)/(ci  Ho) 

_ _ 

[e^  +  e“  Va][^^  +  p]^ 

_ (p,  -  6/A^  £^)[£^  +  h}  e°  p,  +  (1  +  6^°)  £^] _ 

■  Po  £"  Pi)  +  (Pi  -  e2)[£.2  +  ^2  £.a  ^  ^  ^  J.a)  ^2 


J  0(e“  ^)  ,  si  Qt  >  1 
\  ,  si  a  <  1 

£*  [/i^p  —  6£^](e^  +  fi^e^p  +  il  +  be°)  e^] 

(e^  +  /i2  £<»  p]2 

oil  h^p-be^  <  C(e"  Ve) 


r  C?(£"-'),  a>  1 

10(1),  a<l 

h}  e"* 

+  h}  £°  p,\[li^  p  —  6e2j[£2  4.  ^2£apq.  q.  6£®)e2j 

r  0(l/£),  a  >  1 
\  C?(e^-5“),  a  <  1 

On  obtient  done  les  estimations  suivantes: 

•  si  a  >  1  , 

£((A:fc+,  -  Mfc+i)^!  <  Cexp{-c^}  + 

£ 

•  si  a  <  1  ,  puisque  1  —  ^4  <  ©(e**!*"®!), 

E[{Xk+t  -  Mfc+O'l  <  C'exp{4  ^(1  -  a)  \oge}  + 
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2.3.1  Approximations  obtenues  pour  la  covariance  p. 

On  utilise  un  developpement  limite  de  Pa(e)  pour  construire  une  approximation 

p  de  p,  et  on  obtient  done  une  approximation  9  de  par  9  =  -  ,  - . 

+  h^e°p 

On  explicite  quelques  expressions  possibles  pour  I’approximation  de  la  vari¬ 
ance  p; 

Puisque,  par  hypothese,  /i  >  0,  on  a  que 


(a)  si  o  >  2,  le  developpement  de  Taylor  de  Po(^)  i^ous  donne,  apr^  quelques 
calculs, 

Pa(e)  =  2  A  (7  £  -f  O(e^) 

et  done 


<T  6  2  CT‘ 


Si  on  prend,  par  exemple, 
1 


p  =  ajh  e  ,  i.e.  -  p  =  0{e*) ,  alors 


E[{Xk  -  M,)^]  <  Cexp{-cj}  C£3 


(13) 


2. 


p-alhe  +  blh^e^  si  a  ^  2  ,i.e.  p,  -  p  =  V £^)  ,  alors 


E[{X,  -  M,f]  <  Cexp{-c  -f  C  (£^“-‘  V  e®) . 


3. 


p  =  crfhe  +  b/h'^e'^  -f  a^l2e°  si  2  <  a  <  3  ,i.e.  p,-p  =  O(e^),  alors 


E[{Xk  -  Mkf]  <  Cexp{-c-^}  +  Ce^. 
(b)  si  1  <  Q  <  2,  du  developpement  limite  de  Po(£)  il  resulte  que 


P,  =  7  ^  e"  +  0(e*  V  £2“-‘) 

n  I 

done  on  peut  prendre  comme  approximation  de  p,  par  exemple: 

1. 


p  =  ajh  e  ,  i.e.,  p,  —  p  =  <!7(e‘')  et  alors 


(14) 


E[{Xk  -  Mkf]  <  Ccxp{^c^)  T 

£ 
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2.3.2  Approximations  obtenues  pour  le  gain 

Puisque  notre  interet  porte  sur  la  construction  de  schemas  du  type  du  schema 
(9),  on  essaie  de  decrire  une  approximation  du  gain  stationnnaire  B,.  En  fait  ces 
approximations  peuvent  aussi  etre  obtenues  directement  d’apres  le  developp- 
ement  de  /90(e),  sains  passer  pair  I’approximation  de  la  variamce. 

On  a 


0.  =  T 


1  1 
hl  +  Do 


Do  = 


{2b  +  b'^  e")  e*  +  <7*  h?  e®  +  /)a(e) 
Pa{e)  >c(c“  Ve) . 


•  Premier  cas:  a  >  1 


^0  ~  T - TTi  > 

h  cr  £*  * 

puisque  Pa{£)  =  2hae  +  0{£^  V 
Une  approximation  de  Bg  est  done 

d  =  (Tg®-^  I .  (15) 

Proposition  2.10  Si  a  >  I,  avec  le  gain  approchi  donni  par  (15),  alors  on  a 
Vestimation 

£|(Xhi  -  Mu.,  )=|  <  C  exp{-c  ^i±l)  +  C  (£=”*■  V  e") 


Preuve 

Puisque 

et,  dans  notre  cais, 


hp£° 

e^  +  e®  p 


^  h{l-<7he-^)' 

on  obtient,  en  utilisant  le  developpement  limite  de 


1  —  <7^C®“^’ 


p  =  •-  e  +  <r^  e®  +  (7^  A  +  Cl(c^“~*) . 
h 

Par  consequent 

p,  -  p  =  c  (e^  V  e®)  et,  si  a  >  2  ,  p,  —  p  >  0  . 
Done,  de  la  proposition  2.9,  on  deduit  les  estimations  suivantes: 
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r 


pour  1  <  a  <  2, 

E[{Xk+x  -  Mk+i)^]  <  6-exp{-c^}  + 


pour  a  >  2, 


Ept+i  -  M/b+.)’l  <  Cexp{-c^itl}  +  Ce= 

e 


Remarque  2.11  Une  fois  encore  on  trouve  des  estimations  de  —Mk+\ )^] 

d’ordre  inferieur  ou  egal  a  celui  de  E[{Xk+i  — 

Remarque  2.12  Le  schema  correspondant  aux  approximations  (13)  et  (14) 

■)  n’a  pas  d’interet  pratique,  puisque  son 


ae 


a+l 


(i.e.  qui  utilise  le  gain  6  —  ^ 

'  ®  +  0^  /i  ^ 

utilisation  oblige  a  un  calcul  plus  complique  que  celui  du  schema  qu’on  vient 
d’obtenir  alors  que  I’ordre  de  I’erreur  associee  reste  le  meme. 


•  Deuxieme  cas  :  a  <  1 


Do 


,2(I-a) 


Done  une  approximation  du  gain  stationnaire  Ot  est 


Le  schema  (9)  devient  alors 


^k+\  =  ^  J/fc+i 


(16) 


•  Troisieme  cas  :  a  =  1 


Do 


a  h  \/4  4-  0^  h} 
Quand  e  — >  0,  une  approximation  de  0^  est  encore 


e . 


et  on  retrouve  le  schema  (16). 

Proposition  2.13  Si  a  <  1,  poxirun  tel  schema,  on  a: 

1 


! — 
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Preuve 


i)  D’apres  le  lemme  2.3, 

Xk»  -  _  A(1  + 


done 


£((X,+.  -  = 


ft2£“|£“  +  A’£“Pt+ln) 

d’ou  la  premiere  expression, 
ii)  D’autre  part, 

Xfc+i  -  Mt+i  =  —^^72  '^l+i  . 

d’ou  la  deuxieme  expression  dans  la  proposition. 
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2.4  Recapitulatif 

Les  tableaux  suivants  nous  donnent  les  vitesses  de  convergence  a  0  de  I’erreur 
quadratique  moyenne  pour  les  difFerents  filtres  approches  qui  ont  fait  I’objet 
de  cet  etude,  mettant  en  evidence  la  dependance  de  ces  filtres  selon  le  rapport 
entre  la  variance  du  bruit  d’observation  et  le  pas  de  temps. 


•  Pour  a  >  1 


gain  du  filire  approchi 

esiimaiion  de  I’etreur 
E[{Xt+i  -  Mt+,)2] 

0,  (gain  stationnaire) 

Cg-e^+iA 

C(£2“-W£3) 

-  £2  +  <r/i£<«+>+<r2/i2/2e2“ 

C(£3V£''“-3) 

(r£“+i+6//,e«+2 
£2  (7  +  6e“+2 

C(£2“-1  V£®) 

.3  <r£<»+‘+6//.£“+2+*72/,/2e2» 

“  £2  ^  ^  JjO+2  ^  /,2^2e2“ 

Ce^ 

4 


1 


D’autre  part,  on  a: 


E[iXk+i  -  Xk+if]  <  +  Ce  . 


•  Pour  a  =  1 


gain  du  fiUre  approchd 

estimaiion  de  I’erreur 
E[(Xt+i-Mk+in 

6,  (gain  stationnaire) 

1 

h 

Ce 

Ce^ 

■ 

Ce® 

■  4  "2  vr+w/r 

1  +  /i(o- 

/  <r^  <T^fi  bch  . 

¥‘■*■4+2  ’  +  <‘+2v'1+^HV4'‘ 

D’autre  part,  on  a: 

E[{Xk+i  -  Xk+if]  <  +  Ce. 


•  Pour  a  <  1 


join  rfu  ^/<re  approchi 

estimation  de  I’erreur 
E[{Xk^,  -  Mk+iY] 

0,  (gain  stationnaire) 

QgCU+i/t{l-a)\ost 

1 

h 

Ce*-^^ 

e^® 

Cre<‘*+>/t(i-or)iogt  ^  (758-70 

D’autre  part,  on  a: 

E[(Xk+i  -  Xk^if]  <  +  Ce2-“  . 
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Remarque  2.14  Des  calculs  analogues  nous  permettent  aussi  d’estimer  di- 
rectement  Xk  —  Mk. 
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Cas  non  lineaire 


On  suppose  maintenant  que  les  fonctions  b  et  h  sont  non  lineaires. 
On  considere  le  systeme  decrit  par  les  equations 

f  Xk+i  =  Xk  +  b(Xk)  A<  +  (7  sfKi  Wk+t ,  Xo  =  C 

I  »  = 

oil 

(H2’)  ^  est  une  v.a.  de  loi  de  probabilite  pa  telle  que: 

Po  6  /|— |’■po(^c)da:  <  oo  Vr  >  0  . 

J  Po 


(17) 


Pour  des  raisons  liees  a  I’application  d’une  methode  de  changement  de  pro¬ 
bability,  pour  obtenir  les  estimations  de  I’erreur  quadratique  moyenne,  il  sera 
plus  commode  de  consid^er  notre  observation  sous  la  forme  suivante: 


j/fc+i  =  h{Xk)  + 


(18) 


oil  {ufc}  est  un  bruit  blanc  gaussien  standard  independant  de  {uiit}. 

Formellement,  I’observation  yk+i  contient  la  meme  information  que  I’obser- 
vation  yk- 

On  suppose  que  F*'  est  la  tribu  des  observations  jusqu’a  I’instant  tki 


=  <^iyi,y2r--,yk) , 


i.e.  =  F*-^ 

On  veut  etudier  la  “qualite”  de  I’approximation 

Mk+i  =  Mk  +  b{Mk)  +  0  {yk+i  —  h{Mk))  ,  Mo  =  mo  ,  (19) 

correspondante  a  une  ^ape  de  prediction:  Xk  =  E[Xk\Y'^] .  Le  gain  0  sera  donne 
par  une  expression  a  preciser  par  la  suite. 

Des  rfeultats  qu’on  obtiendra  sur  I’erreur  de  prediction  Xk  —  Mk  on  deduira 
immediatement  les  resultats  equivalents  sur  I’erreur  de  filtrage,  puisque,  voir  le 
paragraphe  3.1.3,  on  a  que 

|A'fc_i-AM<cA<, 
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ohXk-i-=E[Xk-i\Y'^]. 

Pour  la  construction  du  filtre  approche  on  va  s’inspirer  de  I’etude  du  cas 
lineaire  et  ainsi  choisir  une  approximation  du  gain.  L’estimation  de  Xk  —  Mu 
I  sera  obtenue  a  partir  de  I’etude  d’une  serie  recurrente,  En  revanche,  pour  estimer 

‘  Xk  —  Mk,  on  introduira  d’abord  deux  changements  de  probabilites  et  on  obtien- 

dra  une  expression  asymptotique  pour  Xk  —  Mk  qui  fait  intervenir  les  derivees 
de  X,  par  rapport  a  la  condition  iriitiale.  Sur  cette  expression,  on  conditionnera 
j  les  differents  termes  par  rapport  a  la  tribu  des  observations  et  on  procedera 

I  a  des  estimations.  Selon  la  valeur  de  a,  les  changements  de  probabilites  choisis 

s  seront  differents  et  on  aboutira,  evidemment,  a  des  expressions  asymptotiques 

i  differentes  pour  Xk  —  Mk.  Dans  une  premiere  section  on  etudie  le  cas  a  >  1 

(section  3.1)  et  en  suite  le  cas  a  =  1  (section  3.2).  Le  cas  or  <  1  offrant  certaines 
difficultes  dans  I’application  de  cette  methode  de  demonstration,  on  se  limitera 
a  faire  une  remarque  proposant  un  filtre  approche  “performant”. 

Hypotheses:  En  plus  des  hypotheses  (Hi)  et  (H2’)  on  suppose  que: 

(H3’)  h  est  une  fonction  a  derivees  bornees  et  ||<tA'||  >  ffh'{x)  >  c;,  >  0  ,  Vx. 

(H4)  6  est  une  fonction  a  derivees  bornees,  |6'(i)|  <  ||6'1|  ,  Vx. 

(H5)  Ai  =  e"  ,  a  >  1. 

Remarque  3.1  La  remarque  2.1  reste  valable  dans  le  cas  non  lineaire,  c’est 
a  dire,  pour  des  raisons  identiques,  la  condition  “crA'(x)  >  0 ,  Vx”  peut  etre 
remplacee  sans  aucune  difficulte  par  “ak'(x)  ^  0,  Vx”.  Ce  qui  est  necwsaire 
est  que  la  fonction  h  soit  supposee  injective. 


3.1  Le  cas  a  >  1 

On  considere  d’abord  le  cas  oii  At  =  e®,  avec  a  >  1. 

On  propose  un  filtre  approche  unidimenslonnel  et,  pour  ce  filtre,  on  estime 

Xk-Mk^iXk-Mk. 

Le  filtre  approche 

En  analogic  avec  le  cas  lineaire,  on  prendra  0  =  a  —  dans  le  schema  (19),  i.e. 
Mk+i  =  Mk  +  b{Mk)At  +  a^{yk+i-h{Mk)),  A/o  =  mo  .  (20) 


I 
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3.1.1  Estimation  de  Xk  —  Mk 

On  obtient  facilement  le  resultat  suivant,  lequel  generalise  le  resultat  dans  le 
cas  lineaire: 

Proposition  3.2  Le  filtre  approche  Mk  etant  donnepar  (SO),  on  a  ^estimation 


E[{Xk  -  Mkf]  <  C  exp{-cj}  +  Ce. 


Preuve 

On  pent  utiliser  les  developperaents  de  Taylor  des  fonctions  b  et  h 


(21) 


et 


b{Xk)  =  biMk)  +  b\C^)iXk  -  Mk) 
h{Xk)  =  h{Mk)  +  h'U^){Xk  -  Mk) , 


pour  obtenir  I’egalite 


AT+i-Mi+i  =  [l-c~h'{^^)  +  b’{<:^)At]{Xk-Mk) 

£ 

+cr\/Ai(tOk+i  -  Vk+i)  . 


(22) 


D’apres  les  hypotheses  (H3),  (H4)  et  (H5)  et,  vu  que,  d’apres  (Hi)  et  (H2’), 
Xk  —  Mk  et  zxjjt+j  —  Vk+i  sont  independants,  on  a,  pour  e  “assez  petit”  (e  < 

At 


E[{Xk+i  -  Mk+i)"]  <  (1  -  CA  —  +  |16'1|  At)^E[{Xk  -  Mk)^]  +2a^At 

£ 


Soit 


/l  =  1-(1  -  cft  — +  i|6'||A0' ,  B  =  2(r^At 


Avec  ces  notations, 

EKXn.,-Mwf]  <  {l-A)‘-*'E[(Xo-Mof]+B'£H-Ay 

1=0 

B 


puisque  1  —  A  <  1. 
Mais 


c^<A<C(—) 

£  £ 


(23) 
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done 


B  s 
-  =  0{e) 

et  (23)  devient: 

E[{Xk+i  -  Mk+i?]  <  C  exp{-c^}  +  Ce 
ou,  plus  precisement, 

E[iX,  -  M,)^]  <C{1-A)^  +  Ce.  (24) 


3.1.2  Estimation  de  Xk  —  Mk 
On  montre  le  theoreme  suivant: 

Theoreme  3.3  Le  schema  (20)  verifie: 

Xk  -  Mk  =  V  e) 

an  sens  oil 

E[\Xk  -  Mjtl]  <  Cexp{-c^}  +  C  V  e) . 

Preuve 

Elle  sera  divis^  en  plusieurs  parties,  utilisant  des  changements  de  proba- 
bilites,  une  version  discrete  du  Theoreme  de  Girsanov  et  la  derivation  par  rap¬ 
port  a  la  condition  initiale. 


Changement  de  probabilites. 

On  introduira  deux  changements  de  probabilites  (etapes  (a)  et  (b)). 

(a)  Le  I**’’  changement  de  probability  affectera  la  loi  de  v. 

On  considwe  la  probabilite  P  definie  par: 


avec 


Lf  =  «p{  HX,-,)]  Vi  -  i  } 

1=1  *  ^  1=1  ® 
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On  utilise  I’equation  de  I’observation 


\/^t  -  f  'v/^^  L/V  \1 

- yk  =  Vk~  I - n{Afc_i)  , 

£  £ 


oil  Xk-i  est  ^k-i-  mesurable. 
On  obtient: 


-  exp{  —  E MAVi) 2/.-  -  -  D  }  . 

^  i=t  ^  1=1 

D’apres  la  version  discrete  du  Theoreme  de  Girsanov  (voir  I’annexe  A), 

sous  P  (probabilite  de  reference),  to*  et  — —  yk  sont  des  Tk- 

£ 

bruits  blancs  gaussiens  independants.  [P  est  equivalent  aP  dans 
chaque  Pit.) 

(b)  Le  2'^"*'  changement  de  probabilite  va  affecter  la  lo;  de  ly. 

On  definit  la  nouvelle  probabilite  P  par; 


A  -  , 


Afc  ‘  =  exp{  Yh  (^( A,-i )  -  )]  ly, 

.•=1  ^ 

^  «=i  ^ 

Soit  {tyt}  le  processus  defini  par 

Wk  =  Wk-  (A(Xit_i)  -  h{Mk-i)]  . 

oil  Xk-i  —  Mk-i  est  Pit-i-niesurable. 

Alors 

\k  =  exp{ -  A{A/,_i )]«>,• 

^  «=i 

A/  ^ 


D’apr^  la  version  discrete  du  Theoreme  de  Girsanov, 

Sous  P,  Wk  et  yk  sont  des  bruits  blancs  gaussiens  independants 
et  Xk  verifie  I’equation: 

Xk+i=(r—[hiXk)-hiMk)]  +  Xk  +  b{Xk)At  +  <Ty/Aiwk+i  .  (25) 


La  densite  de  P  par  rapport  a  P  est  Lkt^k  done: 
V  v.a.  ij)  P-integrable  et  ^jt-mesurable, 


S  ^  log(i»At) , 


Puisque 


(KAV,)  -  6(M._i)]  +  —  [y.'  -  ft(A._0] 

<T  e 


5  =  EfMAi-i)  -  MA/.-i)][(X.-  -  Mi)  -  (X._i  - 

(TE 

t=l 

+—  E(A(A’,-i)  -  A(M(.,))[6(X,-.)  -  6{A/,-i)) 

.=1 


+f  i:  *("<-.)  ft- Em--.). 

'  «=1  i=l 


Remarquons  que  les  3’*"*'  et  4’*”“'  termes  du  second  membre  de  cette  expres¬ 
sion  sont  K^-adaptes  et  disparaitront  done  dans  la  normalisation. 


Derivation  par  rapport  a  la  condition  initiale. 

On  considere  les  variables  qui  in'.ervienncnt  dans  nos  calculs  comme  des 
fonctions  de  Xq  (condition  initiale)  et  des  processus  {wk}  et  {jr*}. 

On  note  V»jt  =  ^{Xo,Wk,yk)- 
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On  a  que  X.  est  differentiable  par  rapport  a  Xq. 

On  veut  deriver  les  expressions  (25)  et  (26).  Definissons  alors  les  processus 


Zk 


A 


Z„k 


axo  ’ 


(Zo  =  1) 


(Zok  =  Z/c) 


et  fixons  k. 

La  derivee  de  X.  par  rapport  a  Xq  est  donnee  r^ursivement  par  I’expression 
(27). 


Lemme  3.4  Pour  n  >  k,  Z„k  verifie 


At 

Znk  =  n0+^  — + 
Zkk  =  1  • 


(27) 


Preuve 

De  (25)  vient  que: 

X.+i  =  <r  ^  [A(X.)  -  HMi)]  +  Xi  +  b{Xi)  At  +  ay/Kt  t&.+a 
done,  en  derivant  par  rapport  a  la  condition  initiale  Xo,  on  obtient: 

Z.-+1  =  (1  •  <r  ^  h'{Xi)  +  At  b\Xi))  Zi . 

Par  recurrence,  on  deduit  I’expression  de  Z,-,,+m : 

t+m-l 

n  (l+<T—k'(Xi)  +  Atb'{Xj)). 

i=i  ' 

Pour  n>  k, 

At 

Z„k  =  ]li^-h(T—h\Xi)  +  Atb\Xi))-^  . 

j=k  ^ 


On  a  la  majoration  suivante: 

Znk<(i+ci,—-m\  (28) 

£ 
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i.e. 

Znk  <  cexp{-c^^}  . 

D’autre  part,  en  derivant  (26)  par  rapport  a  Xo,  on  obtient  I’expression: 

OAq 

+—  -  M.--,)  h'{Xi^r)Zi.r 

^  «=i 

jyi{Xi.,)  -  h{Mi.,)\[Zi  -  Zi.r  -  At  6'(X._i)  Z._,] 

At 

+— i:[6(AVi)  -  6(M.-_i)]  h'iXi^r)  Zi.r  . 


Expression  asymptotique  pour  Xk  -  A/*. 

L’egalite  precedente  nous  permettra  d’ecrire  une  expression  asymptotique  pour 
A,  -  M,. 


On  obtient,  tout  d’abord,  I’egalite 

a  £ 

=  —  -  A/._,)(/«'(A,_,)  2.-1  -  h'{Xi.2)  2, -.2] 

1=1 

At  * 

)  -  M  Ai-i  )l  h'{Xi., ) 

^  <  =  1 

At  ^  d 

+—T\^{Xi-x)  -  -  — log(L,A0  , 

0-£^  C/Ao 


a/ec  Z_i  =  0  . 

On  utilise  les  developpement  de  Taylor  de  2”“^  et  I*"'  ordre,  respectivement,  de 
k  et  b  : 

h{Mi.,)  -  h{Xi_,)  =  A'(AV0(A/.-., + 
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On  obtient  I’egalite: 
Xk  -  Mk 


i)  ~  ^i-i)W{Xi^i)  Zk-i,i-i - —  Zk-i,i-i 

—h'{Xi^2)  Zk-i,i-2] 

At 

+ -  Mi-.)  ‘'(Ci.)  A'W-.) 

(T  Ai  , 

+2eh'{X,.,)  '""K'-') 

log  (LkAk)  Zk  -1.0  • 


h'{Xk.i)dXo 

Mais 

•  pour  i  >  2, 


(29) 


<T  At 


k'(Xi.,)  -  — -  V(X,.,)"  Zi.,.,.,  -  z...,(., 

=  h"(x,.2){x,.^  -  z*-,,i.j  +  i  -  Xm)*  z..,,,.2 


a  At 


112 fcX  \  ,  ^  t'/v  \tinttX  M/v  V  \2 


^  +  5  V(A',_j)  -  A,.,) 

+V(A,.j)  A''(X,.,)(A._,  -  Xi.2) }  Zu-u-t 

+[Ai  V(A,_. )  6'(X«)  (1  -  ^  A'(A,-..))  -  2^  A'(A:,.,  )’  A'CATm)!  Z,..,i-2  , 
puisque 

A'(A,.,)  -  A'(A,_,)  =  A''(A,.;)(A'...  -  A,.,)  +  i  A'”(£,._, )( A..,  -  A,..)^ 


et 


A'(Ai.,)[A'(A,..)  -  A'(A..,)] 

=  -  Xi.^f  +  A'(A,.2)lA"(Ai.2)(Ai.,  -  Am) 

+5  '.'"(&)(->fi...  -  -Vm)"]  . 


Soit 


i  AiA'(A,-.i)y(AM)(l  -  —  A'(Ai.,))  -  ^^A'(Aj.,)H'(AM) . 


31 


On  a  la  majoration 


Af 

\<f>i\<C{-j-VAt),Wi>2. 


•  Pour  i  =  1, 


h\Xi.,)  -  ~  Zk.i,i-1  -  h'{Xi.2)  Zk-X,i^2 

=  {l-(r~h\Xo))h'{Xo)Zk-t,o 
—  4>\  Zk-ifl 


4,,th'{Xo)i\-<T^h'{Xo)) 


L’egalite  (29)  devient: 
Xk-Mk 


TTV  y  -V  -  Xi.2)  h"{Xi.2)  Zk^x,i-:x 


£h'{Xk-x) 


(f>i 


1  ^ 

’’’AWIT)  1  +  cr  A(/£  h'{Xi-i)  +  A( 

+pp^(A'„  -  M„)  A'(A„){1  -  ^  ;,'(X„)  +  Ait'CCo*)) 
~  '''**> 

+27f(^  DA.-.  -  M.-.  f  A'(A.-. )  4"{£i. ) 

<Te  9  fr  A  \  rr 

~li'{Xk.x)  ^  ■ 


On  remarque  que: 

_ ti _ 

1  +  a  At/e  h'{Xi.2)  +  Af  y{Xi^2) 


+  Ai  6'(C^ . )  A'(X(.. )  I  <  C  (=^v Ai)  S  ^ ,  Vi  >  2 . 
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On  cherche  maintenant  des  majorations  dans  pour  les  differents  termes 
de  droite  dcins  (30),  a  I’exception  du  dernier. 

Voyons  d’abord  que; 

£[|(A-i.,  -  <cVAt  {E[(X,.,  -  Mi., )»!)■/»  . 

Puisque,  voir  (24), 

E[\Xi.r  -  <  C [(1  - 


ou 


At 


l-A^(l-c,  — +  1|6'||A0 


et,  en  outre,  d’apres  (28), 


<  (1  +  Cl  — -11411 


on  a  que: 

•  le  I'"’  terme  est  majore  (dans  L^)  par: 

C  \/ee\p{—c—}  +  C  s  ,  pour  A<  =  e®  . 

£ 

Cette  estimation  est  justifiee  dans  I’annexe  B. 

•  le  2’”®  terme  est  majore  (dans  L^)  par: 

^  A/  A/ 

CAt  x:[(l  -  Cl  —  +  II41I  +  C  v^l(l  +  Cl  —  -  II41I  Al) 

i=2  ^  ^ 


-{fc- 


A/  Af 

<  C  Af  (1  +  Cl  —  -  II61I  A()-‘  Ell  -  (Cl  —  -  II61I  AOT' 

£  7^  S 


At 

+C  At  y/e  ^(1  +  c/, - ||6^||  At)~' 

«=i  ^ 

<  CA<(H-ca  — -  ||6'||A0"*(it-l) 

€ 

^  C  +  c, -  II61I  Ai)-' 
+C  At  y/e - ^ 


A< 


1-(H-CA— -||6'||At) 


-1 


<  Ceexpf— c— }  +  , 

e 


puisque 


E[\Xi.i  -  Xi.2\iXi.i  -  <CAt  iE[iXi.^  - 


«+i) 


•  le  3"*®  terme  est  majore  (dans  L^)  par: 

6'e“-^'2exp{-c-}  +  Ce“  , 

£ 

puisqu’il  s’agit  de  la  some  d’un  terme  du  meme  type  du  1’®’’  terme  ici 
etudie  avec  un  terme  du  meme  type  du  2"*®,  multiplies  par  un  facteur 

d’ordre  — . 
e 

•  le  4"*®  terme  est  majore  (dans  L^)  par: 

k  Af  Af 

CpYH^  -  «  —  +  l|6'li  A0‘-‘  +  C  VJ1(1  +  CA  —  -  ||y|| 

t=2  ^  ^ 

<  C  (£  V  £"-‘)  exp{-c  +  C  (£“-*/2  V  £3/2) 

•  le  5”*®  et  6”**  termes  sont  majores  (dans  L^)  par: 

C  (1  +  CA  —  -  11611  <  Cexp{~c^} 

£  £ 

•  le  7"*®  terme  est  majore  (dans  L^)  par: 

c  —  E{(1  -  CA  —  +  ||4'||  A()“<'-'>  +  C£}(1  +  CA  —  -  ||6'||  A()-l‘'-‘l 

£  ,_2  £  £ 

<  Cexp{— c— }  +  C£ 

£ 

done  on  etablit  I’estimation 

+  ^  log  (LkXk)  2a-., 0  =  V  i) ,  (31) 

au  sens  ou 

£[|A:a-A/a+^^ ^ log (1,aA>) 2a-i,,|1  <  Cexp(-c^^}+C(A“-‘'VA) . 
Vu  que  h'  est  bornec  superieurement  et  inferieurement,  on  peut  ecrire  que 
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On  utilise  la  decomposition 


oil  le  2”*^  terme  est  un  produit  de  deux  facteurs  d’ordre  0(y/e)  dans  (cf. 
proposition  3.2)  et,  par  consequent,  d’ordre  0(e)  dans 

Prenons  I’esperance  conditionnelle  par  rapport  a  V"*'. 

On  va  montrer  que: 

£:|j|-l0g(i»Ai)Zi.,,„n=O(l)  (32) 

et  alors  on  a 

Xfc -Mfc  =  V£), 

au  sens  ou 

E[\Xk  -  Mk\]  <  Cexp{-cj}  +  C  V  e)  . 

On  utilise  le  lemme  3.5  pour  prouver  ce  resultat. 


Lemme  3.5  On  a  Vestimation 


au  sens  oil 


£(j|-lo8(i^AOZ>-,.ol?‘)  =  0(l), 


£'(j|-Iog(itAi)Zi-i,ol?*l 


<C. 


Preuve 

On  rappelle  d’abord  que: 

S\  g  ei  g  et  ~  sont  des  fonctions  integrables  par  rapport  a  la 
ox 

mesure  de  Lebesgue,  alors 


En  particulier, 


(33) 


soit  Ip  une  v.a.  differentiable  par  rapport  a  Xq  telle  que: 


filial  +  1^  ^(Xo)\  +  &|]  <  oo  , 

Po  OAo 


drp 


(34) 


soit 


g(x)  -=  po{x)  tj}(x,  {«)*},  {j/it})  , 

oil  ip{x,  {tbfc},  {y/c})  est  une  fonction  differentiable  par  rapport  a  x, 
alors 

rJ 

(35) 


oXq  Po 


On  applique  la  formule  de  changement  de  probabilite, 

a  B(^log(i»Aj)Zt.,,„ijAt|J'‘) 

£l;s|-(iiA»)Z4.,.o|?‘) 

_  uAq _ 


Soit 


Puisque 


£li»A»|?*l 


k  =  £(5|-(i»At)Z».,.„|J-‘] . 


d  0  d 

j^{LkAkZk-i,o)  =  gjr{Lk^k)  Zk-1,0  +  {UAk)  -^Zk-i.o 


on  a  que 


h  =  E[—{LkhkZk-x,o)\Y’^]  -  E\{Lk^k)  ^Zk-^,o\y’^\  . 

OA.Q  (/Aq 

Etant  donne  que  (LjtAjt)  verifie  la  condition  d’integrabilite  (34)  on 

applique  (35)  pour  obtenir: 

£(T^(itA»Zi_,,o)|?‘)  =  -E[(Lu^^)Zt-lAx„)\V'’\ 

oAq  Po 

et  alors 

k  =  -E[  ^{X„){LM  Zi.,.„|f  *)  -  EKLtAi)-^  Zi-ifll?*!  • 

Po  vAo 
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On  a  done 


£[(itAj)  ^(Xo)  £|(i»A»)  Zi.,,o|?*l 

_ _ go _ OAq  _ 

£(i*Ai|?*l  «(i4At|F*l 

Po  Cf^O 

Puisqu’on  avait  suppose  que 

/  |— (a;)rPo(a:)dx  <  oo  ,  po  E  (36) 

>/  Po 

et,  de  (28),  on  a  la  majoration 

Z*-,,0  <  {1+c,  — -||4'||A()-'‘-*> 

£ 

<  Cexp{— 

e 

il  vient  que  le  1’®''  terme  est  d’ordre  Cexp{— c^^}.  Alors,  il  nous  suffit  de 
demontrer  que 

=  C(I) .  (37) 

Or,  de  (27),  on  peut  obtenir  la  derivee  partielle  de  Zno  par  rapport  a  Xq: 

8  At  At 

^Zno  =  -Y^{a  —  h"{Xi)  +  Atb''{Xi)){l  +  <7  —  h\Xi)  +  Atb'{Xi))-'^ 

OAo  ,_0  £  £ 

At 

.  l[{l  +  (T~h'iXj)  +  Atb'{Xj))-'^ 
i-i  ^ 

et 

done 

Ij^rZ^I  <  (-11  ^h"\\  +  ||6”|(  Ai)(l  +  -  -  lli'll  A()-’  t:(l  +  ^  -  lli'll  A()-‘ 

dXo  £  ^  ^ 
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At  At  (1  +  c,.— -||5'i|A<)-* 

<  C -  (1 +  Ill'll  ".. 


l-{l  +  CH  —  -m\At) 


-1 


<  c-, 

ce  qui  prouve  (37)  et  tennine  la  demonstration  du  lemme. 


Remarque  3.6  On  peut  demontrer  un  lemme  plus  general  que  celui-la,  lequel 
peut  etre  enonce  de  la  fajon  suivante: 


Lemme  3.7  Soit  {7^}^.  un  processus  adapts  (dependant  de  e),  differentiable  pa^ 
rapport  a  Xq  (Vi  =  0, 1, ,  K). 

dF 

Si  les  moments  de  {  }  sont  finis  et  si  =  0(£'^),  alors,  pour  un  certair 

9  >  0,  “  ' 


(38) 


3.1.3  Estimation  de  I’erreur  de  filtrage 

Enfin,  pour  obtenir  une  estimation  de  I’erreur  commise  dans  une  etape  de 
filtrage,  il  nous  suffit  de  remarquer  qu’une  approximation  du  filtre  optimal 
^k-\  —  F\Xif_\\V  ]  est  donnee  par  le  schema  (20),  etant  donne  que 

A^.-AV,  =  E\X,\?’^]-E[X,.,\?^] 

=  E{ib{Xk-i)At  +  (Ty/Atwk)\Y’^] 

=  At  E[biXk.r)\?'‘] , 

i.e. 

|A'jk-i  -  Xk\  <cAt .  (39) 

Corollaire  3.8  Le  schema  (20)  verifie: 

E^(|A^jt  -  Afjtl]  <  C exp{-c— }  +  C  w  e) . 

e 
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Remarque  3.9  Dans  le  C2is  h  lineaire,  on  voit  immediatement  qu’il  ne  reste  plus 
que  les  4’"*,  5”*®  et  8”*®  termes  dans  I’expression  (30),  tous  les  autres  disparaissent 
et,  dans  le  lemme  3.5  au  lieu  de  0(1)  on  a  0(e).  On  obtient  ainsi  que,  si  h  est 
lineaire, 

Xk-Mk  =  0(£“-^/2  V  £^/*) , 


i.e. 


E[\Xk  -  MkW  <  Cexp{-c-}  +  C  (£“-^/*  V  . 

£ 


3.2  Le  cas  a  =  1 

On  considere  maintenant  le  cas  ou  le  parametre  vaut  £,  i.e.  a  =  1  dcins 
I’hypothese  (H5). 

Le  systeme  considere  est  done  decrit  par  les  equations: 


=  Xk-^  K^k)  £  +  O’  \/£  Wk+l  )  ^0  =  ^ 

Vk+i  =  h(Xk)  +  \/evk+i. 


(40) 


En  plus  des  hypotheses  enoncees  en  debut  de  cette  section  dediee  au  cas  non 
lin^ire,  on  suppose  que: 

(H6)  \\(Th'\\<2cH. 

Remarque  3.10  Une  telle  restriction  signifie  que  notre  etude  n’est  valable  que 
si  la  fonction  h  “n’est  pas  trop  eloignee”  d’une  fonction  lineaire. 

Remarque  3.11  En  fait  on  peut  considerer  une  hypothese  legerement  moins 
restrictive,  mais  moins  intuitive  quand  a  I’interpretation.  II  suffit,  dans  le  lemme 
3.13,  d’avoir  le  soin  de  chercher  une  minoration  du  gain  un  peu  plus  fine. 

Le  filtre  approche 

On  aimerait  pouvoir  construire  un  filtre  approchee  unidimensionnel,  pour  ce 
probleme.  Cependant  cette  tache  se  revcle  delicate,  ce  qui  nous  ramene  a  la 
recherche  d’un  filtre  du  genre  du  filtre  de  Kalman  etendu. 

On  considere  le  filtre  bidimensionnel  suivant: 


Mk+i  =Mk  +  b{Mk)  e  +  $k  iVk+i  -  HMk)) .  A/o  =  mo  , 


(41) 
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et 


(42) 


avec  mo  quelconque  et  Oq  tel  que  <t  0o  >  0 . 

Pour  b  et  h  lin^ires,  ce  filtre  correspond  a  prendre  le  filtre  de  Kalman  - 
Bucy  et  negliger  les  termes  en  b  At  dans  I’expression  du  gain. 

On  cherche  d’abord  a  estimer  Xk  —  Mk  et  ensuite  Xk  —  Mk-  Les  methodes 
utilisees  pour  obtenir  ces  estimations  sont  essentiellement  les  memes  que  dans 
la  section  precedente.  La  meme  consideration  sur  I’erreur  de  filtrage  qui  a  ete 
faite  dans  le  cas  precedent  est  vraie,  i.e.,  avec  les  notations  du  paragraphe  3.1.3, 
on  s’interesse  d’abord  a  I’estimation  de  Xk  —  Mk  et  puis  on  utilise  le  fait  que 

\Xk-x-Xk\<cK 

pour  deduire  le  resultat  recherche,  I’estimation  de  I’erreur  de  filtrage. 

3.2.1  Estimation  de  Xk  -  Mk 
On  etablit  la  proposition  suivante: 

Proposition  3.12  On  a  I’estimation 

E[{Xk  -  Mk)^]  <  Cexp{-c-}  +  Ce. 

e 

Pour  demontrer  cette  proposition  on  fera  appel  a  un  lenune  (lemme  3.13) 
qu’on  presentera  par  la  suite. 

Preuve 

Etant  donnees  les  hypotheses  de  regularite  sur  les  fonctions  h  et  6,  on  peut  ecrire 
leurs  developpements  de  Taylor  du  1'"  ordre: 

h{Xk)  -  HMk)  =  h'i^^){Xk  -  Mk) 

et 

b{Xk)  -  b{Mk)  =  6'(C*')(AT  -  Mk) . 

On  obtient  alors  I’cxpression  recursive: 

-  Mk  yx  =  (Xk  -  Mk)  +  (6(A't)  -  b{Mk)]  e  -  dk{hiXk)  -  hiMk)] 

+\/ei<Tivk+i  -OkVk+i) 

=  [l-hh'i^^)  +  b'((:^)e](Xk-Mk) 

+\/e{(TWk^i  -  OkVk+i). 
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D’apres  le  lemme  3.13  ci-dessous,  pour  eo  “assez  petit”,  on  a 

Ve  <  eo ,  |1  -  4  h\ek  )  +  b'iC^)e\  <  1  -  A  <  1 . 

D’autre  part,  puisque  {dk)k  est  une  suite  bornee,  on  a  la  majoration 
£[(<710*4.1  -  Ok  Vit4.i)^]  =  <7^  +  E[6l]  +  2  <7  E[0k  Wk+i  t;*4.i]  <  C 
et 

E  {(<Ta;*+i  -  Ok  vk+i)  [(A'*  -  M*)  +  [6(X*)  -  6(A/*)]  £  -  dk[h{Xk)  -  h{Mk)]] }  =  0 

done 

E[{Xk+i-Mk+xf]  <  il-AfE[{Xk-Mkf]  +  Ce 

<  (l-/l)2{*+*)£[(Xo-Mon  +  C'^, 

d’ou  la  proposition. 

■ 

II  nous  reste  a  etablir  le  lemme 
Lemme  3.13  II  existe  eo  >  0  tel  que,  pour  un  certain  >1  >  0, 

V£  <  £o ,  |1  -  h  A'Kf)  +  6'(C*)£  I  <  I  -  /I . 


Preuve 

On  utilise  le  fait  que  5*_i/<7  >  0  et  I’hypothese  (H3’).  De  (42)  il  vient  que 


h  =  cTh'{Mk)[h-i/<T  +  <rn'{Mk.x)] 

a  <7^  h'^{Mk)[0k-il(r  +  <7*  h'{Mk-i)]  +  h'{Mk-\) 

1  ^  J_ 

^  crh'{Mk)  ~  Ch 

et  alors  on  a  que 

h  _ a^h'{Mk)h'{Mk-x) _ 

<7  (T'^h'^{Mk)[\l{<jh'{Mk-i))A<Th'{Mk-\)]+cfh'{Mk-i) 

-  2||<7A1P  +  ||<7/i'||-' 

On  obtient  les  bornes  inferieure  et  superieure  suivantes,  pour  Sk  h'{^k)‘ 


2|KIP  +  ||<7/<'||^ 


Mil 

Ch 


d’ou  on  deduit  que: 
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•  P°“^‘<2||<TV||=  +  ||aMl 

on  a  la  majoration 


-et£<c.  =  Ti-i,( 


r* 


4  -  -^l)* 


l-hh\i^)  +  b\C^)e<l- 


4 


•  pour  /I2  <  —  et  £  <  £2  = 

Ch  Ch  lift'll 


2lK1P  +  ||aA'|l 

c-AiCh 


ll6'ir2l|a/i'lP  +  lla/i'll 

7  +  11^*11  £  1  ~  -^1  > 


,  on  a  la  minoration 

Ike'll 


1  -  h  +  6'(Cf )  £  >  1  -  i!— -  ||5'|1  £  >  ,42  -  1  >  -1 


Ch 


II  suffit  de  considerer 


A  =  min{Ai,  A2)  et  £0  =  mm  (£i,£2). 


3.2.2  Estimation  de  —  Mk 

Le  th«Dreme  suivant  nous  donne  une  estimation  sur  I’erreur  de  prediction  Xk  — 

Mk. 

Theoreme  3.14  On  a  I’estimation 

E[\Xk  -  Mk\]  <  Cexp{-c^-}  +Ce. 

£ 


Preuve 

La  methode  suivie  pour  demontrer  ce  theoreme  est  essentiellement  la  meme 
que  dans  le  c<is  precedent. 


Changement  de  probabilites 

(a)  On  considere  la  probabilite  P  precedemment  definie,  avec  At  =  e. 

(b)  On  considere  la  nouvelle  probabilite  P  definie  par: 
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avec 


i=l  ‘'•-1 

,=l  *'i-l 


m  =  m  — 7=5 —  Ih(Xi^i)  -  A(A4_i)] 


Alors 


A,  =  exp(— 

,=1  ^t-1 

E  a?- )  - '‘(M-i)IM  ■ 

^  ^  .=1  ^i~l 

Sous  P,  Wk  et  Ujt  sont  des  bruits  blancs  gaussieos  indei)endants  et  Xk 

v/e 

verifie  maintenant  I’equation 

Xk+i  =  -5-  m^k)  —  +  Xk  +  f>iXk)  £  +  o’v/e  ii)A:+i  •  (43) 

f'k 

Designant  toujours  par  S  le  logarithme  de  Lk  A*,  puisque 

<Tv/e  CTy/e  <7 

—7:^  (/^(AV,)  -  h{Mi,i)]  +  %  [j/,-  -  , 

y/eOi^i  CTy/e 


S  =  --f  ji-WA',-..)-A(A/i-,)|((X,-A^.)-(A,-,-M-.,)l 
£  .=1  ‘'.-1 

+  12  r-  WXi-, )  -  A(M-i  )][i<A'.-, )  -  KMi-,  )1 

.=1  ‘'i-l 

+-  El-:;  A'lAi-.)  +  A'(Ai-i)  '•(Mi.,)! 

^  <=J 

y,_i  t  ,=i 

Le  dernier  terme  dans  cette  egalite  est  K^-adapte  et  disparaitra  done  dans  la 
normalisation. 
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Derivation  par  rapport  a  la  condition  initiale. 

Avec  les  memes  notations  utilisees  precedemment  on  obtient, 
•  en  derivant  I’expression  (43), 


oil  -j-  h!{Xx']  >  cl  (voir  la  preuve  du  lemme  3.13).  On  a  done  la  majoration 

Znl..<(l+cl-i|6'l|  £)-<«-«  . 

•  en  derivant  I’expression  (44), 


d 


1  A  1 


dXo 


-log(I,A0  =  —^  —  {X,-M.)k'iX,^r)Zi.r 

e  1=1  ^<-1 

+1 E  -  «-■ )  '>'(^'.-1 )  ^.-1 

^  .=1  ‘'<-1 


+  h'iX,.^)  Z._1 

1=1  ^>-1 

“7  E  (MA'.-i)  -  M/V/._0]  {j{Zr-  2^.-1  -  b\Xi.,)e  Z,_,] 


_  A2 


Expression  asymptotique  pour  Xk  -  Mk- 
On  obtient  I’expression  suivante; 

At  -  Mk 

=  T7;?^DA._. 


h'{Xk-^)U' 


Oi-2 


Zk-l,i-2  } 


■£^log(i,.Ai.)Z*-i,o. 


ax„ 
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Mais 

•  pour  i  >  2, 


Vi~2 

-  h  (X,_i)  [^—  -  h  (AT.-i))  y{Xi_2)  £  Zk-\,i-2  +  <i>i  Zk-i,i-2  , 


4i  =  A'(A',.,)  (i  -  A'(AV.))  (1  +  ^  . 

C/,_2  Bi_2 

Etant  donne  que  0^  verifie 

-  A'(M._,)]  (1  +  ^  A'(M,_2)]  ~  =  0 

'-1  t'i-2  y,_2 

et,  en  utilisant  les  developpemenis  de  Taylor  jusqu’a  I’ordre  2  de  la  fonction  h> 
on  a  ’ 

2 

+  '•'(M.-,)) }  (A,.,  - 

■"'  + crb 

•  pour  t  =  1  on  a  tout  simplement 

<*.  =  '■'(-’fo)  (i  -  (.'{A'o)) . 

I'O 

L’expression  de  Xu  -  JVA  devient: 

Xk  —  Mk 


Ok-t  X 


h'(Xk-t)  ~  -  Af,_2) 


^  ^.-1  CX^  +  ^W-i)l[l  +  -f-  h'{Mi_2)] } 

Vi-2 
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+2  9(xt';) 

+77nfS  '  -  «<-.)  )  1 )1  '''(^■->) 

h{Xk-i)  ,_2  Wt-i 

+77^^'  -  jw^-.)  t'K.^.) 

/i(Ai._ij  ,=2  ^i-l 

■^-J^  (A'o  -  Mo)  {i  -  A'(A'o)  +  £  6'(tf )  i  }  V(Ao)  Z*-,,o 

/i(Afc_i)  t/o  po 


Comme  dans  le  cas  precedent  on  cherche  des  majorations  dans  pour  les 
difFerents  termes  a  droite  dans  cette  expression,  a  I’exception  du  dernier. 

On  a  la  majoration  (voir  le  lemme  3.13  et  la  proposition  3.12) 

£{1(A,_, -M._,)(X._2-M._2)i) 

<  eill  -  h'((^)  +  6'(Cf)  £  I  (A',.,  -  M,.,)>) 

+  y/e  E[  1  <7  u;,_i  —  0,_2  Pi_l|  IX, _2  —  M,_2|  ) 

<  (1  -  A)  £((X,_2  -  M-2)1  +  CVe  iE[{Xi.2  - 

<  (1  -  Ay'-'^E[iXo  -  Mof]+Cy/^il  -  Ay-\E[{Xo  -  A/o)'])’/'  +  Cs, 

oil  la  constante  A  dans  le  lemme  3.13  pent  etre  choisie  telle  que  1  —  >1  > 
(l+cj-||6'||  £)-■«. 

On  a  alors  que: 

•  le  r"  terme  est  majore  (dans  L^)  par: 

C  ^  {  C  (1  -  +  Ce  +  Cy/eil-  Ay-^  }  (1  +  -  ||6'|| 

1=3 

<  C exp{-c—  }  +  Ce 
£ 

•  le  2"**  terme  est  majore  (dans  L‘)  par: 

C  Ol  -  +  Ce](l  +cl-  ||6'lle)-^''-'+')  <  Cexp  {-c-  }  +  Ce 


(45) 
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•  le  3*"'  termc  est  majore  (dans  L^)  par: 

c  E[(l  -  +  Ce](l  +  4  -  life'll <  Cexpf-C^  }  +  Ce 

t=i  ® 

•  le  4”*'  terme  est  majore  (dans  L^)  par: 

CeX^[(l->l)('-'4C'>/F|(l+4-||fe'||£)-(*-+')  <  Ceexpl-c^l+Ce^/^ 

t=2  ^ 

•  le  S’"*  terme  est  majore  (dans  Z,')  par: 

C  e  exp  {  — c  }  +  C" 

•  les  S'"*  et  7”*'  termes,  representant  la  contribuition  de  I’erreur  a  I’instant 
initial,  sont  majores  (dans  L^)  par: 

C(1  +  4  -  llfe'lkr^*"’^  <  C7exp{-c^}. 

Finalement,  on  a  que 

E[\Xt- Mu  +  xn¥—\  ^  l«8  (it  At)  2t-i,o  1 1  <  C  exp{  -c  )  +  C  £ . 

h'{Xk~i)  oXo  e 

Comme  dans  la  preuve  du  theoreme  3.3,  prenons  I’esperance  conditionnelle  par 
rapport  a  F*  dans  I’expression  de  Xk  —  Aik- 
Le  lemme  qui  suit  complete  la  demonstration. 


Lemme  3.15  On  a  I ’estimation 

E[^  log(ii.AO  4-..0 1  f*]  =  0(1) . 

Preuve  La  preuve  de  ce  lemme  suit  les  memes  pas  de  celle  du  cas  a  >  1.  II 
suffit  de  remarquer  que  Zno  est  donne  par 

Zno  =  X 

1=0  "'■ 
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et  sa  derivee  par  rapport  a  Xq  est  done 


-T,(t  =)(i  +  X 

ss  •’i 

■  na  +  r^TO  +  ^TO^)''- 

]=i  ^3 


d’ou  la  majoration; 


d 


!^2„o|<C(l  +  4-||61|e)-i^(l+c; 


'dXo 


e)-^  <C. 


t=i 


Remarque  3,16  Dans  le  cas  d’une  fonction  h  lineaire,  on  pent  verifier  que  les 
termes  dans  I’expression  (45)  disparaissent  a  I’exception  des  4”’®,  5”*®,  7”*®  et  S’"® 
termes.  On  obtient  ainsi  une  estimation  d’ordre  0(e)  au  lieu  de  I’ordre  0(1) 
dans  le  lemme  3.15.  On  conclu  alors  que,  si  h  est  une  fonction  lineaire,  on  a 

A'fc-iV4  =  0(£"/^), 


i.e. 


E[\Xk  -  MkW  <  Cexp{-c^}  +  . 

£ 


3.2.3  Estimation  de  I’erreur  de  filtrage 

De  meme  que  dans  le  paragraphe  3.1.3,  on  obtient  immediatement  une  estima¬ 
tion  de  I’erreur  de  filtrage  a  partir  de  I’estimation  de  I’erreur  Xk  —  Mk- 

Corollaire  3.17  On  a  I'eslimation 

El\Xk  -  Mk\]  <  Cexp{-c^}  +  Ce. 


3.3  Le  cas  a  <  1 

Dans  le  cas  oil  At  =  avec  a  <  1,  il  est  bien  connu  que  les  observations 
prises  dans  le  passe  ne  sont  pas  de  grande  utilite  dans  le  calcul  de  I’estimation 
courrante.  On  fait  le  raisonnement  suivant: 
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D’apres  I’equation  d ’observation  on  a 


yk-hiXk)  = 


Wk 


est  un  bruit  blanc  gausssien  de  variance  done  il  est  bien  connu  que,  pour 
Co  >  on  a  la  majoration 


Pi  \yk-hiXk)\>co)< 


,  1,  Co 

I  I AVrv  J  f  -  *  * 

Co\/5x 


Supposons  que  la  fonction  h  est  inversible  et  que  h~^  est  continue  et  a  derivee 
bornee. 

Prenons 

Mk  =  h-^(yk),  (46) 

i.e.  notre  filtre  approche  consiste  tout  simplement  a  prendre  1 ’image  inverse  par 
h  de  I’observation. 

Alors 

done,  pour  ci  >  on  a  que 


En  termes  d’esperances  on  a  tout  simplement  que 


B[(Xk  -  Mk)^}  <  |i(A-‘)'lP  m(Xk)  -  yk)^j  =  ll(h-^)r 


At 


En  un  mot,  le  filtre  definit  par  (46)  est  deja  un  filtre  “performant” . 

Quand  au  resultat  concernant  I’estimation  de  I’erreur  par  rapport  au  filtre 
optimal,  il  semble  difficile  de  trouver  le  bon  changement  de  probabilite  qui  con¬ 
duit  a  I’expression  asymptotique  pour  Xk  —  Mk  dont  les  termes  soient  facilement 
estimes. 
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3.4  Recapitulatif 

Le  tableau  suivant  nous  donne  le  filtre  approche  a  prendre  selon  la  valeur  de  a 
ainsi  que  I’estimation  sur  son  efficacite. 


filtre  approche 

estimation  de  Verreur 

E[{Xu  -  Mu?] 

Q  <  1 

Mk  -  h 

a  =  1 

Mk  =  Mk_i +6(Mt_i)c  +  ^t(yjt  - /»(Mt_i)) 

.  h'{hh)[6k-i+<r^-h!{Mi.x)] 

ce 

ce 

a  >  1 

At 

Mu  =  Mu.i^b{Mu.{)Ai  +  a—{yu-h{Mu.x)) 

ce 

^(ga-l/2  vg) 
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4  Application  numerique 

Supposons  qu’on  veut  resoudre  le  probleme  de  filtrage  suivant: 

A  chaque  instant  t  G  [0,  T]  on  voudrait  estimer  le  processus  X,  donne  par 
I’equation  difFerentielle  stochastique  ci  dessous,  etant  connue  I’observation  Y, 
jusqu’a  I’instant  t 

(  dXt  =  —Xt  dt  +  a  dw\ 

1  dYt  =  {X^-\-Xt)dt-^edw1, 

oil  et  sont  des  processus  de  Wiener  independauts  et  Xo  est  une  variable 
aldatoire  gaussienne  de  loi  A/*(J?o»Po)»  avec 

•  <T  =  1.  •  Xo  =  2. 

•  £  =  1,0.1,0.01,0.001,0.0001  •  po=l. 

La  resolution  numerique  de  ce  probleme  peut  etre  envisagee  de  deux  fagons 
distinctes:  on  peut  discretiser  I’equation  de  ZaJcai  ou  utiliser  un  filtre  approche 
pour  calculer  une  approximation  de  la  solution.  Dans  notre  Ceis,  on  commencera 
par  discretiser  le  systeme  continu  et  ainsi  obtenir  un  systeme  discret.  Pour  ce 
systeme  discret,  on  discutera  I’application  du  filtre  approche  propose  dans  la 
section  3.1  (bien  que  I’hypoth^e  (H3’)  ne  soit  pas  verifiee  dans  ce  cas).  On 
remarque  qu’il  s’agit  d’un  filtre  de  dimension  1,  au  sens  ou  on  a  une  seule 
expression  recursive  a  evaluer.  Get  estimateur  nous  permettra  d’obtenir  une 
approximation  du  filtre  optimal  facile  a  calculer.  Par  ailleurs,  on  avait  demon- 
tre  des  resultats  cisymptotiques  garantissant  la  qualite  du  filtre  (voir  la  section 
3.1).  Les  pas  de  discretisation  a  considerer  varieront  selon  la  valeur  de  e  et  on 
regardera  I’intervalle  de  temps  (0, 1]. 

Dans  la  suite  on  se  propose  de  verifier  numeriquement  la  performance  du 
filtre  approche  par  rapport  a  I’estimation  donnee  pair  I’equation  de  Zakai.  Le 
comportement  du  filtre  par  rapport  a  la  trajectoire  simulee  est  aussi  etudie. 

Les  3  figures  qui  suivent  (figure  1,  figure  2  et  figure  3)  representent  I’erreur 
quadratique  moyenne  (eqm),  obtenue  sur  100  simulations  de  trajectoires  inde- 
pendantes.  Pour  des  raisons  evidentes  on  utilise  des  echelles  logarithmiques  sur 
les  deux  axes. 

On  voit  immediatement  que,  quand  e  est  petit  (e  =  0.1, 0.01,  •  •  •),  la  courbe 
logarithmique  s’approche  d’une  droite,  i.e. 

On  peut  done  constatcr  que  la  proposition  3.2  et  le  theoreme  3.1.2,  enonces 
precedemment,  sont  verifies  numeriquement. 
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iqti(er.qut.a.) 


O.lOOE-04 


p  •  n  t  6  < 

gamma 


-  0.493B-*^00 

«  0.283B4.00 


cpilloa 


Figure  1:  I’erreur  quadratique  moyenne  entre  le  filtre  approche  et  la  trajectoire 
simulee,  en  fonction  de  e. 


0.251E-03 


0.100E>04 


p«Ate  «  O.lOOB'^^Ot 
gamm*  •  0.3ttB<^00 


epiiloR 


Figure  2:  I’erreur  quadratique  moyenne  entre  le  filtre  approche  et  I’estimation 
donnee  par  I’equation  de  Zakai,  en  fonction  de  e. 


fftmm*  ■  0.276B4-00 


epiilon 


Figure  3: 1’crreur  quadratique  moyenne  entre  I’estimation  donnee  par  I’equation 
de  Zakai  et  la  trajectoire  simulee,  en  fonction  de  e. 


Dans  la  figure  4  on  represente,  pour  e  =  0.1,  Devolution  de  I’erreur  quadra- 
tique  moyenne  au  cours  du  temps.  Le  pas  de  discretisation  utilise  a  ete  = 
0.01. 


temps 


Figure  4:  I’erreur  quadratique  moyenne  entre  la  trajectoire  simulee  et  les  esti¬ 
mations  et  celles-ci  entre  elles,  au  cours  du  temps. 

Le  trait  le  plus  epais  correspond  a  I’erreur  entre  les  deux  filtres  et  les  deux 
autres  traits  les  erreurs  de  ces  memes  filtres  par  rapport  a  la  trajectoire  simulee, 
le  trait  le  plus  fin  etant  celui  qui  correspond  a  I’ern'm  relative  au  filtre  optimal. 
La  ressemblance  entre  le  comportement  du  filtre  approche  et  du  filtre  optimal 
est  evidente. 
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Finalement,  dans  la  derniere  figure  (figure  5),  on  represente  une  seule  trajec- 
toire  et  le  filtre  approche,  ainsi  que  I’estimation  obtenue  par  resolution  numeri- 
que  de  I’equation  de  Zakai. 


tempi 


Figure  5:  revolution  des  filtres  approche  et  optimal  sur  une  trajectoire  Simula. 

On  voit  facilement  que  les  filtres  (en  pointille  sur  la  figure)  evoluent  de 
maniere  semblable  et  qu’il  devient  difficile  de  distinguer  les  deux. 

Pour  le  calcu)  de  la  solution  de  i’equation  de  Zakai  on  a  utilise  un  pro¬ 
gramme  FORTRAN  genere  a  cet  effet  par  le  logiciel  ZPB,  travail  en  cours  de 
developpement  au  sein  du  projet  MEFISTO,  a  I’lNRIA  -  centre  de  Sophia  - 
Antipolis. 
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Dans  cette  etude  cisymptotique  de  qualite  des  filtres,  deux  difficultes  im- 
portantes  se  pr&entent,  dans  la  mise  en  oeuvre  des  algorithmes.  EfFectivement, 
quand  e  -+  0,  la  densite  conditionnelle  solution  du  probleme  de  filtrage  s’appro- 
che  d’une  gaussienne  ayant  un  pic  “tr^  prononce”,  ce  qui  nous  a  oblige  a  utiliser 
un  pas  de  discretisation  en  espace  tres  petit  dans  I’algorithme  de  resolution 
de  I’equation  de  Zakai.  Cette  difficulte  conduira  a  la  recherche  d’algorithmes 
plus  performants  pour  la  discretisation  de  I’equation  de  Zakai,  lorsque  le  bruit 
d’observation  est  petit. 

Les  essais  numeriques  ici  pr^entes  ont  ete  realist  sur  un  mini  supercalcu- 
lateur  vectoriel  CONVEX  C210.  Le  temps  de  calcul  necessaire  et  I’espace  de 
memoire  requise  nous  ont  pose  des  contraintes  pour  la  resolution  numerique  du 
probleme  de  filtrage.  D’abord,  dans  la  discretisation  du  systeme  continu  par 
rapport  au  temps,  on  est  oblige  de  prendre  des  pas  de  temps  de  plus  en  plus 
petits,  quand  t  devient  lui  aussi  petit.  Le  pas  de  temps  le  plus  petit  considere 
dans  cet  exemple  est  de  I’ordre  de  10"®.  D’autre  part,  il  y  a  le  fait  qu’on  vient 
de  citer  que  le  nombre  de  points  nec&saires  dans  la  discretisation  en  espace 
de  I’equation  de  Zakai  est  tr^  eleve.  Dans  notre  exemple  il  a  ete  de  I’ordre  de 
2  X  10®  pour  les  valeurs  plus  petites  de  t. 

Pour  la  resolution  de  ce  probleme  on  aurait  pu  aussi  bien  utiliser  d’autres 
filtres  approches,  par  exemple  le  filtre  de  Kalman  etendu  ou  encore  le  filtre  de 
Katzur-Picard,  une  fois  discretise  (cf.  [Milheiro]). 
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A  Version  discrete  du  theoreme  de  Girsanov 

Considerons  une  filtration  ,  k  E  {0, 1,  •  •  •  /’if}. 

Soient: 

•  un  ^/j-bruit  blanc  gaussien  pour  la  probabilite  P. 

•  un  processus  //,._i-mesurablc  (i.e.  previsible)  tel  que: 

]^!v?,P<oo  (P  p.s.) 


{Zk}  le  processus  defini  par: 


k  ,  k 


Zi,  =  cxp{X^  -  -  X]  b.p}  ,  A-  >  1 

1=1  “  t=i 

Zq  —  I  . 


Alors, 


1.  {Zk}  est  une  martingale  discrete. 


Preuve 


P[^k\Pk-i]  =  f^{Zk.i^\Z^k-^] 

''k-i 


=  4-./^lexp{yr.a-r-2M^P}l^^-.]« 

puis(|ue  Zk-i  est  /\_i-m('surable. 

=  4-1. 

prisque  {ie^}  esl  un  bruit  blanc  gaussien  et  .pk  est 
1  -  iP'-surai)le. 


et  Zk  esl  Pk  intcgrable. 


On  considere  la  probabilite  P  definie  par: 

(lP{u;]  =  Z,^{uj)(IP{i^)  . 
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Remarque  A.l 

Du  fait  que  {Zk}  est  une  martingale,  vient  que: 

ii.\  -7 

2.  (Version  discrete  du  Theoreme  de  Girsanov) 

Pour  la  probabilite  P  definie  dans  (47),  le  processus  {tSjt}  defini  par: 


Wk  =  Wk-  (Pk 


est  un  bruit  blanc  gaussien  . 


Preuve 

On  veut  demontrer  que: 


£(exp  {\wk)\:FK-r]  =  exp  (y )  ,  VA  €  HI 


£(exp(A%.--)iJ-,_,)  =  i  ,  VAgE. 


£[cxp(Au)fc  -  —  )ZA-|£fc_i) 

/':[exp(AtL'j:  - -;~)|£fc_i]  =  - -  ^ -  (formule  de  Bayes) 

/'[exp(Atc, -A<,p,.-^)Z,l£fc_,l 


puisque  Zk  est  une  martingale  et  d’apres  la 
remarque  ci  dessus. 

A^  1 

=  £(exp [Xwk  -  Xp>k  -  y ) exp  {^pkWk  -  :^pl)\^k-i] 

=  £{exp((A  +  ipkpok  -  :^(A  +  <,5fc)^)|£jt_i} 

=  1  ,  de  rneme  qu’en  1.  . 
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B  Estimation  du  premier  terme  dans  I’expres- 
sion  (30) 

Dans  le  contexte  du  paragraphe  sur  I’expession  asymptotique  de  Xk  —  Mk,  dans 
la  section  3.1.2,  soil 

-  X._2)  Zk.r„-2 

1=2 

(  Sk  est  Xk-i  —  mcsurable  ) 

So  =  0. 

Notre  but  est  de  majorer  /?|Sjt|. 

On  etabiit  le  resultat  suivant: 

Lemme  B.l  On  a  I'csttmation 

E(Sl,)<Cecxp{-c‘j)  +  CtK 

Preuve 

Notation; 

Dans  la  suite.  P/,_i  de.signera  la  partie  mcsurable  de  Xk  -  Mk  i-c. 

Pk.,  =  {Xk-y-.•P.l)  +  Xl{h{Xk.^)-l>iMk-l)] 
-ZJ^[h{Xk.0-f>{Mk.^}] 

Xk  -  Mk  =  Pk-l  +  CF\/X}{wk  -  Vk)  . 

D'apres  la  defmitioii  de  S^.  ci  de.ssus.  on  a  que 

Sfc+i  =  Sk  +  (A'';.  -  /\/fc)(A'i;  —  Xk-i)  h"{Xk-i)  Zk.k-i  ■ 

•  On  commence  par  estimer  le  facteur  de  Zk.k-\  dans  cctte  expression: 

/'[(.V,.  -  Mkf(Xk  -  Xk.x?h"'^(Xk-^)] 

=  XP  E\{Xk  -  Mk?lHXk.x)h"-^(Xk-x)] 

+  XI  E[{Xk  -  Mk  ?  lot  iXk. , )  ir-iXk- ,  )j 
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+  2aAff^E[{Xk  -  Mk)H{Xk.r)wkh'%Xk.',)] 

<  C  At  E{Pl,)  +  CAP  [l  +  CE[{Xk  -  M,)"]  +  C^[|P,_,|]] 

<  CAfKl-CA^  +  lli'llAO'^  +  Ce], 

puisque 

E{PL,)  <  EHX,  -  M.n  <  (1  -  ^  +  ||6'||  A()“  +  C£  . 

•  On  estime  le  facteur  5*  {Xk  -  Mk){^k  -  h"{Xk-i): 


E[Sk  (Xk  -  Mk){Xk  -  Xk.,)h'\Xk.r)] 

=  At  E[b{Xk-i)h"{Xk-i)Sk  {Xk  -  A/,)] 

+  as/AiE[b{Xk.^)h"iXk.,)Sk{Xk  -  Mk)wk] 

<  C  At  {EiSDY^^EliXk  -  Mkf]y^^  +  CAt  {E{Sl)Y^^ 

Cl  At  A* 

<  i  —  E(S!)  +  C  A1  ((1  -  £»  —  +  lli'll  Al)“  +  Ccl  . 

Z  c  S 

On  obtient  enfin  la  majoration: 


<  (l  +  Ck^-  i|6'!|  AO-'  {(1  +  Ck  —)  E(S^) 

£  S 

At 


+CAt((l  -c  — +  ||6'i|A0"'  +  C£]} 


<  CAKl  +o,^  -||(,'||A()-^‘2i(1  +a— )-  Ili'IIAir^l 
^  .=0  ^ 


At 


[(l-Ck—-i-jjb'jjAiy<'^-UCeJ 


At 

<  CAt{l  +  Ch - \\b'\\At) 


-2 


1  +  C/,  Atf £ 


At 

+  C£At{l  +Ck - lie'll  AO”' X! 

1=0 


[(i+CAAO£-||6'||A02j 
1  +  Ch  At/e 


k  r 

E 

1=1  . 


(1  +CAA</f-||6'||A02 


Or 


(1  +  c,  —  -  ||6'||  Aty  -  (1  -f  c,  — )  >  c—  >  0 
£  e  £ 


{ckAt/£-m\Atyy 

1  +  c/,  At/e 
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(1  +  c,  -  [1  -  (C,  —  -  ||6'1|  At)^  ]^>c,—  >0 


<  C eexp{-c-^}  +  C . 
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2 

Filtrage  lineaire  par  morceaux 
d’un  systeme  en  temps  discret 
avec  petit  bruit  d’observation 


FILTRAGE  LINEAIRE  PAR  MORCEAUX  D’UN  SYSTEME  EN 
TEMPS  DISCRET  AVEC  PETIT  BRUIT  D’OBSERVATION 


DISCRETE  TIME  PIECEWISE  LINEAR 
FILTERING  WITH  SMALL  OBSERVATION  NOISE 


Resume 

Nous  nous  interessons  a  un  probleine  de  filtrage  lineaire  par  morceaux  cn  temps  discret 
avec  petit  bruit  d 'observation.  Nous  presentons  et  comparons  plusieurs  tests  pcrmettant 
de  determiner  les  intervalles  de  linearite  de  la  fonction  d’observation,  notamment  dans 
le  cas  oil  elle  est  symetiique.  Sur  chacun  de  ces  intervalles,  nous  approchons  le  filtre 
optimal  par  le  filtre  de  Kalman-Bucy  correspondant.  Comme  dans  [10],  nous  approchons 
dcs  processus  discreis  par  des  diffusions  pour  estimer  les  probabilites  d’erreur  et  les  temps 
moyens  pour  prendre  une  decision. 


Abstract 

We  are  interested  in  a  piecewise  linear  discret  time  filtering  problem  with  small  ob¬ 
servation  noise.  We  presente  and  compare  different  tests  which  enable  us  to  compute  the 
intervals  of  linearity  of  the  observation  function,  spetially  when  it  is  symetric.  Over  such 
un  interval,  we  ca,"  then  appro.xirnatc  the  optim-al  filter  by  the  corresponding  K'alman- 
13ucy  filter.  As  in  (lOj.  we  approximate  discret  time  processes  by  diffusions,  in  order  to 
estimate  de  probabilities  of  error  and  the  expected  times  for  taking  a  decision. 


1  Introduction 

On  s’interesse  au  probleme  de  filtrage  non  lineaire  suivant: 


r  dXt  =  b{Xt)dt  +  <T{Xt)dUt 
[  dYt  =  hiXt)  dt  +  edVt,  Yo  =  0, 

oil  {Xt ,  i  >  0}  est  le  processus  a  valeurs  dans  IR  non  observe  a  estimer  a 
I’instant  t  au  vu  des  observations  jusqu’a  I’instant  r  du  processus  unidimen- 
sionnel  {Yt ,  t  >  0},  e  un  parametre  “petit”  et  {f/< ,  t  >  0}  et  {Vi ,  t  >  0}  sont 
des  processus  de  Wiener  standards  independants,  a  valeurs  dans  IR. 

II  est  bien  connu  que  le  probleme  est  de  dimension  infinie,  au  sens  ou,  pour 
le  resoudre,  on  a  a  determiner  la  solution  d’une  equation  aux  derivees  par- 
tielles,  par  exemple  I’equation  de  Zakai.  Si  la  fonction  h  est  monotone,  sous 
certaines  hypotheses  generales  de  “regularite” ,  le  filtre  de  Kalman  etendu,  entre 
autres,  est  une  “bonne”  approximation  du  filtre  optimal  (cf.  [Picard],  [KBS], 
[Bensoussan],  [Ji]  et  [Milheiro]).  Le  probleme  avec  h  non  monotone,  sous  une 
certaine  “hypothese  de  detectabilite”  a  ete  traite  dans  [Fleming- Pardoux]. 

Dans  le  cadre  du  filtrage  lineaire  par  morceaux,  on  suppose  que 

•  b{x)  =  /?_xl{i<o}  +  B+xl{i>o) 

•  cr{x)  =  <T_l{x<o}  + 

•  h{x)  =  +  //+xl{x>o}- 

Si  //+//_  >  0,  i.e.  la  fonction  li  est  monotone,  pour  e  =  0,  Xt  peut  etre  par- 
faitcment  connu  et  pour  e  “petit”,  le  probleme  n’ofTre  pas  de  grandes  difficultes. 

Dans  le  cas  //+//_  <  0,  i.e.  h  non  monotone,  bien  que  h(Xt)  puisse  etre 
estime  de  fa^on  precise,  il  n’est  pas  immediat  qu’il  en  soit  de  meme  pour  Xt.  Le 
filtre  de  Kalman  etendu  est  en  general  inefficace,  le  probleme  etant  de  determiner 
le  signe  de  {A'( ,  <  >  0}.  La  determination  du  signe  n’etant  possible  que  dans 
le  cas  oil  la  variance  conditionnelle  est  petite  ('),  on  introduit  Thypothese  de 
“detectablite”  suivante  notce  (HD)  : 

■  Hlal  #  Hlal  (//£?!) 

HID)  ou 

.  //!(?!  =  IIW+  et  D+  5^  B.  {HD2) . 

'Contre-exemple:  b  =  0,li{x)  =  |j:|  ,cr  =  1. 
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Sous  (HDl),  dans  [FJP],  un  filtre  approche  a  etc  propose  et  dans  [Roubaud], 
il  est  montre  que  sous  (HD2)  un  filtre  du  meme  type  peut  etre  utilise  comme 
approximation  du  filtre  optimal. 

L’idee  est  de  construire  un  filtre  approche  a  partir  de  deux  filtres  de  Kalman- 
Bucy,  associes  respectivement  aux  deux  problemes  de  filtrage  lineaire  suivants: 


dX, 

=  B^Xi  dt (T^dUi 

(2) 

dVt 

=  H+Xtdt  +  edUt 

dXt 

=  D.Xtdt  +  (X.dUt 

(3) 

dY. 

=  H.Xtdt  +  edUf 

On  calcule  ces  deux  filtres  “on  parallele”  et  on  utilise  deux  tests. 

•  Un  premier  test  permet  de  detecter  des  intervalles  de  temps  durant  lesquels 
la  trajectoire  de  ne  passe  pas  par  0. 

•  Un  second  test  permet  de  decider  si  Xt  <  0  ou  >  0  sur  ces  intervalles, 
sous  I’hypothese  (HDl)  ou  (IID2). 

Sur  chaque  intervalle  de  monotonie,  nous  approchons  alors  le  filtre  optimal  par 
le  filtre  de  Kalman-Bucy  correspondant. 

Nous  considerons  I’intervalle  dc  temps  fini  [0,  T], 

Dans  ce  rapport  nous  nous  intcrcssons  a  la  resolution  numerique  du  pro- 
bleme  lineaire  par  morceaux  dans  la  situation  oil  la  fonction  d’observation  h 
est  non  monotone.  .N’ous  commen?ons  par  discretiser  le  systeme  continu  (1)  par 
un  schema  classique  de  discretisation  en  temps  avec  pas  At  =  e.  Le  processus 
{A'/..;i,}jtest  approximepar  {.vk}k-  par  {y/Ktuk}k, 

pai-  {(AO"'  (>Wi)Ai  -  par  {l/A-lfc- 

,\ou.s  obtenons  le  inodele  discrct  suivant: 

f  jjt+i  =  3-k  +  shixk)  + s/sa{xk)Uk 

<  ^  (4) 

(  Uk  =  l>(-^k)  +  ^k  1 

avec  les  fonctions  b.  a  et  h  definies  comme  precedemment.  On  emet  les  hypothe¬ 
ses  suivantes: 

(//I)  {uk}k  et  {vk}k  sont  des  bruits  blancs  gaussiens  standards  et  independants, 

(//2)  xo  est  une  variable  aleatoire  reelle  telle  que  E(exp{cox3})  < 
un  certain  co  >  0, 


{H3)  <  0  et  <r_cr+  ^  0 . 


Sans  reduire  la  generalite  du  probleme,  on  suppose  que  h{x)  >  0  Vx  G  IR,  i.e. 
on  suppose  que 

(//3')  <  0 ,  H+  >  0  et  cr_a+  0 . 

Nous  adaptons  I’etude  faite  dans  le  cas  continu  au  cas  discret.  Notre  but  est 
de  mettre  en  oeuvre  les  tests  permettant  de  separer  les  intervalles  de  positivite  et 
de  negativite  de  {xfc}*;  et  de  comparer  leur  performance.  Nous  nous  interessons 
principalement  a  la  situation  ou  I’hypothese  {HD2)  est  verifiee,  le  probleme 
sous  I’hypothese  (HDl)  ayant  ete  traite  dans  [FJSZ]. 

L’etude  de  revolution  de  la  loi  conditionnelle,  obtenue  de  fagon  approchee 
p«Lr  resolution  de  I’equation  de  Zakai  discretisee  (cf.  (Le  Gland]),  nous  aidera  a 
interpreter  les  resultats  obtenus. 

Ce  rapport  est  organise  comme  il  suit; 

Dans  le  paragraphe  2,  nous  presenterons  les  deux  filtres  de  Kalman-Bucy 
(7)  et  (8)  utilises  et  nous  formulerons  des  remarques  preliminaires.  Dans  le 
paragraphe  3  nous  nous  interesserons  a  la  detection  des  passages  du  signal  par 
0.  En  analogic  avec  le  cas  continu  deux  tests  seront  proposes,  I’un  base  sur  les 
accroissements  des  observations  (cf.  {FJPj)  et  I’autre  sur  la  sortie  des  filtres  de 
Kalman-Bucy  (cf.  (Roubaudj).  Sachant  que  sur  un  intervalle  de  temps  [a, 6]  le 
signal  ne  passe  pas  par  0,  avec  une  probabilite  donnee,  nous  considerons  un 
autre  type  de  tests  pour  decider  du  signc  de  {xfc}fc  sur  [a,  6).  Sous  I’hypothese 
(II Dl),  des  etudes  ont  ete  faites  par  [FJSZj.  Ces  auteurs  ont  mis  en  oeuvre  deux 
tests  differents:  I’un  base  sur  la  variation  quadratique,  I’autre  du  type  rapport  de 
vraisemblance  base  sur  les  sorties  des  filtres  de  Kalman-Bucy.  Nous  resumerons 
brievement  ces  procedure.^  dans  le  paragraphe  4.  Sous  I’hypothese  (HD2),  nous 
presenterons  dans  le  paragraphe  5  un  test  de  rapport  de  vraisemblance  base  sur 
les  accroissements  des  observations  et  adapterons  le  test  sur  les  sorties  des  filtres 
de  Kalman-Bucy  du  paragraphe  precedent  au  cas  traite. 

La  mise  en  oeuvre  des  tests  de  detection  et  de  decision  ainsi  decrits  necessite 
la  determination  de  formules  explicites  pour  les  homes.  Bien  qu’une  demons¬ 
tration  rigoureuse  n’ait  pu  etre  faite.  les  formules  proposees  sont  justifiees  de 
maniere  heuristique.  Dans  le  cas  oii  /?_  ^  B+,  nous  emettrons  I’hypothese 
supplementaire: 

(//4)  <  0  et  <  0 . 

Cette  hypothese  n’est  pas  trop  restrictive  pour  notre  propos.  En  effet,  suivant 
le  signe  des  coefficients  de  derive,  nous  avons  deux  comportements  distincts  du 
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siganl.  Si  B+  et  5_  sont  negatifs,  il  devient  stationnaire  des  que  I’intervalle 
de  monotonie  est  assez  grand.  Sinon,  il  pent  “fuir”  vers  I’infini.  En  fait,  nous 
etudions  la  situation  la  plus  delicate  au  sens  ou,  la  probabilite  pour  que  le 
processus  passe  par  zero  est  importante.  Cependant,  si  les  valeurs  1J3_| 
et  \B+  \  sont  grandes,  {x*}  aura  tendance  a  changer  rapidement  de  signe  et  les 
tests  n’auront  probablement  pas  suffisamment  de  temps  pour  decider. 

Le  paragraphe  6  est  consacre  a  I’application  numerique  de  ces  methodes.  Des 
criteres  de  comparaison  entre  les  differents  tests  seront  proposes  et  les  resultats 
des  applications  a  divers  exemples  seront  presentes.  D’autre  part,  la  solution 
de  I’equation  de  Zakai  sera  utilisee  pour  la  justification  du  comportement  des 
filtres  approches. 

Pour  conclure,  on  presente  dans  le  paragraphe  7  quelques  reflexions  sur  le 
comportement  de  ces  filtres  et  on  discute  la  performance  des  differents  tests 
etudies. 

Notation  1.1  Elant  donne  un  processus  {Gji.}  on  ecrira 

oil  q  ^  1R+,  pour  signifier  qu  hi  ciislc  Cj  ,C2  ,C3  >  0  fels  que 

<  ci  ox]){—C2fcs}  +  C3£'^’  ,VA:  G  A''. 

Notation  1.2  On  nolcra  c  on  C  des  conslantes  independantes  de  e,  sans  se 
soucter  de  letir  valeur.  Ces  constantes  pourront  done  prendre  des  valeurs  differ- 
cnlcs  d’une  ligne  d  Vautre. 
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Deux  filtres  de  Kalman— Bucy  en  parallels 


Soit  ^0 ,  la  tribu  des  observations  jusqu’a  I’instant  k,  =  <T{yo,  t/i,  •  •  • ,  yk}- 
La  solution  du  probleme  de  filtrage  associe  au  systeme  (1)  etant  donnee  par 
la  loi  conditionnelle  de  Xk  sachant  nous  commentons  brievement  son  com- 
portement  (cf.  [FJSZ]).  Sous  I’hypothese  de  “detectabilite”  (HD),  la  variance 
de  cette  loi  est  “petite”,  sa  densite  se  concentrant  autour  de  deux  maxima 
locaux,  Pun  negatif  et  I’autre  positif.  Ce  phenomene  se  justifie  du  fait  qu’on 
a  soit  Xk  —  ykIH-  =  \/£fH_Vk  soit  Xk  —  yklH+  =  y/ejH^Vk,  done  la  den¬ 
site  conditionnelle  est  “petite”  excepte  sur  deux  intervalles  d’amplitude  d’ordre 
0{y/£)^  autour  de  yklH^  et  de  yk/H+.  En  emettant  I’hypothese  de  detectabilite 
(HD)  nous  nous  pla^ons  dans  la  situation  ou  lorsque  {ijt}  prend  ses  valeurs 
loin  de  0  durant  un  certain  intervalle  de  temps,  un  des  pics  tend  a  disparaitre, 
la  loi  conditionnelle  s’approchant  d’une  loi  gaussienne.  II  apparait  alors  legitime 
d’approcher  le  filtre  optimal  par  un  filtre  de  Kalman-Bucy.  Intuitivement  on 
s’attend  a  ce  que,  sous  I’hypothese  {HD2),  le  temps  necessaire  pour  faire  “dis¬ 
paraitre”  I’un  des  deux  pics  soit  plus  long  que  sous  I’hypothese  (HDl).  Cette 
idee  sera  illustree  dans  le  paragraphe  6  (voir  les  figures  3  a  6  et  12  a  14). 

Soient  {xf,Q^)  et  (Xk^Qk)  les  filtres  de  Kalman-Bucy  associes  respective- 
ment  aux  systemes  lineaircs 


^k+i  =  ( I  +  e)  Xk  +  ^/e(T+  xik 

XJk  =  //+.Tit  +  y/s  Vk  ,2/0  =  0 

(5) 

•Tfc+ 1  =  ( 1  -b  e)  .rt  -b  “Jt 

(6) 

tjk  =  fl-Xk  -r  y/evk,  yo  =  0, 

avec  conditions  initiates  gaussienne.s 


2^0  =  io  =  ^i^o) 


gj  =  £ 


1  -b  Hlal 


1  Qq  —  - 


1  +  ' 


11s  sont  donnes  par  les  equations  suivantes: 

if,,  =  (l+jfl,|if +  t//,(3nsw.-«+(l  +  £S+)ii) 

o*  -  +  <^1^ 

I  '  (1  +  +  f/l(l  +  c 


o 


et 


'h+i  =  {^  +  eB.)x]; +  ^H_Qi^{yk+i- H.{l+eB_)x-;;) 

o-  (1  +  e  B-)^Q;  +  (Tie 

■  (1  +  Hiai)e  +  Iii{\  +  e  B  ‘ 

Remarque  2.1  Le  fait  d’avoir  considere  une  condition  initiale  gaussienne  pour 
les  systemes  (5)  et  (6)  n’est  pas  restnctif  pour  le  probleme  de  depart.  En  ef- 
fet,  puisque  les  filtres  etudies  se  caracterisent  par  le  fait  d’avoir  une  “memoire 
courte”,  I’influence  de  la  condition  initiale  tend  a  disparaitre  rapidement. 

On  note  les  gains  a  I’instant  k: 

A?  =  -i.*Qt  01  =  '-ii-Q;  ■ 

e  e 

Remarque  2.2  Dans  une  situation  stationnaire,  les  expressions  des  variances 
des  lois  conditionnellos,  Q'^  et  Q",  sont  facilement  calculees.  Elies  sont  d’ordre 
0(e)  et  done  les  gains  stationnaires  K+  et  sont  d’ordre  0(1)  (cf.  [Milheiro]). 

On  introduit  les  processus  dits  “d ’innovation’’ 

=  Uh+i  ~  //+(1  +sB+)xt 

et 

‘'k+\  =  .'At I  -  //-(I  +  £  B.)x;  . 

Les  processus  ainsl  definis  sent  considcres  approximativement  comme  des 
"bruils  blancs"  au  sens  oil  leurs  variances  sont  asyinptotiquement  d’ordre  0(e), 
landis  quo  les  correlat.jiis  et  sont  d’ordre  O(e^). 

On  notcra  edj  ct  ei)'^  leurs  variances  respcctives  el  ex)'^  et  ei)~  les  variances 
stationnaires  c;orrcspondante.s: 

£!?+  =  e(\+(Tllll)-kHli\-keB^)^Q-^, 

ex)-  =  ;{\  -\- ailli)  +  lH{\ +e  B.fQ- ■ 

Dans  la  suite,  on  eludiera  la  procedure  de  test  permettant  de  separer  ies 
intervalles  de  luonolonie  de  la  fond  ion  It  et  de  decider  du  signe  de  Xk-  Sur 
cliacun  de  ces  intervalles,  il  sera  alors  po.ssible  d’approcher  le  filtre  optimal  par 
le  filtre  oe  Kalman- Bucy  correspondant. 
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Tests  de  detection  des  passages  de  {a:jb}  a 
zero 

Nous  devons  premierement  detenniner,  sur  I’intervalle  de  temps  [0,  T],  les  in- 
tervalles  sur  lesquels  {xjt}  cst  different  de  zero,  avec  une  probabilite  donnee, 
proche  de  1.  Nous  presentons  deux  tests  de  detection,  le  premier  s’appliquant 
aux  observations  j/^,  le  second  a  la  sortie  de  I’un  des  filtres  de  Kalman-Bucy. 

Soit  I’intervalle  de  temps  [a, 6],  avec  0  <  a  <  b  <  T.  On  pose  m  =  [(6  —  a)/e] 
et  io  =  [n/e]-  On  considere  les  deux  evenements  suivants: 

A.  =  {xk  <  0  \  k  =  io,io  +  l,- ■  ■  ,io  +  m], 

A+  =  {xk  >  0;  k  =  iQ,iQ  +  i,-.  -Jo  +  m}. 

On  remarque  que 

(.4_  U  /l+)‘^  =  {  XkXk+i  <  0  pour  un  certain  k  ,?o  <  <  lo  + 

3.1  Test  base  sur  les  observations  ijk 

Etant  donne  que  ‘7i(a;)  =  0  ssi  a-  =  0“,  h{xk)  “petit”  implique  par  continuite 
do  I’application  /i,  Xk  proche  de  zero.  Cependant  h(xk)  n’est  pas  observe  mais, 
pour  £  “petit”,  il  est  approche  par  ?/<... 

Soit  Cofcj  >  0  une  constante  a  determiner. 

On  definit  I’evcnement  test  suivant: 

C  =  {li/rl  >  Co6,}- 

On  montre  de  fa?on  siniilaire  a  la  version  en  temps  continu  (cf.  [FJP][proposition 
3.1j)  la  proposition  suivante  (cf.  (f'.JSZjfproposition  2.1]): 

Proposition  3.1  Soient  c^os  >  0  r/  cq  >  0  doiines.  II  existe  0  >  0  tel  que 
Vs  G  (0.  £oj 

Determination  de  la  constante  Cobs- 

Nous  utilisons,  pour  la  determination  de  la  constante  Cota,  un  raisonnement  sim- 
ilaire  a  celui  propose  dans  [FJSZj,  sauf  que  dans  notre  cas  le  coefficient  de  derive 
n’est  plus  suppose  constamt  mais  constant  par  morceaux.  Intuitivement  le  choix 
de  Cobs  (ioit  etre  un  compromis  entre  “conserver  des  intervalles  de  monotonie 
suffisarnment  grands”  et  “reduiie  la  probabilite  d’erreur  du  test”. 
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On  cherche  a  ;e  que  la  probabilite 

P{{xi:X^+l  <  0}  I  {\yk\  >  Col,,,  |2/A:+i|  >  Co6,}) 

soit  faible. 

D’apres  Thypothese  (//3)  on  obtient: 

P{{yk  <  -C06.))  <  P{{vk  <  -Coislyfe))  <  ^  . 

Nous  nous  permetlons  done  de  considerer  uniquement  le  conditionnement  par 
rapport  a  I’evenement  {?/;.  >  Cobs,yk+\  ^  De  plus,  on  a: 

<  0}  O  {uk  ^  (^obs'^yk+l  ^  =  ^Cob,  ^ 

avec  les  definitions  suivantes  pour  ies  evenements  D],  D^: 

D\  ~  {a:*:  <  0,.rfc+i  >  0}  n  {j/fc  >  c,yfc+i  >  c} , 

D]  =  {xk  >  O.xk+i  <  0}  n  {ijk  >  c,yk+i  >  c]  , 
pour  tout  c  >  0 . 

On  considere  separement  Ics  evenements  D\^  et  Dl^. 

•  Sur  Dl^ ,  on  a  Ics  recuncnccs  suivantes: 

Xk+l  =  (1  +  s  D.  )Xk  +  <T_  \/£  Hfc, 

Uk  =  ILxk  -r  \/sVk,  ijk+i  =  H+Xk+i  +  \/?Ui.+  i  . 

On  deduit  Fegalite: 

-^{/Lyk+i  -  //+{1  +  P-)yk)  =  ll+ll-cr^uk  -V  !I-Vk-i  \  -  ^+(1  +£  B.)vk . 

Or.  d’apies  [ll'-V)  on  a  h{xk)  >  0  et,  pour  f  suffisamment  petit,  on  obtient  la 
majoration 

^\Il.~H4l+eD_)\cy<\Zk\, 

oil 

Zk  =  II+ILa.Uk  +  H-Vk-ki  ~  //+(]  ■r£D^)vk 
suit  line  loi  gaussienne  yV(0  .  //+//^cr^  +  //^  +  //|(1  + 

Soit  (Z°}  le  processus  normalise  as.sorie  a  {Zk}-  Pour  un  risque  a  donne  C), 
on  determine  A  >  0  tel  que 

P(\zl\  >  A)  <  a, 

^Generalement  on  prend  a  =  0.05 
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par  inversion  de  la  fonction  de  repartition  de  la  loi  normale  centree  reduite.  On 
obtient  ainsi  I’expression  de  la  constante  ci: 


Cl  =  Xy/e- 


iHlHlal  +  HI +Hl{l+e  B_y 


•SurC.V  un  raisonnement  similaire  au  precedent  nous  amene  a  considerer 
une  constante  donnee  par: 

'HlHial  +  Hi  +  Hi{l  +  eB+y 

\H^  -  H.{1  +  e  B^)\  • 

Dans  les  applications  numeriques.  pour  une  probabilite  d’erreur  a  donnee, 
nous  prendrons  Cohs  —  max(ci,C2). 

Remarque  3.2  Dans  les  calculs  de  la  constante  la  prise  en  compte  de  la 
derive  introduit  uniquciiicnt  dcs  ternies  d’ordre  que  nous  negligerons 

devant  les  tennes  d'ordre  0(y/e}. 


3.2  Test  base  sur  la  sortie  d’un  filtre  de  Kalman-Bucy 

On  se  propose  do  decrire  un  test  capable  de  detecter  les  changements  de  signe 
ovcntueis  du  processus  {a'l,.}  a  I’aide  de  la  sortie  obtenue  par  un  des  deux  filtres 
de  Kalman-Bucy,  par  oxempie  {.rj}-  La  construction  d’un  tel  test  se  justifie  par 
)e  fait  que  ''h{x]  0  ssi  x  —  0"  et  par  la  proposition  suivantc: 

Proposition  3.3  Soil  0  <  a  <  b.  Pour  tout  0  >  0,  il  existe  j3  >  0,  eo  >  0  tels 
que.  pour  z  £  (O.so],  on  dcfinit  io  =  (a/ej  ci  m  =  [(6—  o-)je],  on  a: 

/’({  ?nax  >  ^})  ^  . 

fc=io,  •••■o+m 

En  vue  de  la  demonstration  de  cette  proposition  on  etablit  d’abord  le  lemme 
suivant: 

Lemme  3.4  Pour  tout  r  >  0.  il  exist (  So  >  0  et  c ,  C  >  0  tels  que,  pour  tout 
e  e  (O.coj,  m  =  [rjz], 

E{exi>  max  cxl)<C. 

A=:0.  ,m 


Preuve 

La  preuve  de  ce  lernme  sera  presentee  en  3  etapes  (cf.  [Roubaud]  pour  la 
preuve  analogue  en  temps  continu). 
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Etape  1: 

On  demontre  le  resultat  suivant: 

Lemme  3.5  On  definit  le  processus  {zk}  par 

Zk+i  =  (1  +  £B)zk  +  \/e  (xuk  ,  (9) 

on  zo  est  une  variable  alcatoire  independante  du  processus  {ujt}  et  il  existe  des 
constantes  Cq  ,  Cq  >  0  telles  que 

E(expco2o)  <  Q-  (10) 

Alois,  pour  tout  r  >  0,  il  existe  £0  e/  c,  C  >  0  tels  que,  pour  tout  e  6  (0,£o], 
m  =  [rje:], 

E(exp  max  c  2^)  <  C  . 

k=0,—,m 


Preuve 

Nous  reccrivons  I’expression  de  2^+1  cn  mettant  en  evidence  une  partie 
depcndante  de  zq  et  une  partie  martingale: 

2,+,  =  (1  +  +  (!  +  eDfMk  , 


avcc 


Mk  ~  os/e  ^(1,  +  . 

Puisfiuc  Zq  et  {Mk}  son  independanls  nous  avons: 

zUr  -  (1  +  +  (I  +  eBf’^Ml  <  +  M^) 

D’apres  tlO).  il  est  suffisant  de  montrer  qu’il  existe  eq,  c.  C  >  0  tels  que 
Ve  <  Eo  •  E(exj)  max  r  M?)  <  C  . 

AsO.  -.m - 1 

Or  est  un  processus  gaussieii  centre  de  variance 

2  (l+£^?)-(l  +£/?)-•'• 


done 

si  B  =  0,  E[A'q]<ksa^ 


B(2  +  eB) 
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si  5  >  0,  E[M^?]  <  a' 


si  B  <  0,  E[Mfc2]  <  <r' 


,(l  +  f5r 


B{2  +  eB) 

-(1  +  eBf  +  exp{2keB  +  ke‘^B‘) 


-B{2  +  eB) 

Dans  I’intervalle  fini  [0,t]  il  existe  done  Cq  et  une  constante  C  >  0  tels  que 

Vs  <  £o ,  max  E(A/?]  <  C . 

fc=0,m 

d’ou  on  deduit  cq  et  des  constantes  c,  C  >  0  tels  que 

Vs  <  So  -  E(expcM^_j)  <  C . 

Vu  que  {exp  c  Ml)  est  une  sous-martingale  on  pent  appliquer  la  propriete 
E(  max  exp  cA/j?)  <  E(exp2cM^_j) . 

fc=0,-".m-l 

Alois,  il  existe  So  et  des  constantes  c,  C  >  0  tels  que 


Vs  <  So ,  E(^_max  ^  expcMl)  <  C , 


d’ou  le  lemme. 


Etape  2: 

On  introduit  les  notations: 

B  =  sup(|B+|,i5_|) 

<r  =  sup(|<T4.|,|a_i) 

On  demontre  que  ics  moments  d’ordre  2p,  p  >  1.  du  processus  {xjt}  peuvent 
etre  majores  par  les  moments  d'ordre  2p  du  processus  [zk]  defini  dans  le  lemme 
precedent. 

En  effet.  d'apres  requatioii  d’etat,  on  a 

j=o 

Le  fait  que  {u^}  est  un  bruit  blanc  et  que  ujt  est  independant  de  x*  entraine 
I’egalite 

E|4:,l=  E  + 

}=Q-,]  pair 


3 

% 


■€- 
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On  obtient  la  majoration 


E[4;.i  <  i; .  ■ 


j=o;j  pair 


Or,  de  la  definition  du  processus  {-zjt},  vient  que 


j=Oij  pair 
On  prend  zq  =  l^ol-  Alors 


E(xf|<E(r?j,t  =  o,i,.-.,|rAi. 
Le  lemme  3.5  entraine  I’existence  de  £o ,  c,  C  >  0  tels  que 


Ve<£oi  E(expca:4.)  <  C  . 

h=0,—,m 


Etape  3: 


Nous  utilisons  la  majoration  (11)  pour  obtenir  le  resultat  pretendu,  le  lemme 
3.4. 

D'apres  la  formule  de  Taylor,  on  a: 

exp  c Xj^  1  =  exp c xj  •-  2c Xj  exp  cXj(xj+i  -  Xj) 

-i-2c  cxpcOj(}  +  2c0j)(xj+,  -XjY  . 

On  obtient  rcxpression 

k 

e.xp  cifc+i  =  e.xp CX5  +  2c  ^  cxpic xl_j){xk+i-j  -  Xk-j) 

j=0 

-r2cY,  T^k-jCxp{c0l_j){l  +2c0l_^)ixk+i...  -Xk-j)^. 

j=0 

Puisquc  <  |.rj|;|  +  ixjt.^i|,  on  a  la  majoration 


El^.J»ax_^expc.r?.+  ,j 


<  E[exp  c  x^)  +  2  cc  ^  <'-’^P(c4-j)) 

3=0 
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+2cE 


m-1 


/  k  \ 

^  E exp(cCj)  (^ixk-:)uk-j 

’  '  \i=o 


(12) 


—  A 

+4  c  e  E  E  [exp  (2  c  (xL,  +  4+i-,))  (l  +  4  c  •  +  xL,))  {e  8^%  +  a\l_.)\  . 


Le  processus 


j=o 


est  une  J^<:-martingale,  ou  =  cr(xo,  Uq,  •  •  • ,  Uk). 

On  utilise  une  propriete  cles  martingales: 

Or  Uk  et  Xk  sont  independants  et  Xk  cst  ^k-i  ~  mesurable  done  on  a 

EM/L.j  <  "’^^^^o"'.^^_,E(exp(4cxL,.;)]'^'E[xL, <  c,, 

d'apres  (11). 

II  cxiste  alors  £0  >  0  et  c.  C  >  0  tels  que,  pour  tout  e  G  (0,£o], 


E 


/  r  ^ 

V;=o 


<  C. 


On  utilise  I'incgalite  de  Cauchy-Sclnvarz,  le  resultat  (1 1)  et  I’hypothese  (//2) 
dans  I’inegalite  (12)  pour  deduire  le  Icmme. 


Preuve  (de  la  proposition  3. .3) 

Nous  utilisons  la  notion  de  sous-difFcreiiii-,.!  d"  h. 

Puisque  la  fonction  h  est  convexe.  il  existep^  appartenant  au  sous-differentiel 
de  h  a  un  point  Ok  tcl  que 

/‘ixk+] )  =  h(xk)  +  PA-(.rjt+i  —  Xk) 
et 

iA)<sup(j/4i, !//_))  ^;/. 

On  definit 

=  liixk)  -  . 
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D’apr^  I’equation  du  filtre  et  I’equation  d’etat,  on  obtient  I’expression  re¬ 


cursive: 


Zk+i  =  SZk  +  eb{xk)pkil  -  H+K+)  -f  v/e(7(xit)pjt(l  -  H+K+)uk 

-f  \/e  H+K+Vk+i , 
ou 

<5  =  (1 -f  e5+Hl  - //+A'+)  • 

On  remarque  que  0  <  6  <  1.  On  utilise  cette  recurrence  jusqu’au  terme  Z,,,.  On 
etablit  ainsi  la  inajoration: 


k-.o  T°.o  (^3) 

-f  C  \fe,  ^  +  C  y/e  ^  (S-’lu/t+i.jj . 

;=0  j=0 

On  traite  les  quatre  termes  dans  cette  expression  separemmeni. 

D’apres  I’inegalitc  de  Bienayme-Tchebichev  et  le  fait  quo  B[!Z,j,lj  <  C,  on 
obtient 


D'autre  part,  en  utilisant  Ic  leinnie  3.4,  on  obtient  la  inajoration 

Lc  Icminc  2.3.  dans  [Roubaudj  cutrainc  la  inajoration 

P[  ma.x  iu;|  >  j-  —  <  - e'  ^  , 

\^fc=io.- -.loJ-m-I  4C^£  /  e 


(k-io  0  \ 

pour  £  ‘■sufRsaniment  petit". 

De  meine  pour  le  4""’  terme  a  droite  dans  I’expression  (13): 
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De  I’expression  (13),  en  utilisant  Ics  4  majorations  ci-dessus,  on  obtient  le  re- 
sultat. 


Soil  cpj;  une  constante  strictement  positive  a  determiner. 

On  detinit  I’evenement  test  suivant: 

^  -•  {\^t\  ^  <^Fi<  ;  ~  +  1,-  ••,io  +  m}  . 

Comme  consequence  de  la  proposition  3.3,  on  obtient  le  theoreme  suivant: 

Theoreme  3.6  /'  exi^te  /?  >  0,  co  >  0  tels  que,  pour  e  6  (0,  fo),  oji  a: 

Preuve  La  preuve  de  cc  thioreirte  suit  pas  a  pas  celle  du  resultat  equivalent 
en  temps  continii  (cf.  fRoubaud]  [Theoreme  3.2.]). 


Determination  de  la  constante  cpK' 

On  cherche  une  constante  telle  que  la  probabilite 

P  {{Xk-Vk^l  <  0}  !  {ji-JI  >  Cf-K  ,  >  C/r;,  }) 

•SO it  faible. 

On  considore  pk  coinme  une  approximation  de  h{xk)-  Dans  I’etude  prece- 
dentc,  pour  une  probabilite  d'erreur  a  donnce,  on  a  determine  une  constante  Coba 
permeltant  de  tester  au  vu  de  pk  I’occurence  d’un  changement  du  signe  de  Xk-  Or, 
le  cliangement  de  signe  de  x*  est  directement  lie  au  fait  que  h{xk)  devient  “pe¬ 
tit".  Si  //+xj  etait  une  meilleurc  approximation  de  h(xk)  Que  Pk,  tout  au  moins 
aux  alentours  de  zero,  la  probabilite  d'erreur  du  lest  “xj^  >  CobaPK''  serait  mi- 
noree  par  a.  Cependant,  nous  ne  .savons  pas  si  cette  condition  est  verinee  mais, 
sur  revenement  .4+.  le  filtrc  x;];  est  une  meilleurc  approximation  de  Xm  que 
t/,rt///+.  En  tout  cas,  la  proposition  3.3  nous  assure  que  I’ccart  |n(xjt)  —  H+x'l\ 
est  petit  avec  une  probabilite  proche  de  1.  II  semble  alors  raisonnable  de  penser 
que  le  test  “.rj  >  Co^a! permeltra  de  detecter  convenablement  les  passages 
de  Xk  par  0.  .N'ous  proposons  done  pour  constante 

Cf'l\  —  C(,bal  [//^.l  . 
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Le  meme  type  de  raisonement  est  valable  quand  le  test  est  applique  sur  le 
filt'e  xj^  au  lieu  du  filtre  x^ .  Nous  prenons  alors 

CPK  =  Cobs/  l^-l  • 

Sans  perte  de  generalite  on  considerera  par  la  suite  un  seul  intervalle  [a,  6] 
representant  I’intervalle  de  temps  durant  lequel  il  n’y  a  pas  de  passage  a  0,  avec 
une  probabilite  d’erreur  donn^,  io  et  m  etant  definis  comme  en  debut  de  cette 
si  :tion. 
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4  Sous  I’hypothese  (HDl):  decider  du  signe 
de  Xk 

Dans  cette  section,  nous  pr&entons  deux  tests  permettant  de  choisir  entre 
les  evenements  A+  et  /!_,  dans  le  cas  oil  le  probleme  de  filtrage  lineaire  par 
morceaux  verifie  I’hypothfee  (HDl).  Ces  tests  ont  ete  etudies  dans  [FJSZ].  Le 
probleme  analogue  en  temps  continu  a  ete  traite  dans  [FJP].  Le  premier  test 
considere,  dit  de  la  variation  quadratique,  n’est  applicable  que  sous  I’hypothese 
{HD\).  Le  second  test,  du  tj'pe  rapport  de  vraisemblance  sur  les  sorties  du 
filtre  de  Kalman-Bucy,  pourra  etre  adapte  au  cas  oil  I’hypothese  {HD2)  est 
verifiee  (cf.  paragraphe  5.2).  Le  premier  test  sera  done  introduit  de  fagon  breve. 
Le  deuxieme  test  sera  decrit  de  fagon  plus  detaillee,  puisque  la  meme  id^  sera 
utilise  par  la  suite. 

Dans  [FJSZ],  les  auteurs  ont  mis  en  ceuvre  ces  tests  sous  deux  formes:  test  a 
taille  d’echantillon  fixee  et  test  sequentiel.  D’apr^  leur  etude,  les  tests  sequen- 
tiels  semblent  plus  interessants  que  ceux  a  taille  d’echantillon  fixee,  vis  a  vis 
du  critere  du  “temps  moyen  pour  prendre  une  decision”.  Nous  ne  presenterons 
done  que  les  tests  sequentiels. 

4.1  Test  de  la  variation  quadratique 

On  considere  un  test  d'hypotheses  base  sur  la  variation  quadratique  de  la  suite 
d’observations  {?/*}  pour  decider  entre  les  deux  alternatives  “,4+”  et 
Cette  decision  est  prise  sur  A'  observations  (0  <  N  <  m)  dans  un  intervalle  de 
monotonie  obtenu  par  un  des  tests  introduits  dans  la  section  3.  On  considere  A’ 
comme  etant  un  temps  d’anet. 

Sans  rcstreindre  la  generalitc  du  probleme.  on  suppose  que  B+  =  D-  =  B . 

Remarque  4.1  Le  fait  de  considerer  B^  et  /?_  distincts  a  pour  unique  conse¬ 
quence  i'introduction  de  termes  negligeables  dans  les  expressions,  la  procedure 
restant  essentiellcment  la  nierne. 

Notons: 

Ai  =  t/k*;  -  -r  Bs)  7/k , 

Ti  =  I  4  (1 

Ti  =  IIl^l  +  I  A  (I  +  Bt)\ 

A  lors 

•  si  Xu  >  0  et  xu+i  <  0.  la  variable  aleatoire  A*  suit  une  loi  gaussienne  de 
nioyenne  nulle  et  variance 
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•  si  Xfc  <  0  et  xjk+i  <  0,  la  variable  aleatoire  At  suit  une  loi  gaussienne  de 
moyenne  nulle  et  variance  Tie. 

L’hypothese  de  detectabilite  (HDl)  implique 

y7  u.  T2 

Le  probleme  de  decider  entre  les  alternatives  “.4+”  et  conduit  a  un  test 

d’hypotheses  sur  ces  variances. 

On  definit  le  rapport  logaritliniiquc 

Af(Ai,0,Ti£)’ 

oil  -V(x, represente  la  densite  de  la  loi  gaussienne  de  parametres  fi  et  7^ 
et  la  serie 

*0+n 

s„  =  5:^*, 


oil  id  et  m  sont  definis  dans  la  section  .3. 
On  obtient  les  e.xpressions  suivantes: 


o°T+  '^2  VTi  Ti 


A? 


ei 


-  (n  4-  l)log  j  g  ^  . 

•Soient  /;  >  0  et  /j  >  0  fi.xcs.  On  intioduit  le  temps  d’arrct  A'* : 

A”  =  inf{n  >  0  :  S„  <  — /,  ou  5„  >  /j}  A  m . 

On  decrit  le  test  comine  suit:  si  5.v  >  U  on  acceptc  “.4+"  et  si  5,v.  <  --!i  on 
accepte  "/L";  sinon  on  dira  quo  le  te.st  ne  pcrmct  pas  de  decider. 

11  reste  a  proposer  une  mcthode  de  determination  des  homes  /j  cl  I2.  On 
approclie  le  processus  di.scret  par  un  certain  processus  dc  diffusion 

dont  on  ecrit  les  equations  sous  et  sous  On  obtient  ainsi.  par  la 

.tormule  de  Dynkin.  des  equations  differenticllcs  ordinaires  permettant  de  cal- 
•^uler  les  temps  moyens  pour  prendre  une  decision  EfT^)  et  E(7i)  et  pour  les 
probabilites  d’erreur  p+  et  p_  (cf.  (Fleming- Rislielj). 


1 
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1  « 


On  introduit  les  notations 


•  esi  "la  probabiiile  dc  refuser  a  tort". 

•  P-  esl  "!a  probabiiile  de  refuser  .•t_  a  tort". 

•  7’*  et  7  *  sonl  les  temps  d'arret 

r*  =  inf{/  :  (,'f  >  5  f;  ou  Q  <  — c  /j } .  dans  ie  cas  positif. 


el 


7’*  =  inf{t  :  Q  >  s  I2  ou  (.  <  — s  ,  dans  Ic  cas  ncgatif  . 


Pour  des  probabiiilcs  d’erreur  et  p_  donn«5es  (■*).  les  expressions  (14)  et 
(15)  pcrinettent  dc  calculcr  les  borncs  /j  et  /j. 

■^Gcncralemcnt  on  prend  p+  =  p_  =  0.05  . 
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4.2  Test  du  rapport  de  vraisemblance  sur  les  sorties  des 
filtres  de  Kalman-Bucy 

On  decrit  un  autre  test  du  type  rapport  de  vraisemblance  pour  decider  entre  les 
deux  alternatives  et  .  L’idee  de  ce  test  est  la  suivante: 

Soit  un  entier,  to  <  Intuitivement,  on  s’attend  a  ce  que 

hk  =  E{h{Xk)\yQ)  soit  tres  proche  de  H+xf  sur  A+  et  de  H.x'^  sur  A_,  pour 
h  k  <  t'o  +  m.  On  a 

sur  A+  :  rjk  =  H+{\  +  eB+)x'^_i  +  i/^ , 

sur  v4_  :  yk  =  //_(1  +  +  i/J, 

{^k  }  {^k  }  peuvent  etre  consideres  comme  des  bruits  blancs  gaussiens  dits 

“processus  d’innovation”  (cf.  section  2). 

On  pose 

Zk  =  //+(!+  eB+)xt  -  H^l  + 

et  on  prend  pour  statistique  de  test  le  rapport  logarithmique  suivant: 

^  L^-ail  “  kzzii 

avec  t’l  <  n  <  m. 

La  statistique  />„  se  reecrit  sous  la  forme  suivante: 

1  i|+n  1  »i+n 

Sur  A^  -.L  =  E  2*  +  7  E  ■ 

k=il  ^  k=i, 

Sur  A_  :  L„  =  ^k  +  -  ^kt'k+i  ■ 


On  note  le  terme  de  precision  du  test: 

1  il+n 

r..  =  -Y,zI 

^  k=ii 

Pour  Rn  “suffisamment  grand”,  le  signe  de  permet  de  choisir  entre  les  alter¬ 
natives  “>1_”  et 

Soient  /j  >  0  et  /j  >  0  des  homes  fixees.  On  introduit  le  temps  d’arret  N" : 
=  inf {n  >  t'l  :  Ln  >  h  on  Ln  <  -li}  /\  m . 

On  applique  la  regie  de  decision  qui  suit. 
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Regie  de  decision: 

•  Si  Ln>  >  I2,  on  decide 

•  Si  <  — /i,  on  decide  “y4_”. 

•  Sinon,  on  ne  pent  pas  decider. 

II  nons  reste  a  calculer  les  homes  It  et  I2  pour  une  probabilite  d’erreur  donnee 
a  et  a  estimer  les  temps  moyens  pour  atteindre  une  decision.  Le  raisonnement 
presente  par  la  suite  constitue  une  justification  purement  heuristique,  sans  la 
pretention  d’etre  une  preuve  rigoureuse. 

Si  on  se  place  dans  le  cadre  de  I’hypothese  d’apres  les  equations  des 

filtres  de  Kalman-Bucy  (7)  et  (8),  on  a: 

+{H+h\il  +  eR+)  -  +  eB-))4+t.  (16) 

Vu  que  les  filtres  en  question  sont  a  memoire  courte,  d’apres  cette  expres¬ 
sion,  I’hypothese  d’inclependance  sur  A'*'  de  et  entraine 

asymptotiquement  “I’independance”  de  Zk  et  a  partir  d’un  certain  instant 
it  >  io. 


Remarque  4.2  Pour  s  “assez  petit”,  on  a  la  majoration 

On  obtient  alors  que 

•  si  R+  =  i?_  =  B,  dans  une  situation  asymptotique, 


ElZ^;/l+j 


'■  I  -  (1  +  eB^il  - 


•  si  R+  <  0,  alors  ;  A+]  admet  une  solution  d’equilibre: 


E(xf  ;  .4+)  = 


Kl 


|5+((2  +  eB+) 
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On  utilise  I’inegalite 


E[Zi  *; ;  At]  <  ;  ■4+1  +  ^  Elxf  ;  A+l) 

pour  obtenir  I’estimation 

mir ; 

<  [(1  +  -  H.K.f  +  s/lH+\B+  -  B_|(l  -  +  e5+)(l  +  eB.)]nzl ;  /1+] 

+£3/2  (1  +  £B+){\  +  eB_)\B+  -  B_l(l  -  i/_  A'_)//+£[x+2  ;  A+] 

+£  [H+K+il  +  £B+)  -  //_  A'_(l  +  £B.)]H-^ 

+£2(1  +  £5+)2(i?+  -  5_)2//2i;[i+2  ;  A+]  . 

Remarque  4.3  Pour  £  “assez  petit”,  on  a  la  majoration 

(1  +  eB.fil  -  II.K.f  +  \/£i/+|R+  -  R_|(l  -  //_A'_)(1  +  £R+)(1  +  eB.) 

<  1  -  //_iC_(2  -  //_  A'_)  +  c  x/?  <  1 . 

Dans  une  situation  stationnaire,  on  a  I’estimation 

E[Z2;/4]  <  A+£  +  C£3/2, 

^  //./C(2-//_A:_)  ■ 

Do  plus,  pour  un  certain  n^,  on  pent  coiisidcrer  Zk  comme  une  combinaison 
lineaire  de  processus  =  Zkt'k+i  approximativement 

un  processus  melangeant.  En  consequence,  on  approche  £A„  par  le  processus  de 
diffusion  {Ct;  t  =  en}  solution  de  I’equation  difTerentielle  stochastique 

(/(,'(  =  A+l'Zdt  +  y^£A+t?+  dW^  . 

Comme  pour  le  test  de  la  variation  quadratique,  on  obtient  la  probabilite  d’er- 
reur  p+  =  P({Ct;  =  -e/i}p4)' 

_  1  - 

~  f’Hh  _  ’ 

oil  est  le  temps  d’arret  defini  par 

Tl  =  inf{t  ;  Q>£l2onQ<  -e  /j} 
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et  rj^  =  l/t?+  . 

Si  on  se  place  dans  le  cadre  de  I’hypothese  par  un  raisonnement  iden- 

tique  an  precedent,  on  obtient  la  probabilite  d’erreur  p_  =■  P({Ct*  = 

Elle  est  donnee  par; 

_  1 

~  Qf)-h  _  ’ 

avec  =  1/t?". 

Pour  des  valeurs  de  et  p_  donnees  (“*),  on  determine  les  bornes  de  precision 
/i  et  I2  et  on  estime  les  temps  moyens  pour  prendre  une  decision: 


E(r;)  = 
E(r:)  = 


Remarque  4.4  Dans  [FJSZ],  les  auteurs  constatent  qu’en  pratique  le  temps 
moyen  pour  .  rendre  une  decision  est  plus  long  pour  le  test  de  rapport  de 
vraisemblance  <tue  pour  le  test  de  variation  quadratique. 


■^En  general  on  prend  p_  =  p+  =  0.05 . 
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5 


Sous  I’hypothese  {HD2):  decider  du  signe  de 

Xk 

Sous  I’hypothese  {HD2):  =  (//+(t+)^  et  jB_  ^  B+,  le  test  de  la 

variation  quadratique  n’est  pas  applicable.  En  remplacement  nous  considerons 
un  test  du  type  rapport  de  vraisemblance  sur  les  observations  {j/jt}.  Puis  nous 
adapterons,  sous  cette  hypothese,  le  test  du  rapport  de  vraisemblance  sur  les 
sorties  des  filtres  de  Kalman-Bucy,  presente  dans  le  paragraphe  precedent  sous 
rhypothese  (HDl). 

Les  calculs  qui  suivent  utilisent  le  fait  qu’une  situation  stationnaire  esl 
necessairement  atteinte,  raison  pour  laquelle  nous  avons  formule  I’hypothese 
suplementaire  {HA),  decrite  dans  I’introduction.  Nous  I’utiliserons  par  la  suite. 


5.1  Test  du  rapport  de  vraisemblance  sur  les  acroisse- 
ments  des  observations 

Ce  test  est  deduit  de  I’etude  du  cas  e  =  0. 

5.1.1  Le  cas  e  =  0 

On  considere  le  systeme  suivant: 

I  =  Xk  +  Af  b{.rk)  +  \/Ai(7{xk)  Uk, 

-  ^ 

[  Vk  =  h{xk). 

L’equation  d’etat  de  ce  systeme  generalise  celle  du  svsteme  discret  (4)  a  un  pas 
de  temps  At  . 

Le  probleme  est  de  determiner  le  signe  de  Xk  pour  k  ~  ?0)  +  b  •  •  • !  h  + 

Le  test  d’hypotheses  est  fait  sur  les  observations  {i/k,  +  I  <  k  <  io  +  m}. 

On  considere  la  filtration  (^l'^)k,  ou  io  <  k  < 

io  +  m. 

Sur  A+  :  ijk+i  =  (1  +  AtD+)yk  d  x/aIFu^  .  (18) 

oil  {uf;  k  >  Jo}  est  un  j^o'^tuit  blanc  gaussien  pour  une  probabilite  P'*'  telle 
quo 

P-^{D)  =  P{D),\/DcA+. 

Sur /!_  :  f/k+i  =  {I  +  AtB^)yk+ x/AtFuj; .  (19) 

oil  {wF;  k  >  I'o}  est  un  .F,^-bruit  blanc  gaussien  pour  une  probabilite  P~  telle 
que 

P-(D)  =  PiD),W  C  A_. 
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t’o+m 

Le  processus  k  >  io+l}  est  ^*"‘-niesurablc.  (i.e.  previsible)  et  ^  < 

k=’io 

+OC  P-p.s. 

On  utilise  une  version  discrete  du  theoreme  de  Girsanov  (cf.  [Milheiro])  et, 
vu  que  ^  B+,  on  prend  pour  statistique  le  loga'rithme  du  rapport  de  vraisem- 
blance: 

dP+  . 

Pn  =  ln 


dP-  Vo+"  ’ 

•o 


oul<n<m  —  1.  On  a  la  formule: 


I'o+r.  1  t'o+n 

(B+-B.)  ■£  Y.  ill 

^■=«0+l  K=>0  +  1 


Sur  I’evenement  A+,  la  statistique  se  reecrit  sous  la  forme: 

^  i  i_t  *  I.— 


k=io+'l 


A:=io+l 


Sur  I’evenement  /!„,  on  a; 


^n  =  --At^ - ^  y^  +  VAt - = -  ^  VkU^. 


r2 


^=10+1 


^=to+l 


(20) 


(21) 


(22) 


On  pose 


^  A:=io  +  l 


On  appclle  R„  le  terme  de  precision.  Pour  grand,  Ln/R„  est  proche  dc 
±  1/2  et  la  probabilite  d’erreur  est  faible.  Si  i/„  >  0  on  decide  et  si  Z-„  <  0 
on  decide  ‘M.", 

On  etablit  le  resultat  suivant: 

Theoreme  5.1  Soil  r  >  0.  Vn  tel  que  io  <n  <m  -  1,  on  a: 

P({Ln<0}n(J„>r}nA^})  <  e-^l\ 
P({i^„>0}niK„>r}nA_})  <  e-^l\ 


Preuve  La  demonstration  est  identique  pour  les  deux  inegalites.  Considerons 
par  exemple  la  premiere. 
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On  introduit  la  martingale  exponentielle 

(  _  D  P  io+k 


^  *  i=io+l  "  1=10+1 

avec  A  <  0  et  on  raisonne  comme  dans  le  cas  continu  (cf,  [Roubaud]). 


^‘1’ 


On  decrit  deux  procedures  de  mise  en  oeuvre  de  ce  test  sous  forme  sequen- 
tielle. 


Premiere  procedure: 

Pour  r  fixe,  on  introduit  le  temps  d’arret  N‘  : 

N'  =  mi{n>l  :  Rn>r}Am. 

•  Si  I/AT.  >  0  on  decide 

•  Si  Ln>  <  0  on  decide  “,4_ 

•  Sinon  on  ne  peut  pas  decider. 

Nous  devons  done  proposer  une  methode  de  determination  de  cette  constante 
r.  On  considere  le  probleme  sur  A+  et  sur  A-  et  on  choisit  r  =  max(r_,r+). 


Sur  Pevenernent  /!+>  Pose 


_  n  n  lo+n 

m;  =  v/At-  ----  ■  -■  Y.  ■ 

k=to  +  l 

La  statistique  L„  sc  reccrit  sous  la  forme: 

Ln  =  +  A/,f  . 

On  introduit  le  temps  d’arret  iV*  : 

=  inf{n  >1  :  /?„  >  r+}  A  m . 


On  a  que  /?;v*  est  approximativement  egal  a  . 

0.1  cherche  a  determiner  une  constante  r+  telle  que,  pour  une  probabilite 
1  —  a  donnee,  le  terme  de  precision  soil  significatif  dans  le  calcul  de  Dans 
d’autres  termes,  on  cherche  r+  lelie  que: 


>  ^r+}j>l+)  <  Q-. 


(23) 
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D’apres  (18),  on  reecrit  sous  la  forme  suivante: 


no— I 


=  (1  +  AtB+ryk+u.,u.  +  \/^r  (1  +  • 

•=o 

L’hypothfee  {H4)  nous  permet  de  considerer  t/jt  comme  une  combinaison  lineaire 
des  variables  aleatoires  gaussiennes  et  independantes  ujjlnp+i, '  ‘ = 
jjku'l-  T-e  processus  — 1  ^  ^  io+m}  est  alors  approximativement  u.n 

processus  melangeant.  De  plus,  la  valeur  stationnaire  de  E[y^ ;  A+]  est  approchee 
par 

.  p 

*  -B+(2  +  A(B+)' 

D’apres  le  theoreme  fonctionnel  de  limite  centrale  [Billingsley][th.20.l], 

^/Ki 

v^+  k=io+l 

oil  L  =  71  At  et  est  un  processus  de  Wiener  standard. 

On  deduit  de  ce  qui  precede  que  suit  approximativement  une  loi  normale 


jV  O.nAt 


-B+{2  +  AtB+)) 


Vu  que  /?,v*  ^  r+,  on  obtieut  la  relation  suivante  entre  r+  et  : 

(Z?+  -  B.f 


^  A'*  Ai 


-B42  +  AtB^)- 


Alors  I’incgalite  (23)  devient; 

P 

On  choisit 


/  I  \ 


'•+  =  -1  A\ 

avee  A  obtenu  par  la  table  de  la  loi  normale  centree  reduite. 
Sur  I’eveii'-iinent  .-L,  on  prend 

r_  =  4  A^ 

et  on  a  la  relation  suivante  entre  r_  et  A'*  : 

,  {B.  ~  B.f 

r_«vVAt  '  - 


^  - 
f 
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Exemple  5.2  Pour  —  —1,  B+  =  —0.25  610;  =  0.05  on  obtient  r  =  10.89. 
Pour  les  estimations  des  temps  moyens  on  obtient:  AtN^  =  9.66,  AtNZ  = 
38.52 . 

Deuxieme  procedure: 

II  semble  que  dans  le  test  precedent  le  temps  d’attente  pour  prendre  une  decision 
soit  “trop”  long.  En  efFet  les  bornes  ot  r_  sent  determinees  independamment 
I’une  de  I’autre.  Pour  reduire  ceite  duree,  nous  proposons  une  autre  procedure. 

Pour  des  constantes  /i  >  0,  4  >  0  a  determiner,  on  introduit  le  temps  d’arret 
N': 

N*  =  inf{n  >  io  t.q.  L„  >  ou  Ln  <  — /i}  A  m  . 

Regie  de  decision  : 

•  Si  La'*  >  /2  on  decide  “/!+”. 

•  Si  Z//V  <  —h  on  decide  M_”. 

•  Sinon,  on  ne  peut  pas  decider. 

On  utilise  a  nouveau  ^approximation  des  processus  discrets  par  des  diffusions 
pour  determiner  les  constantes  /i  et  l2- 

Sur  A+,  on  approche  par  t  =  Atn},  oil  est  le  processus  de 
diffusion  defmi  par  I’cquation  differentielle  stochastique 

iW. . 

On  introduit  le  temps  d’arret  T+  : 

r;  =  inf{i  >  0  :  Cf  >  4  ou  C(  <  -/i }  • 

On  defmit  la  prob?.biliic  /;+  (a:)  =  P<:„=  r({CT’  =  — 4}M+)-  En  utilisant  la  for- 
mule  de  Dynkin,  on  obtient  Tequation  differentielle  ordinaire  suivante: 

P+  +  74  =  0 

/4(-/i)  =  1  .  p+ih)  =  0. 

La  resolution  de  cette  equation  nous  donne  la  probabilite  de  “refuser  A+  a  tort 
”,  p+  =  7>x(0)  (cf.  [Fleming-Rishelj(cliapV.th7.1]): 

1  -  c-'» 

-  cb  -  e-‘^  • 
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Par  un  raisonnemenfc  identique  sur  /4_,  on  obtient  I’expression  de  la  proba- 
bilite  d’erreur  p_  =  p_(0) : 


V-  = 


!>h  _  p—h 


De  meme,  on  calcule  les  temps  moyens  pour  prendre  une  decision  en  resolvant 
deux  equations  difFerentielles  ordinaires  (cf.  [Fleming-Rishel][chapV.  th7.1]). 
Plus  precisement  E(2’^)  est  donne  par  la  solution  de  I’equation 

.  - 2n - - 2n - 5+ +  1  =  0 

.  (J+[-h)  =  9+{l2)  =  0 

au  point  0. 

On  obtient  E(T'*)  de  fagon  similaire: 

-2B+(2  +  AtB+) 


E(r;)  = 
E(r)  = 


-2Z?_(2  +  Af5_) 

[B^-B.f 


[h  —  P+{k  +  4)1 , 

[h  ~  P-{h  +  h)]  ■ 


Remarque  5.3  Dans  les  applications  on  prend  generalement  =  p_  =  a. 
Ainsi  un  simple  calcu!  nous  donne 


/|  I2  =  In 


a 


a 


el  on  note  quo 


E(r:)«|^E(r;). 


Exemple  5.4  Pour  /?_  =  —  1,  =  —0.25  ct  o  =  0.05.  on  obtient  =z  = 

2.94.  Les  temps  moyens  sont  alors  E{Tl',  =  4.70  et  E{Tl)  =  18.79. 

Remarque  5.5  Pour  comparer  les  tcm()s  moyens  dcs  deux  procedures,  on  s’in- 
tcrcsse  au  signe  de  I’expression: 


2A^-(I-2o)ln 


1  -  a 


o 


avec  A  tel  que  — ;==  [  e  =  a. 

\/27  2-00 

On  verifie  numeriquement  que  cette  expression  est  positive,  VO  <  a  <  1/2. 
On  cn  deduit  que  pour  une  probabilite  d’erreur  donnee,  les  temps  moyens  de  la 
premiere  procedure  sont  plus  grands  que  ceux  de  la  deuxieme.  Cette  difference 
est  d’autant  plus  grande  que  a  est  petit. 
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En  preparation  de  I’etude  du  cas  e  >  0,  nous  reecrivons  la  statistique 
sous  une  autre  forme.  En  utilisant  I’egalite 

h  {vk+i  -  yk)  =  -r  ((yt+i  -  hY  -  yl+i  +  yl) 


on  obtient 

T  _  ~  /~2  ~2  \  ~  ^  •  \2 

•^n  ““  217^  3^10+1/  2  ^  ^  \3/A:+1  Vk} 

.2 

2  ^  ’ 

^  ‘  Jt:=to+l 

Lemme  5.6  Sur  A+  U  A^,  on  a  ^estimation: 

I'o+n 

E  At  =  c»(n/aI)  . 

i=IO  +  l 

Preuve  Demontrons,  par  exemple,  ce  resultat  sur  I’evenement  A+. 
Sur  A+,  on  a  : 

yk+i  =  ilk  +  At  B+  B+  Xk  +  \fKia+  H+  Uk , 


d’ou  I’expression: 


lO+n  t'o+n  lo+n 

E  (»+.-»)’  =  AiV/'Bi  E  4  +  E  “* 

fc=io  +  l  k=io+l  k=to  +  i 

lo+n 

+2  At^^^  III  a+ B+  E 

^■=*0+1 


On  obtient  I’estimation: 


lo+n  lo+n 

E  E  {yk+i-ykY\  A^  <At'^nlBl  E  E(x2)  +  Atnr^ 

t=iO  +  l  k=io  +  l 


L’egalite 


lo+n  io+n 

E  (  E  «*)(  E  =0t 

i=«o+l  A:=io+l 


entrame 


f 

(<0+n  \  2 

E  (yk+i-ykf~r^nAt]  ;A+ 

kstQ  +  l  j 


<  E  d  +  Hla+B+  £ 


XkUk 


*=io+l 


1'*°'"+'  J  iv»-.'.-n  /i=ini  /, 


-Af%=r-2Ai“„r’//jz?’  '£“  e(x2), 

*0  + 1 

On  remarque  que 


^  ^  “a)  ~  E{“a)  +  Ji(»  —  1)  =  7?(n  +  2) 

.  A=io  +  I  A=10+I 

puisque  {n;.}  est  un  bruit  blanc  standard.  Ainsi  on  a; 

Af^r'E  (  Y,  wt]  -  r-*  =  r-’ 


/.*— *0+1 


On  utilise  l-inegnme  (»  +  /,)»  <  2u=+2i^  pour  traitor  le  premier  lermede  (24)  et 
nnega  .tc  <le  Caucl.y  -  Sclitvara  pour  le  troisieme.  Vu  ,„e  les  moments  d’ordre 
-  et  4  do  n.  sont  lints  sur  I  inlervallc  |«.  (,|  ot  que  n  <  (t  -  a)JAt .  on  oblient  la 
majoration  suivante: 


(lO  +  n  \  2" 

H  (tft+i  -  -  r*  n  A<  I  <  c„.6  A/ . 

A=iO  +  l  J 

ou  c„.6  est  une  constantc  positive  dcpciidantc  de  Pintcrvalle  {a.6J  fj.ve. 
Sur  r-evenemeni  .4_.  la  demonstration  est  siinilaire. 


5.1.2  Le  cas  £  >  0 

On  dcfinit  la  statistique  L'  comme  il  suit: 

r'-  -  2 

.-V  Hi 


*2  pa  (yio+n+i  yl+i)  —  (Bi- —  B^) 
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(26) 


—e 


Bl-Bl 

2P 


io+n 

E  yl 


^.■=<0+1 


On  considere  At  =  e  ei  on  s’interesse  a  I’ecart  —  Z-„. 


Theoreme  5.7  Sur  les  evenements  C  et  C  definis  dans  la  section  3,  on  a 
I ’estimation 


L* -L,.  =  0(v/?). 


Preuve 

On  demontre  d’abord  que  —  L„  cst  d’ordre  0{y/£)  sur  U  /l_. 


L‘n-Ln  = 


-  B. 


+ 


2P 

B+-B^ 

2r2 


{ivl+n  -  2/.o+n)  -  (2/fo  +  l  -  yl  +  l)) 


lO+n 

E  (^^+1  -  n  e 

j_fc=(a+l 


r>2  p2  lo+n 

E  iyl-yl) 


^‘=*0+1 


L’egalite 

implique 


yl~yl  =  \/s  vk  (2  li{xk)  +  \/e  Vk) 


^ri  p2  (^io+n  +  1  /^^C'^i'c+n+l )  ^I'o  +  l  ^K^io+O) 


+ 


B+  ~  B. 
2r2 


10 +n 


{h{xk+i)  -  h{xk)f  -V^ne 

fc=  <0  + 1 


B+  ~B.,  ,  2  X  3/2  Bi  -  , 

+  ^  -  0-r2  +  I'.o+l)  -  -—p—  E 

A;=:»o  +  l 


2  Bl  -  B'^  2 


2r2 


A?=io  +  I 


Sachant  que  {u*}  est  un  bruit  blanc  gaussien  standard  et  que  tous  les  mo¬ 
ments  de  Xk  sont  finis  pour  k  <  m,  on  deduit  a  I’aide  du  lemme  5.6  que  Zr'  —  L„ 
est  d’ordre  0(\/?)  sur  .4+  U  A_. 
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Pour  raisonner  sur  C,  on  remarque  tout  simplement  qu’on  a  la  majoration 

+El(i‘„-i„f  iCn(A+UA.)1. 

De  I’inegalite  de  Cauchy-Schwarz,  on  obtient 

E[  (L'  -  L„Y  ;  C)  <  P[C  n  {A^  U  A.yy/^  +  Ce. 

Du  fait  que  E((iy*  —  Z/n)']  <  C  et  de  la  proposition  3.1,  on  deduit  le  resultat 
sur  Pevenement  C. 

De  meme,  du  theoreme  3.6,  on  deduit  le  resultat  sur  Pevenement  C. 

m 

Etant  donne  que  les  homes  detenninees  dans  le  paragraphe  5.1.1  sont  d’ordre 
0(1).  il  est  done  legitime,  d’apres  le  resultat  precedent,  de  les  utiliser  pour  tester 
Lf, 

D’apres  la  remarque  5.5,  nous  avons  choisi  d’implemcnter  le  test  selon  sa 
deuxieme  version. 


5.2  Test  du  rapport  de  vraisernblance  sur  les  sorties  des 
filtres  de  Kalman  -  Bucy 


On  montre  que  la  statistique  ulili.see  dans  Ic  paragraphe  4.2  sous  Phypothese 
(IIDl),  permet  encore  de  decider  entre  les  alternatives  “/I4.”  et  “.4_”  sous 
Phypothese  {H D2). 

Soil  /'o  <  <  >0  +  ni .  On  a 


+n 


i.  =  -  E  -  TT  E  (("+('  +  -  {«-(!  +  . 

^  K=n 

Zk  =  //+(!  -fe/i+)Xfc  -  //-(I  +cB.)3\  . 

On  rappclle  la  definition  du  terme  de  precision  du  tost: 


1  U+n 

«.  =  7  E  4 


k-tt 


.Sur  /-l+j  a  parlir  de  Pexpression  recursive  de  Zk  (equation  (16)),  on  obtient; 


E[Zl, ;  /4]  <  (1  +  erj..yi\  -  H.KJ^E[Zl ;  >1+] 

+e2(]  +  sB^y(B+  -  B.yjilE[x+'^ :  A+] . 

+£  [11+  /f + ( 1  +  e i?+ )  -  //_  A'_  ( 1  +  £  ))^t?+ 

+2£(1  +eB+){l+eB..)\B+-  B^\{1  -  ILK.)Il+E[Zkxt-,A+]. 
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On  utilise  I’inegalite 
E[Zkxt]A+]  < 


2e(H- eJ5_)(l  +  £B+)1B+  - 

£{\ £B^)\B+ —  B-\H^ 

2HZk1 


E[Zl ;  A^] 

•  E[x+2 ;  A+] . 


Sous  {HD2),  la  difference  H+K+  —  H_K_  est  d’ordre  0(e). 

On  utilise  le  fait  que 

0  <  (1  +  eB+)^l  -  H.K_f  +  <  1  -  H.K.  +  Ce  <  1. 

Sous  I’hypothese  (7/4),  on  obtient,  de  la  meme  fa^on  que  dans  le  paragraphe 
4.2, 

EfZ^  j  A+j  ^  A.).  +  C  , 


oil 


15+1(2  + ^5+) 


Sur  A_,  on  fait  un  raisonnement  similaire. 

Remarque  5.8  Le  processus  (Zkj  est  d’ordre  0(v/e)  sous  I’hypothese  (77771) 
et  d’ordre  0(e)  sous  I’hypothese  (77/72). 

De  meme  que  pour  le  test  precedent,  on  considere  deux  types  de  procedures. 


Premiere  procedure: 


Pour  r  fixe,  on  introduit  le  temps  d’arret 


A'*  =  inffjj  >1  :  R„  >  r}  A  ?n  . 

On  calcule  L„  jusqu’a  ce  qu’on  ait  7?„  >  r.  Par  un  raisonnement  similaire  a 
cclui  du  paragraphe  5.1,  on  choisit  r  =  max(r_,r+)  avec 

r_  =  4A^7-  et  r+  =  4A^i?+. 


La  constante  A  est  obtenue  par  inversion  de  la  fonction  de  repartition  de  la  loi 
A^(0, 1)  calculee  pour  la  probabilite  a.  Les  estimations  des  temps  moyens  pour 
prendre  une  decision  sont  donnees  par 


eN’ 


eNl 


,-B..(2  +  eB.) 
(7?+-5_)2  ’ 

^  ^2~5+(2  +  eB+) 

HTTsIp  '- 
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} 


i 

I 


Remarque  5.9  Les  temps  moyens  pour  prendre  une  decision  associes  a  cette 
procedure  sont  identiques  a  ceux  du  test  precedent  quand  on  considere  le  meme 
type  d’approche  (cf.  paragraphe  5.1). 

Deuxieme  procedure: 

Soient  /j  >  0  et  /2  >  0  des  homes  fixees. 

On  introduit  le  temps  d’arret  N’ : 

N”  =  inf{n  >  i]  :  Ln  >  h  ou  <  —li]  A  m . 


Raisonnons  sur  I’evenement  par  exemple. 

Nous  savons,  d’apres  ce  qui  precede,  que  sous  I’hypothese  (//4)  le  processus 
=  ZkVk+i  est  approximativement  un  processus  melangeant.  Par  un  raison- 
nement  similaire  a  celui  du  paragraphe  4.2,  pour  i  =  en,  on  approche  Ln  par  le 
processus  de  diffusion  {Ct}  defini  par: 

dQ  =  A+/2  dt  +  \/i)+k^  dW+  . 

On  obtient  la  probabilite  d’erreur  “probabilite  de  refuser  A+  a  tort”: 

_  1  - 

e’Hh  _  g->;+b  ’ 

oil  ?/+  =  I/O'*'  ainsi  que  le  temps  moyen  pour  prendre  une  decision  sur  A+  : 


Eir;) 


h  _  {h  +  h) 

—  />+  - 

A+  A+ 


De  meme,  on  obtient  la  probabilite  d’erreur  p_.  “probabilite  de  refuser  A. 
a  tort”: 

_  1  -  c-’*-'' 

~  (’t-h  _  ’ 

oil =  l/0~  et  le  temps  moyen  pour  prendre  une  decision  sur  : 


E(r)  =  2 


/. 


-P- 


{L  +  L) 

A_  . 


Remarque  5.10  Sous  {HDl)  les  temps  d’attente  moyens  sont  d’ordre  0(e) 
tandis  que  sous  (1102)  ils  sont  d’ordre  0(1).  On  pent  s’attendre  done  a  ce  que, 
dans  ce  deuxieme  cas,  il  y  ait  des  intervalles  ou  aucune  decision  ne  sera  prise. 
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Remarque  5.11  Sous  I’hypothese  {HD2),  I’ecart  est  d’ordre  0(e), 

done  les  bornes  li  et  I2  sont  “tres  proches”.  On  a  une  situation  similaire  a  celle 
decrite  dans  la  remarque  5.5,  i.e.  les  temps  moyens  d’attente  sont  plus  longs 
dans  la  premiere  procedure  que  dans  la  deuxieme. 

Le  test  du  rapport  de  vraisemblance  sur  les  sorties  des  filtres  de  Kalman  sera 
mis  en  oeuvre  selon  la  deuxieme  procedure. 


6  Application 

Nous  rappelons  que  notre  objectif  est  d’etudier  le  comportement  des  filtres  pro¬ 
poses  par  rapport  a  la  solution  optimale  du  probleme  de  filtrage  et  de  comparer 
les  differentes  procedures  permettant  de  les  obtenir.  Nous  devons  done,  en  pre¬ 
mier  lieu,  definir  des  criteres  de  comparaison. 

6.1  Les  criteres  de  comparaison 

Une  premiere  idee  est  de  considerer,  pour  une  probabilite  d’erreur  fixee, 

•  I’amplitude  moyenne  de?  intervalles  detectes,  pour  comparer  les  tests  de 
detection 

•  le  temps  moyen  necessaire  pour  prendre  une  decision,  pour  comparer  les 
tests  de  determination  du  signe. 

Pour  comparer  les  tests  en  termes  de  probabilite  d’erreur,  il  semble  naturel  de 
s’inieresser  aux  pourcentages  suivanis: 

•  pourcentage  d’intorvalles  detertes  corrects. 

•  pourcentage  de  decisions  correcles  (i,e,  nombre  de  decisions  concctes  par 
rapport  au  nombie  de  decisions  prises). 

Nous  qualifions  une  decision  clc  correcle  lorsque  le  signe  choisi  correspond  a 
celui  de  la  trajcctoire.  ccci  sur  un  intervalle  de  monotonie  correct. 

Les  temps  moyens  j)our  atleindre  une  decision  sont  calcules  d’abord  theori- 
quement  (voir  les  paragraphes  4  et  5)  et  sont  ensuite  confrontes  aux  resultats 
obtenus  par  simulation. 


6.2  Bxempies 

L’etude  du  probleme  sous  I’hypothcse  {If  D2)  etant  notre  principal  c^entre  d’inte- 
ret,  nous  consideron-s  des  excmples  verifiant  cette  hypothesc  ct  nous  appliquons 
les  difFerents  tests  decrits  dans  le  paragraphe  5.  Nous  etudions  I’influence  de  la 
fonction  cl’observation  sur  I’efficacite  des  tests,  dans  i’exempie  6.2  et  celle  de  la 
derive  dans  I’exemple  6. .3.  Nous  nous  intercssons  ensuite  a  un  exemple  verifiant 
I’liypothese  (HDl),  presentc  dans  [FJSZj.  Ceci  nous  permettra  d’illustrer  que 
dans  ce  cas,  les  decisions  sont  prises  bcaucoup  plus  rapidement  que  sous  {HD2). 
Nous  justifierons  ce  fait  par  I’etude  du  comportement  de  la  solution  de  I’equation 
de  Zakai. 
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Dans  ie  paragraphe  5,  sous  {HD2),  nous  avons  propose  une  methode  de 
determination  des  homes  pour  les  tests  du  signe  et  d’estimation  des  temps 
d’atteiite  sous  I’hypothese  {H4):  <  0  et  <  0.  II  est  interessant  d’etudier 

les  resultats  obtenus  dans  le  cas  ou  I’hypothese  {H4)  n’est  pas  verifiee.  C’est  ce 
que  nous  nous  proposons  de  faire  dans  I’exemple  6.6. 

Pour  chaque  exemple,  nous  commentons  les  resultats  statistiques  sur  le  com- 
portement  des  tests.  On  pent  consulter  ces  resultats  dans  I’annexe  A. 

Remarque  6.1  On  rappelle  que  la  determination  des  homes  utilisees  dans  les 
tests  de  decision  sur  le  signe  repose  sur  la  possibilite  d’atteindre  une  situation 
stationnaire.  Une  fois  les  passages  a  zero  detectes,  nous  devons  attendre  un 
“certain”  temps,  no,  avant  de  cumuler  les  satistiques  de  decision  des  tests.  Ce 
temps  est  cependant  assez  court  puisque,  comme  nous  I’avons  deja  remarque,  a 
I’interieur  d’un  intervalle  de  monotonie,  le  systeme  a  tendence  a  accompagner 
rapidement  revolution  du  systeme  lineaire  associe.  Nous  avons  choisi  pour  ces 
exemples  no  =  6. 

6.2.1  Exemples  sous  I’hypothese  (HD2) 

Exemple  6.1  On  considere  le  systeme 

j  Xk+i  =  Xk  +  0.01  b{xk)  +  \/0.0T Uk  ,  To  ~  A/'(-5,0.1) 

1  Vk  =  l-t;!,-!  +  \/0.01  Vk , 

avec 

bix)  =  ° 

[  —0.25  a: .  t  >  0  . 

On  s’interesse  tout  d’abord  a  une  trajectoire  simulw  sur  I’intervalle  de  temps 
(0. 10). 

Dans  ce  cas,  les  resultats  des  difFerentes  procedures  etant  similaires,  nous 
illustrons  par  la  figure  1,  uniquement  ceux  des  tests  de  detection  et  du  rapport 
de  vraisemblance  sur  les  observations. 

Dans  cette  figure,  la  trajectoire  du  signal  a  estimer  est  tracee  en  trait  plein  et 
le  filtre  propose  en  pointille.  Les  intervalles  de  monotonie  detectes  sont  signales 
en  noir  dans  le  cas  ou  une  decision  est  prise  sur  le  signe  et  en  gris  s’il  n’y  a  pas  de 
decision.  Nous  rappelons  que  sur  un  intervalle  de  monotonie,  si  le  test  applique 
decide  “signe  positif”,  (respectivement  “signe  negatif”)  on  trace  la  sortie  du 
filtre  de  Kalman-Bucy  (7)  (respectivement  la  sortie  du  filtre  (8)).  Si  le  test  ne 
prend  pas  de  decision,  les  sorties  des  deux  filtres  sont  representees. 
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II  est  interessant  de  regarder  le  comportement  de  la  solution  optimale  pour 
ce  probleme.  Dans  la  figure  2  on  represente,  en  plus  de  la  trajectoire  (en  trait 
plein),  I’estimation  optimale  obtenue  par  resolution  numerique  de  I’equation  de 
Zakai  (en  trait  pointille),  entouree  par  sa  region  de  confiance  correspondante  a 
une  probabilite  de  95%.  Un  generateur  automatique  de  programmes  FORTRAN 
a  ete  utilise  pour  obtenir  cette  solution. 

Dans  les  figures  3  a  6,  on  trace  la  densite  conditionnelle  a  divers  instants 
de  I’intervalle  [0,10].  Le  trait  vertical  correspond  a  la  valeur  du  signal  simule 
a  I’instant  t.  On  note  qu’au  moment  ou  la  decision  sur  le  signe  est  prise,  la  loi 
conditionnelle  approche  une  ioi  gaussienne  et  que  pres  d’un  passage  a  zero,  la 
loi  conditionnelle  a  deux  “bosses”. 

On  pent  constater,  par  exerr  pie,  qu’a  Tinstant  t  =  6.5  les  tests,  ayant  deja 
detecte  un  intervalle  de  monotonie.  ne  sont  pas  encore  capables  de  decider  sur 
le  signe.  Cette  decision  sera  pri.'e  aux  environs  de  <  =  8.48.  En  regardant  le 
comportement  de  la  solution  de  I’equation  de  Zakai,  on  remarque  I’existenco 
de  deux  “bosses”  a  partir  de  I’instant  t  •—  3,  approximativement.  A  I’instant 
t  =  8.48,  instant  oii  les  tests  de  decision  sn/  le  signe  sont  capables  de  prendre 
une  decision,  une  des  “bosses”  de  la  densite  conditionnelle  est  deja  negligeable 
(voir  la  figure  5),  ce  qui  explique  le  fait  qu’une  decision  puisse  etre  prise.  A  cet 
instant,  le  filtre  propose  et  I’estimation  optimale  ont  un  comportement  similaire. 

On  s’interesse  maintenant  aux  resultats  numeriques  obtenus  pour  une  tra¬ 
jectoire  'imulee  sur  I’intervalle  de  temps  (0, 100)  (voir  les  tableaux  1  et  2). 

On  constate  que 

•  Le  test  de  detection  sur  le  sorties  du  filtre  de  Kalman-Bucy  permet  de 
detector  des  intcrvalles  de  monotonie  plus  longs  que  le  test  sur  les  obser¬ 
vations;  cependant  I’encur  commise  c'^t  aussi  plus  grande. 

•  Le  fait  d’augmenter  la  borne  de  detection  (par  exemple  en  diminuant 
la  probabilite  d’erreur  a)  entrainc  evidemment  la  detection  d'inlervalles 
dc  monotonie  plus  petits.  II  semble  que  nous  controlons  moins  bien  la 
probabilite  d’eneur  du  test  sur  la  sortie  des  filtres  de  Kalman-Bucy  que 
celle  du  test  sur  les  observations.  Voyons  par  exemple  ce  qui  se  passe 
quand  on  prend  aj  =  5%  (voir  le  tableau  2).  Le  pourcentage  d’erreur  de 
detection  est  superieur  a  celui  qui  avait  etc  impose. 

•  Les  resultats  des  tests  de  decision  sur  le  signe  bases  sur  le  rapport  de 
vraisemblance  sur  les  observations  ou  sur  les  sorties  du  filtre  de  Kalman 
sont  comparables.  Toutefois  le  temps  moyen  pour  prendre  une  decision 
semble  etre  legereraent  plus  petit  dans  le  premier  cas  que  dans  le  deuxieme. 
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•  Le  fait  de  diminuer  la  probabilite  d’erreur  pour  la  decision  sur  le  signe 
entraine  naturellement  une  reduction  du  pourcentago  d’intervalles  avec 
decision  et  des  temps  d’attente  plus  grands  pour  prendi  ''  cef  decisions. 

Remarque  6.1  Les  temps  moyens  pour  prendre  une  decision  calcules  tbeori- 
quement  sont  beaucoup  plus  eleves  que  ceux  verifies  en  pratique. 

Exemple  6.2  On  modifie  la  foncMon  d’observation  dans  I’exemple  6.1.  On  con- 
s,  .iere: 

Ii{x)  —  0.4  |a;| . 

On  obtient  les  resultats  presentes  dans  le  tableau  3. 

Le  fait  de  diminuer  le  coefficient  \H^\  —  |//_|  semble  avoir  pour  effet  I’aug- 
nientation  des  pourcentages  de  detection  d’intervalles  corrects.  Autrernent  dit, 
lorsqu'on  diminue  le  coefficient  [Lf+j  =  |//_|,  les  constantes  des  tests  de  detection 
augmentent.  On  detecte  plus  dc  pa.ssages  a  zero  el  les  longueurs  des  intervalles 
de  momotonie  diminucnt.  Ce  comportement  des  tests  est  satisfaisant  vu  qu’une 
diminution  du  coefficient  [/f+j  =  |//_|  est  equivalente  a  une  augmentation  du 
bruit  d ’observation. 

On  verifie  %alement  une  “legere”  diminution  des  temps  d’attente  pour  pren¬ 
dre  une  decision.  Ceci  s’explique  du  fait  que  les  statistiques  des  tests  du  signe 
sont  calcules  loin  des  passages  a  zero. 

Dans  ce  cas,  le  test  de  detection  des  passage,s  a  zero  base  sur  les  sorties  d’un 
des  filtres  de  Kalman-Bucy  est  meilleur  que  cclui  base  sur  les  observations. 

Exemple  6.3  On  considere  les  memes  parametres  que  dans  I’exemple  prece¬ 
dent  sauf  pour  la  derive  b(x)  : 


b{x) 


—5x ,  X  <0 
—0.25  a: ,  i  >  0  . 


Les  resultats  sont  a  consulter  dans  le  tableau  4. 


On  remarque  que  le  pourcentage  de  decisions  est  maintenant  plus  eleve  et 
que  le  temps  moyen  pour  prendre  une  decision  diminue  considerablement.  Ceci 
illustre  parfaitement  les  calculs  theoriques  pr&entes  dans  le  paragraphe  5.  Le 
fait  de  diminuer  le  rapport  \B+\/{B+  -  j3_)^  (respectivement  \B^\/{B^  -  R_)^) 
a  pour  effet  d’augmenter  le  terme  de  precision  dans  les  statistiques  L  et  L.  On 
obtient  done  de  meilleurs  resultats,  notamment  pour  ie  temps  d’attente  “cote 
positif”  (respectivement  “cote  negatif”). 
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6.2.2  Exemple  sous  I’hypothese  (HDl) 

Exemple  6.4  Nous  changeons  la  fonction  d’observation  de  rexeniple  6.1: 

—X ,  X  <0  , 

2x ,  X  >0 . 

Le  systeme  considerc-  verifie  I’hypothese  (HDl). 

Par  les  figures  7  a  10,  nous  illustrons  les  resultats  des  difFerentes  procedures 
appliquees  a  une  trajectoire  simulee  sur  I’intervalle  de  temps  [0, 10].  Nous  re- 
marquons  que  le  test  de  la  variation  quadratique  decide  plus  souvent  que  le  test 
sur  les  sorties  des  filtres  de  Kalman-Bucy. 

Le  comportement  de  I’estimation  optimale  et  sa  region  de  confiarice  a  95% 
representes  dans  la  figure  11,  montrent  clairement  que  dans  ce  cas,  la  loi  condi- 
tionnelle  approche  beaucoup  plus  rapidentent  une  loi  gaussieniie  apres  un  pas¬ 
sage  a  zero  que  dans  I’exemple  6.1  (voir  la  figure  2).  Ceci  justifie  le  fait  que  nous 
decidons  du  signe  de  Xk  plus  souvent  sous  I’hypothese  {HDl)  que  sous  {HD2). 

La  comparaison  entre  les  figures  12  a  14  et  les  figures  3  a6  (voir  I’exemple 6.1) 
illustre  a  nouveau  la  remarque  precedente.  Dans  ce  cas,  une  des  deux  “bosses”  de 
la  loi  conditionnelle  s’estompe  rapidement  aux  environs  des  instants  de  decision 
(t  =  6.86  et  t  =  6.88),  tandis  que  dans  I’exemple  6.1  on  ne  retrouve  un  tel 
comportement  qu’aux  alentours  de  I’instant  t  —  10. 

Les  resultats  numeriques  des  tests  appliques  a  une  trajectoire  simulee  sur 
I’intervalle  de  temps  {0, 100]  (voir  le  tableau  5)  mettent  egalement  en  evidence 
que  les  temps  d’attente  pour  prendre  une  decision  sous  I’hypothese  {HDl)  sont 
ncttement  plus  courts  que  sous  I’hypothese  {HD2). 

Remarque  6.5  Quelques  remarques  sur  les  resultats  des  procedures  appliquees 
a  plusieurs  trajectoires  simulees  sur  un  intervalle  de  temps  [0, 100]. 

•  Les  tests  detectent  moins  de  passages  a  zero  lorsque  I’ecart  \H^  —  H^\ 
augmente  et  le  taux  d’erreur  est  superieur  au  taux  souhaite.  De  fa?on 
plus  nette  que  sous  {HD2),  nous  controlons  “moins  bien”  la  probabilite 
d’erreur  du  test  de  detection  sur  la  sortie  d’un  des  filtres  de  Kalman-Bucy 
que  celle  du  test  sur  les  observations. 

•  Pour  une  difference  \(7+H+  —  =  0.5,  les  tests  de  determination  du 

signe  ont  des  taux  d’erreur  superieurs  a  ceux  demandes.  Le  test  de  la 
variation  quadratique  a  tendance  a  decider  plus  souvent  que  le  test  du 
rapport  de  vraisemblance  sur  les  sorties  des  filtres  de  Kalman-Bucy,  mais 
semble  generalement  moins  fiable.  (rf,  [F.ISZ]). 
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•  En  augmentant  I’ecart  —  les  tests  de  determination  du  signe 

decident  plus  souvent  et  prennent  moins  de  fausses  decisions. 


6.2.3  Un  exemple  sous  I’hypothese  {HD2)  ne  verifiant  pas  (^^4) 
Exemple  6.6  On  considere  le  systeme  suivant: 


^k+i  =  Xk  +  0.01  b{xk)  +  \/0.01  Uk  ,  xo  ~  A/'(-5,0.1) 
Vk  =  |xfc|  4-  \/0.01  Vk, 


avec 


b(x) 


—X  X  <  0 
0.2x  X  >  0 


Remarque  6.7  Le  coefficient  de  derive  B+  etant  positif  et  etant  negatif, 
il  y  a  “peu’'  de  passages  a  zero  et  le  signal  a  tendance  a  “fuir”  vers  I’infini  des 
qu’il  est  du  cote  positif. 


Nous  etudions  un  cas  ou  1’  hypothese  {HD2)  est  verifi^  mais  pas  (//4). 
Cependant  on  applique  les  tests  presentes  dans  la  section  5.  On  rappelle  que 
les  bornes  utilisces  dans  ces  tests  ont  ete  determinees  sous  I’hypothese  {Hi)  et 
ne  sont  plus  justifies  dans  cet  exemple.  II  est  done  interessant  de  regarder  les 
resultats  qu’elles  engendrent  dans  une  telle  situation. 

Dans  la  figure  15,  nous  tra^ons  le  filtre  approche  construit  a  partir  du  test 
de  detection  et  du  test  de  rapport  de  vraisemblance  bases  sur  les  observations. 
Les  r^ultats  obtenus  par  les  autres  procedures  sont  similaires. 

En  comparant  le  filtre  propose  et  I’estimation  optimale  (voir  la  figure  16),  on 
remarque  que  leurs  comportements  sont  semblables  sur  I’intervalle  de  monotonie 
detecte  [0,3]  et  aussi  a  partir  de  I’instant  t  =  7,  environ.  Ceci  est  satisfaisant  vu 
que  sur  I’intervalle  de  monotonie  (5.80, 10],  les  instants  de  decision  sur  le  signe 
sont  t  =  7.15,  pour  le  test  sur  les  observations,  et  <  =  7.20  pour  le  test  sur  les 
sorties  des  filtres  de  Kalman-Bucy. 

On  se  permet  de  conclure  qu’au  moins  pour  certains  cas  ne  verifiant  pas 
I’hypothese  (Hi),  on  peut  encore  utiliser  les  bornes  calculees  sous  cette  hy¬ 
pothese  et  obtenir  des  resultats  “corrects”. 


Remarque  6.8  Dans  cette  situation  ou  I’hypothese  {Hi)  n’est  pas  verifiee, 
nous  ne  savons  pas  estimer  theoriquement  le  temps  d’attente  d’une  decision. 
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7  Conclusion 


On  a  decrit  des  tests  permettant  de  determiner  les  intervalles  de  linearite  de  la 
fonction  d ’observation,  pour  un  systeme  lineaire  par  morceaux  en  temps  discret 
avec  petit  bruit  d ’observation.  Deux  types  de  tests  sont  necessaires:  un  premier 
test  permet  de  detecter  les  intervalles  ou  il  n’y  a  pas  de  passages  a  zero;  un 
deuxieme  test  permet  de  decider  du  signe  pris  par  la  trajectoire  sur  cet  intervalle. 
Sous  I’hypothese  (HD2),  difFerents  tests  on  ete  proposes.  Pour  la  detection  de 
I’intervalle  de  monotonie,  on  a  le  choix  entre  un  test  base  sur  les  observations 
et  un  test  base  sur  les  sorties  d’un  des  filtres  de  Kalman-Bucy.  Pour  I’etape 
de  decision  du  signe,  on  decrit  un  test  de  rapport  de  vraisemblance  sur  les 
observations  et  un  autre  test  de  rapport  de  vraisemblance  sur  les  sorties  des 
filtres  de  Kalman-Bucy. 

L’efficacite  des  tests  de  detection  depend  de  la  determination  de  la  valeur 
d’une  borne,  le  “cutoff”,  line  formule  permettant  d’obtenir  cette  valeur  a  ete 
justifiee  theoriquement  pour  le  test  base  sur  les  observations.  Ceci  n’a  pas  ete 
de  meme  pour  le  test  base  sur  les  sorties  du  filtre  de  Kalman-Bucy.  11  nous  sem- 
blc  done  naturel  que  ce  deuxieme  test  se  comporte  moins  bien  que  le  premier. 
Quant  aux  tests  de  determination  du  signe,  on  a  verifie  que  leur  comnortement 
est  semblable.  Des  formules  permettant  de  determiner  les  bornes  et  les  temps 
moyens  d’attente  pour  prendre  une  decision  ont  ete  obtenues  theoriquement 
sous  I’hypothese  (H4).  Les  probabilites  d’erreur  fixees  sont  a  comparer  avec 
celles  obtenues  sur  le  probleme  simule  ainsi  que  les  temps  moyens  d’attente. 
On  remarque  de  maniere  generale,  que  les  temps  moyens  d’attente  calcules  the¬ 
oriquement  sont  plus  longs  que  ceux  obtenus  en  pratique.  On  note  egalement 
que  I’ordre  de  grandeur  des  temps  moyens  d’attente  theoriques  et  empiriques  est 
plus  grand  si  on  se  place  sous  I’hypothese  (HD2)  que  sous  I’hypothese  (HDl). 

On  peut  envisager  d’etendre  I’etude  ici  presentee  a  des  systemes  oil  les  fonc- 
tions  b,  a  ei  h  possedent  /  points  critiques,  avec  /  >  1.  On  peut  aussi  considerer 
des  systemes  de  dimension  n  >  1.  avec  des  bruits  correles. 

Une  extension  au  cas  de  fonctions  monotones  par  morceaux  existe  dans 
[Fleming- Pardou,x],  sous  I’hypothese  (HDl).  Des  essais  ont  ete  faits  par  [FJSZ] 
dans  le  cas  d’une  dynamique  non  lineaire.  II  reste  a  justifier  son  analogue  sous 
I’hypothese  (HD2). 
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A  Resultats  numeriques  et  representations  gra- 
phiques  des  examples  present es  dans  la  sec¬ 
tion  6 

A.l  Resultats  num4riques 

Nous  introduisons  les  notations  suivantes: 

TDO:  test  de  detection  sur  les  observations. 

TDK:  test  de  detection  sur  la  sortie  d’un  des  filtres  de  Kalman-Bucy. 

TLO:  test  du  rapport  de  vraisemblance  sur  les  observations. 

TLK:  test  du  rapport  de  vraisemblance  sur  les  sorties  des  filtres  de  Kalman 
-Bucy. 

ad :  niveau  d’erreur  pour  le  test  de  detection. 

ofj :  niveau  d’erreur  pour  le  test  de  determination  du  signe. 

C :  “cutoff”  du  test  de  detection. 

—h  ,  I2 :  intervalle  de  confiance  du  test  de  determination  du  signe. 

ET* ,  ET^  :  temps  moyens  theoriques  d’attente  d’une  decision. 

Pi :  pourcentage  de  pas  de  temps  detectes. 

P2 :  pourcentage  d’intervalles  de  monotonie  corrects. 

P3 :  pourcentage  de  pas  de  temps  avec  decision. 

P4 :  pourcentage  de  decisions  correctes. 
fl  moyennes  empiriques  des  temps  d’attente  d’une  decision. 


Test  de  detection 

TDO 

TDO 

TDK 

TDK 

Test  du  signe 

TLO 

TLK 

TLO 

TLK 

C 

0.143 

0.143 

0.143 

0.143 

—li  I2 

-2.94  2.94 

-7.70  7.67 

-2.94  2.94 

-7.70  7.67 

E7’:  et; 

18.8  4.71 

18.9  4.69 

18.8  4.71 

18.9  4.69 

pi 

84.5% 

84.5% 

86.3% 

86.3% 

Pi 

96.3% 

96.3% 

93.2% 

93.2% 

P3 

40.6% 

40.6% 

41.7% 

41.7% 

Pi 

100% 

100% 

100% 

100% 

1.3  1.5 

1.3  1.5 

2.1  1.6 

1.7  1.6 

Tableau  1:  resultats  numeriques  pour  I’exemple  6.1  (a^  =  5%,  a,  =  5%) 


Test  de  detection 

TDO 

TDO 

TDK 

TDK 

Test  du  signe 

TLO 

TLK 

TLO 

TLK 

C 

0.17 

0.17 

0.17 

0.17 

—h  h 

-2.94  2.94 

-7.70  7.67 

-2.94  2.94 

-7.70  7.67 

et:  et; 

18.8  4.71 

18.9  4.69 

18.8  4.71 

18.9  4.69 

Pi 

82.4% 

82.4% 

84.2% 

84.2% 

P2 

97.9% 

97.9% 

97.5% 

97.5% 

P3 

39.0% 

39.0% 

41.4% 

41.4% 

Pi 

100% 

100% 

100% 

100% 

f_f+ 

1.3  1.5 

1.3  1.5 

2.1  1.5 

1.7  1.5 

Tableau  2:  resultats  numeriques  pour  I’exemple  6.1  (aj  =  2.5%,  a,  =  5%) 


Test  de  detection 

TDO 

TDO 

TDK 

TDK 

Test  du  signe 

TLO 

TLK 

TLO 

TLK 

C 

0.121 

0.121 

0.302 

0.302 

—ii  I2 

-2.94  2.94 

-4.38  4.36 

-2.94  2.94 

-4.38  4.36 

et:  et; 

18.8  4.71 

18.9  4.70 

18.8  4.71 

18.9  4.70 

Pi 

70.8% 

70.8% 

75.4% 

75.4% 

P2 

98.4% 

98.4% 

98.8% 

98.8% 

P3 

35.0% 

34.9% 

40.1% 

38.8% 

Pi 

100% 

100% 

100% 

100% 

f _  f+ 

0.0  1.1 

0.05  1.2 

0.0  1.4 

0.05  1.5 

Tableau  3:  resultats  numeriques  pour  I’exemple  6.2  {aj  =  5%,  or,  =  5%) 


Test  de  detection 

tdo 

TDO 

TDK 

TDK 

Test  du  signe 

tlo 

tlk 

TLO 

TLK 

C 

0.144 

0.144 

0.144 

0.144 

— /i 

-2.94  2.94 

-7.69  7.54 

-2.94  2  94 

-7.69  7.54 

et:  et* 

2.350.117 

2.400.115 

2.35  0.317 

2.400.115 

pi 

85.5% 

85.5% 

87-0% 

87.0% 

P2 

95.3% 

95.3% 

92.5% 

92.5% 

P3 

66.8% 

66.8% 

69.3% 

68.4% 

P4 

95.8% 

100.0% 

90.0% 

94.7% 

T^f+ 

0.0  0.28 

0.12  0.29 

0.0  0.30 

0.12  0.33 

Tableau  4:  resultats  numeriques  pour  Texemple  6.3  (aj  =  5%,  a,  =  5%) 


Test  de  detection 

TDO 

TDO 

TDK 

TDK 

Test  du  signe 

TLO 

TLK 

TLO 

TLK 

C 

0.165 

0.165 

0.165 

0.165 

—li  I2 

-4.70  2.37 

-15.80  7.80 

-4.70  2.37 

-15.80  7.80 

et:  et; 

0.44  0.13 

2.46  0.44 

0.44  0.13 

2.46  0.44 

Pi 

86.6% 

86.6% 

88.2% 

88.2% 

P2 

93.5% 

93.5% 

87.5% 

87.5% 

P3 

71.2% 

63.4% 

76.1% 

65.4% 

P\ 

82.4% 

87.0% 

71.1% 

66.7% 

T.7\ 

0.38  0.29 

1.60  0.28 

0.40  0.37 

1.60  0.30 

Tableau  5:  resultats  numeriques  pour  I’exemple  6.4  (aj  =  5%,  a,  =  5%) 
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A. 2  Representations  graphiques 


I 

i 


\ 


(testdobs-testlrobs) 


signal  -f  estimate 


Figure  1:  Sous  {HD2),  exemple  6.1:  tests  de  detection  des  passages  a  zero  et  de 
decision  du  signe  bases  sur  les  observations. 
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5 


signal  estimate 


Figure  2:  Sous  {HD2),  exemple  6.1;  solution  dc  I’equation  cle  Zakai  (le  niveau 
cle  gris  correspond  a  une  region  de  confiance  a  95%). 
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Figure  3;  Sous  (HD2),  exemple  6.1:  densite  conditionnelle  a  I’instant  t  =  6.5 


jJA. 


Figure  4:  Sous  {HD2),  exemple  6.1:  densite  conditionnelle  a  I’instant  t  =  7 
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Figure  5;  Sous  {HD2),  exemple  6.1:  densite  conditionnelle  a  I’instant  t  =  8.48 


Figure  6:  Sous  {HD2),  e.xemple  6.1:  densite  conditionnelle  a  I’instant  t  =  10 
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Figure  7:  Sous  (///)!),  exemple  6.4:  test  <Je  detection  des  passages  a  zero  base 
sur  les  observations  et  test  de  decision  du  signe  basesur  la  variation  quadratique. 
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(testdobs-testirko)  signal  +  estimate 


Figure  8:  Sous  (HDl),  exemple  6.4:  test  de  detection  des  passages  a  zero  bcise 
sur  lea  observations  et  test  de  decision  du  signe  base  sur  la  sortie  des  filtres  de 
Kalman-Bucy. 


(tcstdko-testqv  ) 


signal  +  estimate 


(testdko-testirko)  signal  +  estimate 


Figure  10:  Sous  (JIDl),  exemple  6.4:  tests  de  detection  des  passages  a  zero  et 
de  decision  du  signe  bas&  sur  ia  sortie  des  filtres  de  Kalman-Bucy. 
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Figure  IJ:  Sous  {HDl),  exemple  6.4:  solution  de  1 ’equation  de  Zakai  (le  niveau 
de  gris  correspond  a  une  region  de  confiance  a  95%). 
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Figure  12;  Sous  (HDl),  exemple  6.4:  densite  conditionnelle  a  I’instant  t  =  6.60 


Figure  13:  Sous  (HDl),  exemple  6.4:  densite  conditionnelle  a  I’instant  t  =  6.86 
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Figure  14:  Sous  {HD\)^  example  6,4:  densite  conditionnelle  a  I’instant  t  =  6.88 
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Figure  15:  Exemple  6.6:  test  de  detection  des  passages  a  zero  base  sur  les  ob¬ 
servations  et  test  de  decision  du  signe  base  sur  la  variation  quadratique. 
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Filtres  approches  pour  un 
probleme  de  filtrage  non  lineaire 
d’une  diffusion  de  dimension  2 
mesuree  par  un  processus  de 
dimension  1  faiblement  bruite 


III 


FILTRES  APPROCHES  POUR  UN  PROBLEME  DE  FILTRAGE 

NON  LINEAIRE  D’UNE  DIFFUSION  DE  DIMENSION  2  MESUREE 
PAR  UN  PROCESSUS  DE  DIMENSION  1  FAIBLEMENT  BRUITE 

APPROXIMATE  FILTERS  FOR  A  NONLINEAR  FILTERING 
PROBLEM  OF  A  TWO-DIMENSIONAL  DIFFUSION  MESURED  BY 
A  ONE-DIMENSIONAL  PROCESS  IN  A  LOW  NOISE  CHANNEL 


Resume 

On  fait  une  etude  asymptotique  d’un  probRme  de  filtrage  non  lineaire  en  temps  con- 
tinu,  lorsque  la  variance  du  bruit  d’observation  tend  vers  0.  On  considere  un  signal 
de  dimension  2,  dont  une  seule  des  composantes  est  bruitee,  mesure  par  un  processus 
d’observation  de  dimension  1  faiblement  bruite.  Sous  certaines  hypotheses  de  detectabil- 
ite,  on  demontre  quo  I’erreur  de  filtrage  converge  vers  0  et  on  donne  une  majoration 
pour  cette  vitesse  de  convergence.  Les  demonstrations  font  intervenir  des  resultats  sur  la 
stabilite  des  sytemos  differentiels  lineaires  non  autonomes. 


Abstract 

We  study  the  asymptotic  behaviour  of  a  nonlinear  continuous  time  filtering  problem,  as 
the  variance  of  the  observation  noise  tends  to  0.  We  consider  the  case  of  a  two-dimensional 
signal  from  witch  only  one  of  the  components  is  nrlsy,  measured  by  a  one-dimensional 
observation  in  a  low  noise  chanel.  Under  some  detectability  assumptions,  we  show  that 
the  filtering  error  converges  to  0  and  we  give  an  upper  bound  for  this  converging  rate.  The 
proofs  involve  some  results  on  the  stability  of  linear  non  autonomous  differential  systems. 


1  Introduction 


Du  a  sa  vaste  appliratiou  dans  les  dotuaincs  de  I’engciiicrie,  In  luobleiiu'  du 
liltrage  d’uii  signal  aleatoiic  a  pait.ir  d’une  obsorvation  bruitde  a  fail.  I’«jbj«'l, 
il'tdiKU’s  lie  la  part  de  plusieurs  auteurs.  En  particulier,  le  cas  d’uiu?  i)bs<!rva 
lion  faiblouient  bruitee  a  cte  abondainincnt  traite,  pendant  ces  dernieriw  an- 
nbes.  On  pent  trouver,  dans  la  litterature,  des  travaux  consacres  a  la  rcdicrehe 
de  (litres  approches  asymptoticiucnient  cfficaces  quand  le  bruit  tend  vers  0. 
Darmis  ces  travaux  on  reniarque  un  premier  groupe  on  le  cas  d’un  systeme 
unidimensionnel  a3'ant  une  fonction  d’observation  injective  est  traite  (cf.  [KBS], 
[Picard-86],  [Bensoussan]).  L’etude  du  cas  multidimensionnel  parait  plus  tard 
avec  [Picard-88].  Dans  le  cas  d’une  fonction  d’observation  non  injective,  des  ef¬ 
forts  ont  ete  faits  pour  rcsoudre  certaines  classes  de  problemes.  D’abord  avec 
[Fleming-Pardoux],  qui  ont  traite  le  cas  unidimensionnel  avec  une  fonction 
d’observation  monotone  par  morceaux  et  plus  recemment  [YBS]  et  [Picard-89]. 
Dans  [YBS],  les  auteurs  etudient  formellement  le  filtrage  d’un  signal  de  dimen¬ 
sion  2  mesure  par  un  processus  d’observation  unidimensionnel.  Dans  [Picard-89], 
I’auteur  considere  le  probleme  gen4-al  du  filtrage  d’un  signal  multidimension¬ 
nel  au  vu  des  realisations  d’un  processus  d’observation  multidimensionnel  et 
il  enonce  et  demontre  des  resultats  sur  I’efficacite  du  filtre  de  Kalman  etendu 
quand  le  bruit  est  faible. 

Dans  ce  travail,  on  considere  un  systeme  pour  Icquel  la  dynamique  est  mo- 
delisee  par  un  processus  de  dilTusion  de  dimension  2,  dont  seulement  une  des 
composantes  est  bruitee,  et  I’observation  est  modelisee  par  un  processus  de 
dimension  1,  faiblement  bruite,  mesurant  I’autre  composante. 

On  considere  un  espace  de  probabilite  (fi ,  JT,  ,  P). 

On  s’intercsse  au  probleme  suivant.  On  veut  estimer,  a  I’instant  t,  un  signal 
{,Y,}  de  dimension  2  donne  par  I’equation  d’lto 


(1) 


[  (lx2{t)  -  f2{xi{t),X2(t))dt  + cr{xi{t),X2{t))dw{t) , 

avec  condition  initiale  A'q  =  (.i’i(0),.T2(0)),  au  vu  des  observations  <t{7/j,0  < 
s  <  t}  donnees  par  I’equation 


dy(t)  =  li{xi{t))dt  + edtu{t),  (2) 


ou  (u)}  et  {u’}  sont  des  processus  de  Wiener  standards  independants  a  valeurs 
dans  TZ  et  6  est  un  parametre  suppose  “petit”,  e  >  0. 


On  cherche  a  demontrer  que,  sous  certaines  hypotheses  de  regularite  sur 
les  fonctions  /  =  (/l,/2),  a  et  /i,  ce  signal  peut  etre  estime  par  un  filtre  de 
dimension  finie. 

Dans  [YBS],  les  auteurs  considerent  ce  probleme,  entre  autres.  Ils  obtiennent, 
dans  une  situation  stationnaire,  des  developpements  asymptotiques  pour  le  filtre 
optimal  et  I’erreur  quadratique  moyenne.  De  cette  etude,  ils  deduisent  un  filtre 
sous-optimal  de  dimension  finie  asymptotiquement  performant. 

Nous  allons  demontrer  que  I’erreur  de  filtrage  est  asymptotiquement  nulle 
quand  e  — ♦  0  et  qu’il  existe  un  filtre  approche  de  dimension  finie  pour  ce  pro¬ 
bleme.  Notre  point  de  depart  sera  le  filtre  approche  propose  par  [YBS].  Nous 
obtiendrons  une  estimation  de  la  vitesse  de  convergence  confirmant  le  resultat 
de  ces  auteurs.  La  methode  de  demonstration  utilisee  est  inspiree  par  le  contenu 
de  [Picard-89]. 

Ce  travail  est  organise  comme  suit.  Dans  une  premiere  etape  on  etudie  le  cas 
oil  la  fonction  /i  est  lineaire  par  rapport  a  x-i  et  la  fonction  h  est  aussi  lineaire. 
Ceci  est  fait  dans  la  section  2.  Dans  la  section  3,  on  etudie  le  cas  plus  general 
d’un  systeme  modelise  par  les  equations  (1)  et  (2).  Pour  cela  on  se  ramene,  par 
une  transformation  de  coordonnees,  au  cas  traite  dans  la  section  precedente.  On 
termine  notre  etude  avec  quelques  commentaires  sur  le  probleme. 

On  utilisera,  par  la  suite,  la  definition  suivante: 

Definition  1.1  Soil  At  une  famille  (indexee  par  e )  de  processus  mesurabies 
bornes  a.  valeurs  dans  I’espace  des  matrices  2x2.  Soil  Id  solution  (a  valeurs 
matriciclles)  de  Vequation  differentielle 

1 1 :  '» 

Soil  On  dil  que  At  est  exponentiellement  stable  avec  iaux  d’ordre 

\l\fl  s’il  exisle  des  constantes  c.  6'  >  0  Idles  que 

Vs  <  < ,  IlCt^ll  <  C e.xp{-----^-  }  , 

y/e 

oil  II .  II  represente  la  norme  ||/1||  =  sup|u|_,  |-4iz|  e<  | .  |  est  la  norme  euclidienne 
dans  TV. 

Notation  1.2  On  notera 


n(0 

r^it) 


2 


un  processus  a  valeurs  dans  'R?. 

Les  constantes  strictement  positives  independantes  de  e  intervenant  dans  les 
calculs  seront  notees  c  ou  C,  sans  expliciter  leur  valeur. 

Si  A  est  une  matrice  on  notera  A*  sa  transposee. 

Par  abus  de  notation,  on  supprimera  quelques  fois  la  dependance  en  temps 
dans  les  objets  mathematiques  ayant  un  indice  inferieur. 


2  Le  cas  particulier  d’une  derive  iineaire  par 
rapport  a  et  d’une  fonction  d’observation 
Iineaire 


On  considere  le  cas  particulier  du  systeme 
'  dxi{t)  =  fi{xi{t),X2it))dt 

<  dx2{t)  =  f2{xi{t),X2{t))dt-\- (j{Xi{t),X2{t))dw{t) 

dy(t)  =  H x\{t)  di  +  e  dw{t) , 


ayant  comme  hypotheses 
(HI)  H  >  0  et  est  une  fonction  bornee 
(H2)  a  et  A  sont  des  fonctions  C*  a  derivees  bornees  et 

:r —  =  ^12  >  U  . 

0X2 


(4) 


On  considere  le  filtre  {lut}  donne  par  I’equation  d’lto 


( 

dmi{t) 


dm2{t) 


fiimi{t),W2{t))dt  +  )l~~[dyit)  -  Hmi{t)dt] 
f2{xiit),X2{t))dt  +  ^[dy{t)  -  Hmi(t)dt], 


(5) 


avec  condition  initiale  Mq  et  oil  <t  est  une  constante  strictement  positive. 

Pour  cr  =  c*'  >  0.  le  filtre  propose  dans  (YDS)  coincide  avec  ce  filtre  quand 
on  prend  a  =  a. 

On  cherche  a  estiiner  I’erreur  de  ce  filtre,  i.e.,  pour  p  >  1,  on  s’interesse  a 
des  majorations  de  E(  II.Y,  -  Soil  {Z,}  le  processus 


l/\/?  0 

0  1 


(A'r-.U,). 


On  commence  par  ecrire  ce  proce.ssus  comme  solution  d’une  equation  dilFeren- 
tielle,  dont  la  derive  .sera  definie  par  une  certaine  matrice  .4.  De  I’etude  sur  la 
stabilite  de  cette  matrice  on  deduit,  on  utilisant  quelques  proprietes  elementaires 
du  calcul  matriciel,  Pcstimation  de  I’erreur  du  filtre. 
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2.1 


f 

1 


% 


2.1  Determination  d’une  equation  differentielle  stochas- 
tique  pour  le  processus  {Zt} 

Df^s  equations  (4)  et  (o),  on  obtient  I’equation  differentielle  stochastique 

dz,  =  °'j  'j-TW 


dt 


(  1/v'^ 

V  0 

“) 

i  j2aFios  j_  ^ 

"V  H 

a{Xt)dwt  —adxvt  y 

On  applique  la  forrnule  de  Ta3dor  a  la  fonction  f: 

/(A^.)-/(M)  =  F(6,p,)(^.-M,), 

oil  {(fj}  et  {pd  sont  des  processus  a  valeurs  dans  TV  qui  dependent  de  {A’d  et 
de  {Mi}  et 


/^K,) 

Oh,  ,  Oh:  , 

\ SI''''*/ 


On  obtient  alors  que  le  processus  {'h)  est  solution  de  I'cquation  differentielle 
stochastique 

/ 


dZi  ~  A(  Z(  dt  T 


2  af‘ \2  j. 


\ 


(6) 


(7(Xt}dw;  —  adult  j 


ou 


/ 


1  ^ 
.4,  = 


hafi^fl 


V“T 


\ 


(7) 


V  Fn{it,Ht)\/^  -  Hdjy/i  F'i2{it,tLi)  ) 
est  un  processus  a  valeurs  dans  Pespace  des  matrices  2x2. 


2.2  Etude  de  la  stability  du  processus  matriciel  {.4^} 

L’elude  de  la  stabilile  des  systemes  lineaires  non  autonomes  sera  fondamental 
dans  Tobtention  des  estimations  pour  I’erreur  de  filtrage. 

On  etablit  d’abord  le  lemme  suivant: 


■| 
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Lemme  2.1  Soit 


2a  H  —\/2crF\2H 
-y/TFT^  2Fn 


une  matrice  symetrique  et 


1 

2  V  e 


Alors  B  est  definie  positive  et 

BAt  +  A’tB  <  -XByt. 

Preuve 

Le  fait  que  B  est  une  matrice  definie  positive  est  evident.  II  suffit  de  remar- 
quer  que  traceB  >  0  et  detB  =  2a  FnH  >  0. 

Soit  Dt  =  BAt  A-  A’B  A-  XB.  On  veut  demontrer  que  Dt  est  semi-definie 
negative,  pour  tout  t. 

Un  calcul  elementaire  matriciel  nous  donne  les  expressions  de  la  trace  et  du 
determinant  de  £),: 


traceDt  =  A  a  H  an{i)  —  2  J2a  FizH  Oziit)  A  • 


ha  FnH 


-2  \j2  a  FnH  ai2(t)  +  4  Fn  «22(0  +  y  ~~~ — 


detDt  —  ‘\a  H  a\i{t) —  2\l2a  FnH  a2i{t)  A  aH 


2aFi->H 


—2  \j2a  P 12H an{t)  +4  Fn<^22{t)  +  \  — —  Fn 


—  ■“  v  2  ^  P12H  flti  (t)  A2  Fn  <^2i(0  A  2a  11  on{t)  —  v/2  a  FnH  <222(0 


aPnH 


I 


Etant  donne  I’ordre  de  grandeur  des  elements  de  la  matrice  At  et,  en  utilisant 
les  hypotheses  (Hi)  et  (H2),  on  obtient 


iraceDt  =  -{a  H  +  Fn)  +  0{\)  <  0 

Lf2  p2  1 

deWt  =  -  +0{-^)>0yt. 

£  yje 

■ 

On  pout  maintenant  demontrer  la 

Proposition  2.2  Le  processus  At  est  exponentiellement  stable  avec  taux  d’ordre 
1/V^- 

Preuve 

Soient  B  et  A  definis  par  le  lemme  2.1.  Soit  L  la  forme  quadratique  associee 
a  la  matrice  B  : 

LiX)  =  X‘BX . 

Soit  C  la  solution  de  I’equation  (3).  On  note  Ct,  la  i'"'  colonne  de  Ct»-  On  a  done 
que  est  solution  de  I’equation 

t  =  AtCs- 

.Alors  Lt  =  L{Q,)  est  une  fonction  de  Lyapunov,  puisque 

i)  L{0)  =  0,  L{X)  >  0  si  X  5^^  0 

ii)  ^IC.)  =  C:ABa,*a;b)q.. 

D’apres  le  lemme  2.1,  on  a  la  majoration 

f^(a)<-H;;Bc;.  =  -AL(a). 

Vu  que  B  est  une  matrice  bornee,  on  a 

et,  puisque  B  est  uniformement  clliptique  (B  >  cl),  on  a 


ce  qui  est  I’estimation  recherchee. 


I  1 


2.3  Estimation  de  I’erreur  de  fiitrage 

On  etablit  le  resultat  suivant: 

Th^oreme  2.3  Soit  p  un  entier,  p>l.  Si  E[l|Xo  —  A/olp'’]  <  C"  <  +oo  alors 
on  a  les  estimations 


La  preuve  de  ce  theoreme  necessite  2  lemmes  qu’on  enonce  et  demontre. 

On  definit  {QJ  comme  etant  le  processus  a  valeurs  dans  I’espace  des  matrices 
(s}'metriques  definies  positives)  solution  de  I’equation 

{  Qt  =  AtQt  +  QtA\  +  / 

I  Qo  =  gos/el- 


Lemme  2.4  La  matrict  {Qt]  est  bortiee,  d’ordre  0{y/s). 

Preuve 

Puisque  Qt  est  la  solution  de  I’equation  (8),  on  connait  sa  forme  explicite: 


Qt  —  CtQoCt  +  / 

Jo 


oil  {Ct}  est  la  solution  de  I’equation  (3). 

D’apres  la  proposition  2.2, 

vs<t, 

De  la  majorat  ion  de  la  norme  de  I’integrale  on  obtient  la  majoration  suivante 
pour  1?  norme  de  Qt  : 

IIQ.II  <  Cx/?. 


On  considere  le  processus  =  (2,  *. 

Ce  processus  est  solution  de  I’cquation  de  Riccati 

f  Vt  =  -VtAt  -  A;Vt  - 

I  =  Qo‘ . 
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Iiemme  2.5  La  matrice  Vt  est  bomee  et  elHptique  uniformement  par  rapport  a 
t.  Plus  precisement, 

C_ 

V^‘ 


Vt  >  ■—  I  et  traceVt  < 


Preuve 

Du  lemme  2.4  vient  immediatement  que 


v;> 


'  i>~i. 


iie.ii  -v/? 

De  I’equation  (9)  on  obtient,  en  utilisant  les  proprietes 


\trace{VtAi)\  <  ||i4j  ||  fracePj 

trace(Vf^)  >  -  (traceVj)^  (cf.  [Saniuk-Rhodes]) 

et  le  fait  que  <  Cjy/e,  une  inequation  differentielle  pour  la  trace  de  VJ  : 

—{traceVt)  <  -—(traceVt)  -  ^  {traceVf)^ . 

La  resolution  de  cettc  inequation  permet  d’obtenir  la  majoration: 

traceVt  < - - FTTT^  ’ 

c  y/e  +  ( 1/ (trace Vq)  —  c  ^/e)  e~^  v  ^ 

d'ou  le  lemme. 


On  est  pret  a  demontrer  le  theoreme  2.3. 

Preuve  (du  theoreme  2.3) 

Soit 

/ 


\  <r(Xt)  -a  / 


// 


D’apres  la  formuie  d’lto, 

diZ;VtZt)  =  -Z;V^Ztdt  +  trace{u;vjj,)dt  +  2Z;VtUt  ^  ^  • 
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Soit  p  >  1.  D’apres  le  lemme  de  Gronwall,  les  moments  d’ordre  p  de  ZfVtZt 
sont  finis.  On  applique  la  fornmle  d’lto  a  : 

j^iZ;VtZtY  =  -p{Z:VtZtY-^Z:V,^Ztdt  +  p{Z:VtZtY-Hrace{U;VtUt)dt 

+2p{p-i){z;v,z,Y-'\ut\Ut  +  2p{z:v,ZtY-Hz:VtU,)  (  , 

oil 

On  obtient  alors  I’expressiori 

^e[(z,-k2,)']  =  -pE{{z;v,z,)'-‘z;v,^z,\+pEi{z;v,z,r-Hrace{u;v,u,)] 
t2p{p-  i)El(z;v,z,r'  l-.p), 

D’apres  le  lemme  2.5,  on  a 

z:vYz,  >  a,Z;V,Zt 

avec  Of  =  c/^/e  et  I’in^alite  de  Cauchy- Schwartz  entraine 

krP  <  trace{U^VtL't) . 

On  obtient  I’incgalite 

jE[{Z;V,ZtY]  <  -pa,E{{Z:V,Z,Y]  +  Pi2p  -  l)E[{Z:V,ZtY''traceiU:\m]. 

Puisque  est  bornee  d’ordre  0(1),  on  a 


trace(U’VtUc)  <  Coa^. 


On  utilise  la  propriete 

V  suite  Cp  >  03  suite  Cp  :  x^~^y  <  CpX^  +  Cpi/ ,  V.t,7/  >  0 . 

On  prend  Cp  =  (1  —  c)/[(2p  -  1)  C'o],  oil  c  est  une  constante,  c  <  J. 

On  obtient  I’incgalite 


dt 


n(z:v,z,Y]  < 


~pE((Z;K2,)')+^E 


'tract[U:ViUi)Y 
.  Co  a,  ) 
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est  borne. 


7traceiU:VtUt)Y 
\  Coote  j 

En  resolvant  cette  inequation  difFerentielle  on  obtient  la  majoration 


oil,  on  le  rappelle,  E 


E[{Z:VtZtY]  <c,  +  E[iZ;VoZor]e-^P^^'^, 

oil,  vu  que  E[||.Yo  -  Mo|p'’]  est  suppose  fini,  E[(ZqVoZoY]  < 
Du  lemme  2.5,  en  utilisant  le  fait  que 


£P 


Ei(xr(t)-m,(t)Y»j  +  E[(x2(t)-m2(t)Yn  <  E[(Z;ZtY] 


£P 


on  deduit  le  theoreme. 


Remarque  2.6  En  plus  de  I’ordre  de  grandeur  de  I’erreur  de  filtrage  on  a 
obtepu  que  TefFet  de  la  condition  initiale  disparait  a  une  vitesse  exponentielle. 


3  Le  cas  general 

On  considere  le  systeme 

^  dxiit)  =  Mxi(i),X2(t))dt  +  a{xi{t),X2{t))div{t)  (10) 

.  dyit)  =  h{xi{t))dt  +  edw{t), 

ayant  comme  hypotheses 

(H3)  h  est  une  fonction  C^,  h'  est  bornee  strictement  positive,  /i"  est  bornee  et 
/i“‘  est  Lipschitzienne 

(H4)  /  est  une  fonction  /i  est  inversible  par  rapport  a  T2,  d ’inverse  Lips¬ 
chitzienne. 

On  note  g  cette  inverse 


X'J  =  ff(xi,X2)  ^  $2  =  fl(XuX2) 
et  g  la  fonction  definie  par 

g(xi,X2}  =  g  {h~^ [xx), X2I li  {h~^ [xx))^ 

(H5)  h' fx  est  une  fonction  Lipschitzienne 
(He)  la  fonction 

<^{Xx,X2)h'{Xx)~[xx.X2) 

dx2 

est  bornee 
(H7)  la  fonction 


h  (.Tl)/f(Xi,.T2)  -f  {Xl,X2)  fl{Xx.X2) 


+  '‘'(^1 )  ^  X,)  A(^„  I,)  +  fe-(x,  )  0  (x.,  X2) 


est  differentiable  et  ses  derivees  partielles  sont  bornces. 


On  va  demontrer  que,  pour  ce  systeme,  le  theoreme  2,3  enonce  dans  la  section 
precedente  est  encore  vrai  pour  p  =  1.  En  effet,  par  une  transformation  de 
coordonnees,  on  peut  se  ramener  a  un  systeme  du  type  traite  dans  la  section  2. 
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Des  estimations  de  I’erreur  pour  ce  systeme  transforme,  on  va  pouvoir  deduire 
des  estimations  analogues  pour  le  systeme  initial. 

On  considere  la  transformation 

f  xi{t)  =  h{xi{t)) 

\  12(0  =  h'{Xi{t))  fi{Xi{t),X2(t))  . 

Le  nouveau  systeme  est,  d’apres  la  formula  d’lto, 

dxi{t)  =  X2{t)dt 

'  dx2{t)  =  f{xi{t),X2{i))dt-\-a[xi{t),X2{t))dw{t)  (11) 

dy{t)  —  Xi{t)  dt  +  £ dib{t) , 

ou 


f{xi,X2)  =  li"(h~\xi))f^{h  ‘(xi),5(a:i,X2)) 

+h'[lr\xi))^^{h-^xi),g(xux2))fi{h-\xi),g{xi,x2)) 

+h'{lr^{xi))^(li-^{xi),g{xuX2))f2{h-'^{xi),g{xuX2)) 


.  cr‘^ih'Hxi),g{xi,X2)) 
+ 


kVr\xx)) 


dxl 


(/r*(xi),^(5i,x2)) 


et 

a'(.f,,.f2)  =  (r(lr^{xi),g{xi,X2))h'{}r^{xi))~{h~‘^(xi),g{xi,X2)). 

0X2 

On  est  dans  les  conditions  de  la  section  precedente.  Pour  le  systeme  (11),  on 
considere  le  filtre 


dm]{t)  =  m2{t)  dt  +  J  —  [dy{t)  —  mi{t)  dt] 

^  "  (12) 

dw2it)  =  f{m]{t),in2{t))dt-\- ^[dy{t) -mi{t)dt], 


avec  condition  initiale  jVIq  et  ou  a  est  une  constante  strictement  positive. 
Pour  le  systeme  initial  (systeme  (10)),  on  considere  naturellement  le  filtre 


mi(i)  =  h-^xi(t)) 

m2{t)  =  g{fiii{t),m2{t)) . 


(13) 


On  a  le 
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Theoreme  3.1  Si  E[  |jA'o  -  A/olP]  <  C"  <  +oo  alors  on  a  les  estimations: 

E[(.'ci(<)-mi(0)^]  <  + 

Elixiit)  -  m2{t)f]  <  +  Ce^l‘^  . 


Preuve 

Le  systeme  transforme  (11)  verifie  les  hypotheses  enoncees  dans  la  section 
precedente  et  EfllvVo  -  A/o|P]  <  C,  car  la  fonction  h'fi  est  supposee  Lipschitzi- 
enne.  On  pent  done  appliquer  le  theoreme  2.3: 

E((ii(t)  -  m,(t))2j  <  + 

E[(.T2(<)  -  m2(0)')  <  +  05^/2 

£ 

E[(x,(0-mi(t))-*j  <  +  Ce^  . 

D’autre  part,  puisqu’on  suppose  que  les  fonctions  h~^  et  g  sent  Lipschitziennes 
et  h'  est  born^,  on  a  les  majorations 

E{(x, -m,)']  <  CE[(x, -m,)') 

E((T2-m2)^|  <  CE((.f2  -  m2)^)  +  CE[i2(3:i  -  mi)^] . 

Vu  que  Efxj]  <  C  <  +oo,  I’inegalite  de  Holder  permet  d’obtenir  le  theoreme. 


Remarque  3.2  On  pent  penser  a  appliquer  a  ce  probleme  le  meme  genre  de  de¬ 
marche  que  celle  utilisee  dans  la  section  2.  On  obtient  alors  un  resultat  analogue 
au  theoreme  3.1,  mais  seulement  si  on  fait  une  restriction  sur  la  variation  des 
derivees  des  fonctions  li  et  fi  :  on  cxige  que  les  valeurs  prises  par  ces  fonctions 
restent  dans  un  certain  voisinage. 

Plus  precisement,  sous  les  hypotheses 

(H8)  h  est  une  fonction  C‘,  h'  est  bornee  slrictement  positive  et  a  est  une 
fonction  bornee  strictement  positive 

(H9)  /  est  une  fonction  (7*  a  derivees  bornees  et 

^  J 

^(.T,,l2)  >  0,  V{X,,X2). 
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Si  on  note: 


9/i 


<f>x 


dx 


-(*1,2:2) 


/l'(Xi) 

h  =  infj;/i'(x), 


H  =  sup^  h'{x) , 


fu  =  mf(:„,^2)^(2-’l,*2),  I'n  =  SUP(^,,.J)  — (Xi,X2) 

et  on  considere  le  filtre  approche 


2  a 


dmi{i)  =  /i(m,(0,m2(0)d/ +  J — (pM  [dy{t)  -  h{mi{t))dt] 


drriiit)  =  f2{mi(t),m2it))dt  + -[dii{t)  ~  h{mi{t))dt], 
oil  a  est  une  constante,  o'  >  0,  on  pent  dcmontrer  le  resultat  suivant: 


SI 


^  HF\2  ^  (\/l+4(l  +  v^)-l)2 

-  hfu  -  4 


(‘) 


alors 


(14) 


E((x,(0-ni,(0)'1  <  Cpe-^P^^  +  Cpe3p/2 

E((.r2{f)-m2(0)"1  <  + 


La  metliode  qui  a  pcnnie  d’obtenir  !a  condition  (14)  n’est  pas  exhaustive. 
On  pout  facilement  trouver  d’autres  conditions  du  metnc  genre  pour  Icsquelles 
le  resultat  est  encore  vrai. 


{Jl  +  4(1  +  \/2)  -  1)- 
- - : - ss  1. 


282 
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4  Commentaires 


On  a  etudie  I’erreur  de  filtrage  pour  le  cas  d’un  signal  de  dimension  2  dont  une 
seule  des  composantes  est  bruitee  et  dont  seule  I’autre  composante,  faiblement 
bruitee,  est  observee.  On  a  obtenu  des  estimations  asymptotiques  pour  cette 
erreur  qui  nous  donnent  I’ordre  de  grandeur  quand  I’intensite  du  bruit  de  mesure 
£  tend  vers  zero.  Une  etude  asymptotique  de  I’erreur  du  filtre  approche  par 
rapport  au  filtre  optimal  est  aussi  envisageable  (cf.  [Picard-89]). 

On  pent  envisager  le  meme  type  d’etude  pour  le  cas  ou  c’est  la  composante 
bruitee  du  signal  qui  est  observee.  Pour  ce  cas  la,  dans  [YDS],  les  auteurs  pro- 
posent  un  filtre  de  dimension  2  obtenu  formellement  d’apres  un  developpement 
asymptotique  du  filtre  optimal. 
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